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FOREWORD 


This  report  presents  the  conclusions,  recommendations,  and  typical  results 
available  from  the  efforts  to  date  on  the  project  entitled  "Aerodynamic 
Analysis  of  the  Self- Suspended  Flare.  "  This  project  was  conducted  by  the 
Government  and  Aeronautical  Products  Division  of  Honeywell  Inc.  for  the 
Research  and  Development  Department  of  the  Naval  Ammunition  Depot 
(NAD),  Crane,  Indiana,  under  Contract  N00164-69-C-0662.  The  effort 
described  was  essentially  conducted  during  the  period  from  Tune  1969 
through  June  1970,  with  intermittent  computer,  documentation,  and  present¬ 
ation  efforts  in  the  interim. 

The  cognizant  Project  Engineer  for  NAD/Crane  was  Mr.  J.  J.  Riester.  The 
principal  investigator  for  Honeywell  was  Mr.  G.  D.  Stilley,  with  support 
from  other  aerodynamics,  computer  programming,  wind  tunnel,  and  flight 
test  photography  personnel,  as  required.  The  studies  were  also  supported 
by  in-house  facilities  consisting  of  a  subsonic  wind  tunnel,  a  CDC  6600  com¬ 
puter  terminal,  a  Honeywell  Computer  Network  teletype  computer  terminal, 
several  cameras,  and  by  a  semi-automatic  motion  analyzer  for  film  data. 

The  efforts  were  also  supported  on  a  mutual  team  work  basis  by  free-flight 
test  facilities  at  Hurricane  Mesa,  which  were  the  responsibility  of  a  separate 
contractor,  Denver  Research  Institute,  under  direct  contract  with  NAD/Crane. 

Some  informal  consultive  support  was  obtained  from  Dr.  J.  D.  Nicolaides 
and  Dr.  C.  W.  Ingram  of  the  Aerospace  Engineering  Department  at  Notre 
Dame  University  and  from  Dr.  E.  E.  Covert  from  the  magnetic  support 
wind  tunnel  facility  at  the  Massachusetts  Institute  of  Technology  (MIT). 

This  technical  report  has  been  reviewed  and  is  approved. 
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It  is  recognized  that  the  name  "Frisbee"  is  a  copyrighted  trademark  of  the 
Whamo  Corporation.  Therefore,  the  more  precise  and  legally  correct  way 
to  describe  the  domed,  cupped  configurations  studied  herein  would  be  with 
the  terms  "frisbee-like  configuration,  "  etc.  However,  in  the  interest  of 
brevity  and  clarity,  especially  on  figures,  this  report  will  define  the 
word  "frisbee"  to  mean  "frisbee-like  configuration.  "  Therefore,  "solid 
frisbee"  would  refer  to  a  frisbee-like  configuration  whose  top  is  domed 
similarly  to  the  "Frisbee"  bui  whose  bottom  surface  is  filled  in  flush,  etc. 
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I.  INTRODUCTION 


A.  OBJECTIVE  OF  THE  SELF-SUSPENDED  FLARE  PROJECT 

The  objective  of  the  Self-Suspended  Flare  project  conducted  by  NA D/Crane 
was  to  develop  an  air  launched  flare  concept  which  would  utilize  gyroscopic 
stabilization  of  the  case  of  a  disk  shaped  flare  to  retard  the  descent  of  the 
flare  sufficiently  to  provide  intermediate-term  illumination,  [without  the 
parasitic  weight  and  volume  of  conventional  parachute  flares]  (see  Figure  1). 

It  was  not  intended  to  provide  long-term  battlefield  illumination  and  the 
associated  very  slow  descent.  Preliminary  analysis  indicated  that  reasonable 
descent  velocities  could  be  achieved  for  reasonable  times  and  composition 
weights  if  such  a  flare  could  be  kept  oriented  and  if  the  proper  combustion 
characteristics  could  be  obtained  in  a  spinning  body.  Preliminary  free- 
flight  and  pyrotechnic  laboratory  experiments  indicated  that  this  could  be 
done. 

The  analysis  of  the  aerodynamic  feasibility  of  this  concept  was  assigned  to 
the  Government  and  Aeronautical  Products  Division  (formerly  Ordnance 
Division)  of  Honeywell  Inc. ,  which  had  the  required  capability  and  experience 
in  flight  and  fuzing  of  unconventional  ordnance.  The  pyrotechnic  and  case 
study  and  development  was  primarily  conducted  in-house  by  the  Research  and 
Development  Department  of  NAD/Crane,  whose  primary  mission  was  pyro¬ 
technic  development.  Design,  fabrication,  and  operation  of  the  launcher  was 
assigned  to  the  Denver  Research  Institute,  which  had  been  providing  con¬ 
tinuing  research  support  to  NAD/Crane. 

B.  OBJECTIVES  OF  THE  AERODYNAMIC  ANALYSIS  OF  THE  SELF- 
SUSPENDED  FLARE 

The  overall  objectives  of  aerodynamic  analysis  of  the  self-suspended  flare 
were  as  follows: 
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PATH  OF  FLIGHT 


Figure  1.  Self- Suspended  Flare  Concept 


.  Establish  by  analytical  means  the  aerpdynamic  feasibility  of  the  self- 
suspended  flare  concept. 

.  Define  the  optimum  aerodynamic  configuration. 

.  Obtain  aerodynamic  data  for  use  in  further  development. 

C.  INTERIM  GOALS  FOR  THE  AERODYNAMIC  ANALYSIS  OF  THE 
SELF-SUSPENDED  FLARE 

Because  of  the  complexity  of  the  task  and  the  magnitude  of  the  effort  that 
could  have  been  required  to  guarantee  achievement  of  these  objectives, 
compared  to  the  time,  funds,  and  flare  combustion  data  available,  interim 
goals  were  identified  which  would  both  have  a  higher  probability  of  success 
and  provide  a  firmer  basis  for  planning  of  more  extensive  follow-on  efforts. 

D.  SCOPE  OF  EFFORTS,  AERODYNAMIC  ANALYSIS  OF  THE  SELF- 
SUSPENDED  FLARE  (CONTRACT  N001 64-69-C-0662) 

The  planned  scope  of  the  efforts  directed  at  these  interim  goals,  which  were 
considered  compatible  with  the  overall  Self-Suspended  Flare  project 
development  status  and  funding  level,  encompassed  the  following  major 
tasks,  in  addition  to  the  necessary  coordination  and  documentation: 

* 

.  Trajectory  simulation. 

.  Nonspinning  wind  tunnel  tests  of  three  basic  configurations  (each  to 
be  solid  and  hollow). 

.  Establishment  of  aerodynamic  criteria. 

As  the  Hurricane  Mesa  tests  materialized,  the  coordination/liaison  efforts 
became  more  significant  than  normal. 

On  the  basis  of  the  preliminary  results  obtained  and  the  concurrent  scheduling 
of  relatively  independent  Hurricane  Mesa  tests,  these  tasks  were  informally 
diverted,  in  concurrence  with  the  NAD/Crane  Project  Engineer,  to  include 
full  scale  participation  in  the  Hurricane  Mesa  tests  of  November  1969. 
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This  was  expected  to  increase  the  yield  from  both  projects.  The  objectives 
of  this  joint  test  effort  were  to  maximize  the  data  obtainable  from  the  flight 
tests  and  if  possible  to  obtain  aerodynamic  or  dynamic  data  not  readily 
obtained  from  the  existing  wind  tunnel  models  and  simulation  and/or  to 
validate  the  existing  data. 

Subsequent  to  the  Hurricane  Mesa  tests,  the  scope  of  Contract  N00164-69-C-0662 
was  extended  to  extract  and  utilize  the  test  data  and  to  closely  examine  the 
aerodynamic  phenonema  involved.  These  additional  tasks  were  as  follows: 

.  Analysis  of  Hurricane  Mesa  test  data. 

.  Incorporation  of  Hurricane  Mesa  test  data  in  the  flare  analysis. 

.  Spinning  wind  tunnel  model  test  of  one  configuration. 

While  the  long  range  objectives  have  not  yet  been  achieved,  the  more  reason¬ 
able  interim  goals  have  essentially  been  met  in  terms  of  criteria  and  a  data 
base  for  further  work. 

This  final  summary  report  for  Contract  N00164-69-C-0062  is  arranged  in 
the  format  of  Summary  (Section  II),  Conclusions  and  Recommendations 
(Section  III),  Technical  Discussion  (Section  IV),  Bibliography  (Section  V), 
and  Appendices.  The  Technical  Discussion  itself  consists  of  a  narrative 
semichronological  summary  of  the  Aerodynamic  Analysis  of  the  Self- 
Suspended  Flare  and  the  rationale  applied  and  typical  results  of  these 
efforts  in  the  areas  of  trajectories,  analytical  criteria,  configuration 
aerodynamic  data,  flight  test  results,  and  math  models. 


II.  SUMMARY 


A.  THEORETICAL  DEVELOPMENT 


The  following  are  the  accomplishments  of  the  Aerodynamic  Analysis  of  the 
Self- Suspended  Flare  in  theoretical  development: 


.  A  relatively  versatile  trajectory  simulation  has  been  developed  and 
typical  results  are  presented  herein. 

.  Results  correlate  qualitatively  with  available  free-flight  tests  and 
analytical  solutions. 

.  Preliminary  analytical  criteria  relating  aerodynamic,  inertial,  and 
delivery  characteristics  have  been  derived  and  correlated  with  the 
well-known  flight  of  toy  frisbees  and  the  flight  of  clay  pigeons  on 
the  range  or  tossed  from  Hurricane  Mesa. 

.  Key  features  for  those  devices  are  the  hollow  domed  shape,  light 
weight,  and  delivery  orientation. 

.  Cursory  applications  of  the  trajectory  model  to  other  ordnance  delivery 
applications  show  promise  for  these  applications. 

.  The  trajectory  simulation  is  available  for  study  of  perfornance  of  ' 
inert  disks  or  burning  flares  (linear  burning,  constant  thrust). 

.  The  simulation  permits  the  assumption  of  either  steady  gyroscopic 
precession  under  the  aerodynamic  moments,  or  gyroscopic  rigidity, 
depending  on  the  purpose  of  the  study  and  the  available  data. 

.  For  the  low  burning  pressures  originally  anticipated  for  the  flare,  the 
analysis  indicated  negligible  effect  on  the  aerodynamic  forces. 


B.  TESTING 


The  following  are  the  accomplishments  of  the  Aerodynamic  Analysis  of  the 
Self- Suspended  Flare  in  testing: 


Preliminary  wind  tunnel  tests  have  been  conducted  on  typical  solid 
and  hollow  candidate  shapes  and  trends  observed. 

Configuration  parameters  warranting  further  wind  tunnel  test  with 
existing  techniques  have  been  identified. 


-5- 


Limitations  of  these  techniques  and  the  required  improvements  or 
alternates  have  been  identified,  considered,  and  deferred  for  future 
development  program(s). 

Results  of  preliminary  free-flight  tests  (conducted  by  another  contractor) 
have  correlated  qualitatively  with  theoretical  results,  in  addition  to 
indicating  areas  for  pyrotechnic  research. 

Free-flight  tests  to  date  have  not  yet  adequately  simulated  the  air- 
launch  environment. 

Theoretical  trajectory  extrapolations  to  air-launch  conditions  indicate 
that  either  the  optimum  design  has  not  yet  been  tested  or  identified  or 
that  the  theory /model /data  require  additional  free-flight  validation. 
Configurations  which  appeared  to  manifest  acceptable  precession  and 
negligible  trajectory  displacement  under  Hurricane  Mesa  test  conditions 
appear  to  precess  excessively  for  simulations  of  higher  speeds  or 
higher  drop  altitudes. 

Wind  tunnel  tests  indicate  that  domed  flare  configurations  generate 
some  of  the  key  aerodynamic  characteristics  of  clay  pigeons  or  toy 
frisbees. 

The  experience  gained  to  date  in  instrumenting  and  analyzing  the 
Hurricane  Mesa  tests  indicates  that  instrumenting  formal  airdrops 
will  be  difficult.  Detailed  dynamic  or  aerodynamic  information  will 
most  likely  not  be  obtainable.  However,  basic  trajectory,  performance, 
and  possibly  lift/drag  parameters  will  be  available. 
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III.  CONCLUSIONS  AND  RECOMMENDATIONS 


A.  CONCLUSIONS 


The  following  are  the  conclusions  of  the  Aerodynamic  Analysis  of  the  Self- 
Suspended  Flare: 


Interim  goals  of  establishment  of  trajectory  models  and  aerodynamic 
requirements  have  been  met  and  validated  on  typical  configurations. 

Adaptation  to  denser  ordnance  items  could  take  either  the  course  of 
configuration  and  delivery  as  scaled  analogies  to  frisbee  tosses  with  tilted 
domed,  cupped  configurations  or  by  neutralizing  the  aerodynamic 
characteristics  with  relatively  thick  cylinders,  or  by  a  combination 
thereof. 

Math  models  are  available  to  perform  preliminary  application  require¬ 
ment  studies  in  parallel  with  the  completion  of  feasibility /configuration 
efforts. 

As  may  well  have  been  expected  within  the  scope  of  the  program  to 
date,  insufficient  data  are  yet  available  to  conclusively  establish 
concept  feasibility  or  infeasibility  and  the  optimum  configuration 
required. 

Such  feasibility  and  configuration  establishment  will  entail  additional 
effort  balanced  between  analysis,  computer  simulation,  wind  tunnel 
testing,  pyrotechnic  development,  and  free-flight  testing  of  inert 
and  live  flares. 

In  view  of  some  of  the  results  and  limitations  of  the  aerodynamic 
studies  to  date,  informal  flight  tests  of  typical  configurations,  at  the 
earliest  opportunity,  would  be  extremely  valuable  to  provide  the 
following: 

—  Interim  concept  validation/configuration  selection  prior  to 
committment  to  live  flare  demonstration. 

—  Validation  of  available  math  models  and/or  aerodynamic  data. 

—  Identification  of  technical  problems  for  both  flare  development  and 
flight  test  instrumentation. 

—  Adequate  emphasis  in  forthcoming  studies  between  applications, 
flight  dynamics,  and  pyrotechnic  development. 
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B.  RECOMMENDATIONS 


The  following  are  the  recommendations  for  future  development: 


.  From  the  available  data,  the  most  promising  configurations  appear 
to  be  the  3:1  to  2:1  diameter-to-thickness  ratio  right-circular  cylinder 
and  the  relatively  thick  domed /cupped  clay  pigeon-like  shape,  (see 
Figures  2,  and  3)  which  also  show  some  possible  alternates. 

.  Because  of  the  interaction  between  aerodynamics  and  launch  conditions, 
the  choice  of  configuration  and/or  launcher  design  must  take  into 
account  the  importance  of  flare  orientation  at  launch  to  provide  desired 
steady-state  glide  characteristics. 

.  The  next  step  in  flare  development  should  be  preliminary  flight  tests 
using  available  models  (large  6:1  and  4:1  cylinders)  and  the  most 
promising  and  easily  made  candidates  such  as  3:1  or  2:1  cylinders  - 
small  heavy  cylindrical  or  clay  pigeon  shapes  (perhaps  ballasted 
actual  or  plastic  clay  pigeons), using  the  present  launcher  for  the  large 
models  or  a  modified  clay  pigeon  launcher  for  the  small  models. 

.  The  near-term  subsequent  steps  in  aerodynamic  development  should 
consist  of  the  following: 

-  Systematic  iteration  of  this  family  of  configurations  with  the  avail¬ 
able  wind  tunnel  techniques  and  math  models  to  separate  effects 

of  corner  radii,  thickness,  and  cavity. 

-  Extension  of  the  aerodynamic  criteria  and  math  models  as  required. 

-  Instrumented  flight  tests  for  gross  performance  data  whose  nature 
would  be  partially  contingent  on  the  results  of  the  informal  flight 
tests. 

-  The  development  of  improved  wind-tunnel  techniques  should  probably 
be  deferred  because  of  the  complexity,  time,  and  cost  involved  and 
the  unknown  effects,  although  a  Magnus  moment  balance  measure¬ 
ment  should  be  obtained  as  soon  as  economically  possible. 
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Figure  3.  Possible  Alternate  Aerodynamic  Configu*  ations, 
Telf- Suspended  Flare 
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IV.  TECHNICAL  DISCUSSION 


This  section  presents  a  description  of  the  technical  problems  and/or  the 
efforts  directed  at  their  solution,  followed  by  presentation  of  the  typical 
results  of  these  efforts. 


A.  SUMMARY  OF  TECHNICAL  CONSIDERATIONS  AND  EFFORTS 


The  efforts  expended  to  date  on  the  Aerodynamic  Analysis  of  the  Self- Suspended 
Flare  have  encompassed  overlapping  tasks  of  dynamic  analysis,  aerodynamic 
analysis,  wind  tunnel  testing,  and  free-flight  testing.  The  somewhat  arbitrary 
distinctions  are  that  the  dynamic  analyses  meant  trajectory  and  stability  and 
precession  analyses  utilizing  the  physical  and  aerodynamic  coefficients  for  a 
given  configuration,  or  generation  of  analytical  solutions.  The  aerodynamic 
analyses  meant  determining  the  aerodynamic  coefficients  of  a  given  configura¬ 
tion  from  estimates  based  on  theory,  available  data,  or  new  wind  tunnel  data, 
and  synthesizing  promising  configurations  for  further  study.  Wind  tunnel 
testing  included  the  actual  fabrication  of  models  and  wind  tunnel  testing  and 
data  reduction  for  static  and  dynamic  aerodynamic  coefficients.  Since  the 
contractor  was  not  responsible  for  conducting  free-flight  tests,  the  flight  test 
efforts  included  observation  and  photography  of  tests,  recommendation  of 
instrumentation  and  test  conditions,  and  analysis  of  photographic  data  for 
trajectory  parameters. 

These  efforts  were  directed  toward  the  ultimate  objectives  of  establishing 
the  feasibility  of  the  concept  and  specifying  the  optimum  flare  configuration 
and  delivery  conditions.  The  relatively  austere  allowable  scope  of  the  effort, 
especially  as  originally  planned  independent  of  flight  tests,  did  not  provide  a 
guarantee  that  these  objectives  could  be  met  in  this  phase  of  Self-Suspended 
Flare  development.  It  was  hoped,  however,  that  the  interim  goal  of  estab¬ 
lishing  aerodynamic  criteria  could  be  achieved. 

These  aerodynamic  criteria  would  be  intended  to  provide  satisfactory  flight 
behavior  if  and  when  met  by  some  available  future  configuration/delivery 
combination.  Of  course,  if  it  were  then  immediately  obvious  that  such 
criteria  could  never  be  met,  then  infeasibility  would  be  indicated.  The  goals 
of  criteria  establishment  has  been  achieved  and  validated  for  familiar  flying 
objects,  but  the  available  data  to  not  yet  permit  conclusions  about  the  ability 
to  meet  them  in  practice  with  the  Self-Suspended  Flare. 


-12- 


The  considerations  and  efforts  will  be  described  in  a  loosely  chronological 
narrative  because  of  the  considerable  overlap,  shifts  in  emphasis,  and  cross 
interpretation  involved. 

1.  Problem  Identification/ Definition 

The  right-hand  rule  of  gyroscopic  precession  means  that  if  a  torque  is  applied 
to  a  symmetrical  spinning  object,  the  object  will  rotate  steadily  about  an  axis 
at  right  angles  to  the  axis  about  which  the  torque  is  applied.  This  rotation 
will  be  at  a  rate  proportional  to  the  applied  torque  and  inversely  proportional 
to  the  angular  momentum  about  the  spin  axis.  This  steady-state  rotation  is 
referred  to  as  gyroscopic  precession,  and  the  steady-state  rotations  of  the 
flare  will  thus  be  referred  to  as  precession  in  this  report. 

This  right-hand  rule  ignores  the  transient  accelerations  of  the  spinning  object 
which  usually  have  a  negligible  effect  on  the  amplitudes  of  the  motion  for  most 
rapidly  spinning  bodies.  The  high  frequency  transients  are  quite  frequently 
referred  to  as  nutations. 

The  gyroscopic  precession  of  a  top  involves  coning  of  the  spin  axis  about  the 
vertical  because  the  forcing  moment  due  to  gravity  is  always  in  the  plane 
which  contains  the  vertical  and  the  spin  axis  and  which  rotates  with  the  spin 
axis.  Therefore,  instead  of  tipping  over  at  right  angles  to  the  torque,  the 
top  keeps  chasing  the  gravity  torque.  This  plane  could  be  called  the  gravity 
plane  (see  Figure  4). 

The  gyroscopic  precession  of  a  disk  flying  roughly  edge-on  involves  rolling 
about  the  velocity  vector  because  the  aerodynamic  overturning  moment  due 
to  an  angle  of  attack  with  respect  to  the  velocity  vector  is  in  the  plane  which 
contains  the  velocity  vector  and  the  spin  axis.  It  rotates  with  the  spin  axis 
as  it  precedes  about  the  velocity  vector  and  maintains  the  angle  of  attack. 
This  plane  will  be  called  the  aerodynamic  plane  (see  Figure  4). 
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Figure  4.  Gyroscopic  Precession  Principles 
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Because  the  disk  in  free-flight  is  "suspended"  at  its  center  of  gravity,  it 
does  not  have  a  corresponding  moment  due  to  gravity  which  would  directly 
cause  the  more  desirable  coning  about  the  vertical. 

Both  of  these  effects  can  be  demonstrated  with  toy  or  laboratory  gyroscopes 
which  are  suspended  by  the  end  of  the  spin  axis  or  at  the  center  of  gravity, 
respectively.  It  was  recognized  that  in  addition  to  any  initial  inclination  to 
the  velocity  vector,  the  curvature  of  the  path  under  gravity  could  impose  an 
additional  angle  of  attack  leading  to  roll  precession. 

Somehow,  toy  frisbees  and  clay  pigeons  have  been  able  to  bypass  this  appli¬ 
cation  of  the  right  hand  rule  of  gyroscopics.  This  was  apparently  attributable 
to  the  cupped  shape  of  these  objects  or  to  some  dihedral-like  response  to 
sideslip  caused  bytilting.  Therefore  the  problem  was  to  determine  if  this 
behavior  could  be  identified  and  extrapolated  to  an  initially  solid  body  launched 
at  high  speed  for  a  finite  flight  time. 

2.  Technical  Approach 

The  approach  to  dynamic  analysis  of  the  flight  of  the  flare  is  the  standard  one 
in  flight  dynamics;  separation  of  the  slowly  varying  parameters  which  define 
the  trajectory  from  the  higher  frequency  motions  or  oscillations  about  the 
trajectory.  The  trajectories  would  then  be  calculated  by  digital  simulation, 
while  the  motions  of  the  body  with  respect  to  the  trajectory  would  be  treated 
by  stability/ response  analytical  solutions  at  selected  points  on  the  trajectory. 
If  insufficient  data  were  available  to  quantitatively  predict  the  stability/ 
response  of  a  specific  case,  the  analysis  would  take  the  form  of  stability/ 
response  criteria  to  be  met  by  the  missing  data  for  satisfactory  behavior. 

The  two  separated  motions  could  later  be  combined  in  a  6-degree-of-freedom 
simulation  if  desired. 

The  separation  was  necessary  and  reasonable  because  the  spin  rate  introduces 
nutational  motions,  normally  excited  by  impulsive  inputs,  with  frequencies 
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on  the  order  of  twice  the  spin  rate.  This  frequency  would  call  for  fantasti¬ 
cally  large  amounts  of  computer  time  for  the  flight  times  involved,  and  the 
simulation  would  fail  if  complete  aerodynamic  data  for  a  stable  body  were 
not  provided.  By  the  same  token,  the  high  frequency  motions  would  normally 
have  such  low  amplitudes  for  a  stable  body  that  they  would  have  negligible 
effect  on  the  trajectory.  T'.e  latter  observation  was  validated  by  the  stability 
analysis. 

These  dynamic  analyses  would  draw  on  aerodynamic  data  as  available  from 
literature,  theoretical  estimates,  and  wind  tunnel  tests  of  representative 
configurations.  To  serve  as  a  starting  point,  the  configurations  selected  for 
initial  study  were  a  16  inch  diameter,  2  inch  thick  right  circular  cylinder 
yielding  an  8:  1  diameter  to  thickness  ratio,  and  solid  and  hollow  versions  of 
the  familiar  clay  pigeon  and  the  toy  frisbee. 

This  cylinder  at  that  time  appeared  to  be  the  prime  candidate  for  tactical 
deploy  ment.  The  frisbee  and  clay  pigeon  presented  possible  variations  in 
corner  radius  and  thickness  from  the  cylinder,  which,  of  course,  was 
expected  to  be  the  preferred  shape  for  pyrotechnic  design.  Since  they  were 
part  of  the  inspiration  of  the  concept  and  their  flight  characteristics  were 
supposedly  well-known,  the  results  obtained  were  expected  to  assist  in 
validation  of  the  techniques  developed  for  the  Aerodynamic  Analysis  of  the 
Self- Suspended  Flare. 

Initially,  the  approach  did  not  include  formal  flight  testing  because  launcher 
design  and  operation  was  the  responsibility  of  another  contractor  and  the 
schedule  was  not  necessarily  compatible.  As  the  study  proceeded,  this 
launcher  became  available  and  the  level  of  the  aerodynamic  contractor's 
participation  in  the  testing  increased.  This  division  of  responsibility  limited 
the  in-house  "flight  test"  to  observations  of  informal  clay  pigeon  and  toy 
frisbee  flights.  In  retrospect,  a  parallel  program  of  more  formal  in-house 
small  scale  flight  tests  would  have  been  valuable  if  the  scope  and  funds  had 
permitted. 
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3 .  Aerodynamic  Force,  Moment,  and  Axes  Definitions 


The  aerodynamic  coefficients  presented  in  this  report  are  defined  in  terms  of 
the  aerodynamic  axis  systems  and  Figure  5  (and  nomenclature  of  Appendix  B). 

The  aerodynamic  axes  have  their  origin  at  the  center  of  gravity  of  the  disk 
and  are  expressed  in  terms  of  the  aerodynamic  plane  in  a  manner  somewhat 
analogous  to  symmetrical  projectiles,  except  that  the  disk  is  considered  at 
rest  or  at  zero  orientation  and  angle  of  attack  when  edge -on.  Because  of 
rotational  symmetry,  this  aerodynamic  plane  is  independent  of  arbitrary 
transverse  axes  of  the  body.  The  aerodynamic  plane  is  defined  to  be  the 
plane  containing  the  positive  spin  axis  (designated  Z^)  and  the  velocity  vector. 
The  aerodynamic  roll  axis,  (X^)  is  perpendicular  to  the  spin  axis  in  this 
aerodynamic  plane.  The  aerodynamic  pitch  axis,  (Y^)  then  completes  the 
right-hand  rule  triad  of  axes.  and  may  be  considered  to  define  the 
disk  plane. 

The  aerodynamic  plane  corresponds  to  the  pitch  plane  as  the  disk  is  tilted  in 
the  wind  tunnel  for  measurements  of  forces  and  moments.  The  angle  of 
attack,  1^,  is  defined  as  the  angle  between  the  velocity  vector  and  the  disk 
plane  (or  axis)  measured  in  the  aerodynamic  plane  perpendicular  to  the 
disk  plane.  The  aerodynamic  forces  expressed  in  standard  body  axis  form 
for  projectiles  or  missiles  are  the  axial  force  A,  parallel  and  opposite  to 
the  X^  axis,  while  the  normal  force,  N,  is  normal  to  the  disk  plane  and 
parallel  and  opposite  to  the  positive  spin  axis.  The  normal  force  may  be 
considered  to  act  at  the  center  of  pressure,  which  is  X£p  forward  of  the 
center  of  gravity  along  the  X^  axis. 

The  velocity  axis  forces,  lift  (L)  and  drag  (D)  are,  as  usual,  perpendicular 
and  parallel  respectively  to  the  velocity  vector  in  the  aerodynamic  plane. 

L  and  D  or  N  and  A  can  be  resolved  back  and  forth  as  most  convenient,  but 
the  present  day  integral  wind  tunnel  balances  are  parallel  and  perpendicular 
to  N  and  A  and  give  them  more  directly. 
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Figure  5.  Aerodynamic  Axes,  Forces,  and  Moments 
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If  there  were  a  side  force  (Y),  which  would  have  to  be  due  to  spin  for  these 
disks  with  axes  of  symmetry,  it  would  be  defined  as  parallel  to  the  aero¬ 
dynamic  pitch  axis,  Y^.  This  force  has  been  ignored  so  far  but  could 
become  significant.  It  would  be  expressed  in  the  Magnus  force  coefficient 
form  (CY  )  with  the  coefficient  a  function  of  the  angle  of  attack. 

The  pitching  moment  (M)  is  considered  to  be  a  torque /couple  in  the  aero¬ 
dynamic  plane  or  a  torque  vector  along  the  positive  Y^  axis,  positive  over¬ 
turning.  In  addition  to  the  static  moment,  due  to  the  angle  of  attack 

of  the  disk,  the  pitching  moment  description  includes  a  damping  moment 
(M(Q)  =  Mq  Q^)  due  to  an  angular  rate  about  the  Y^  axis.  This  pitch  damping 
moment  is  also  a  function  of  the  angle  of  attack.  The  effect  of  the  spin  on 
the  pitching  moment  is  considered  negligible. 

The  standard  aerodynamics /ballistics  approach  to  expressing'the  aerodynamic 
forces  and  moments  in  terms  of  shape  dependent  dimensionless  coefficients 
which  are  independent  of  size  and  mass  and  velocity  (except  as  dependent  on 
Mach  number  or  Reynolds  number)  has  been  followed  here.  The  wind  tunnel 
data  are  all  presented  in  this  form.  The  reference  area  for  this  process  has 
been  taken  as  the  disk  "cross-sectional"  area  (S  =  n  r  )  in  square  feet  and 
the  reference  length  as  the  diameter  (d)  in  feet.  Defining  the  dynamic  pres¬ 
sure  as  q,  the  force  coefficients, of  the  form  Cp  and  moment  coefficients, 
of  the  form  (^,  are  related  to  their  respective  forc  es  (F)  and  moments 
(M)  by  the  following  typical  expressions: 

CF  =  F/qS 
c  M  =  M/iSd 
_  2 

where  q  =  1/2  oV  (pounds  per  square  foot) 

o  =  atmospheric  density  (slugs  per  cubic  foot) 
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These  coefficients  which  depend  on  angular  rates  are  further  defined  in  a 
standard  form  of  dimensionless  angular  rate  in  the  form 


M. 
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The  rolling  moment  CL  is  defined  about  the  aerodynamic  roll  axis  X^,  which 
forms  the  intersection  between  the  disk  plane  and  the  aerodynamic  plane.  It 
is  considered  to  be  made  up  of  a  static  moment  due  to  spin  and  angle  of  attack, 
called  the  Magnus  moment,  and  the  roll  damping  due  to  roll  rate  about 
the  aerodynamic  roll  axis.  The  Magnus  moment  coefficient  is  normalized 
to  spin  rate  in  the  from  C.  ,  which  is  itself  a  function  of  angle  of  attack,  and 

the  roll  damping  is  of  the  standard  form  C.  .  These  coefficients,  as  is  the 

pitch  damping,  are  made.dimensionles  s  by  the  standard  aforementioned  approach 
of  normalizing  them  to  the  dimensionless  pitch,  roil,  or  spin  rates  (or  reduced 
frequency)  —  d/2V,  P^d/2V,  uid/2V.  In  the  case  of  roll,  the  parameter  is 

sometimes  called  the  equivalent  roll  helix  angle,  and  for  spin  it  is  usually 
called  the  peripheral  speed  ratio  and  sometimes  expressed  as  (Vr/V).  The 
use  of  the  d/2  in  these  reduced  frequencies  reflects  both  aircraft  practice  in 
the  case  of  the  roll  helix  angle  and  the  relationship  of  the  Magnus  effect  or 
airfoil  circulation  to  the  surface  velocity.  However,  many  ballisticians  do 
not  use  the  factor  of  2  in  their  definitions,  so  that  care  must  be  taken  in 
comparing  and  or  using  data  from  different  sources. 


The  yawing  moment  C/V  is  defined  about  the  positive  spin  axis  Z^,  and,  in 

this  case,  is  considered  to  consist  of  any  spin  torques.  For  a  symmetrical 

disk  with  no  driving  vanes,  the  aerodynamic  torque  coefficient  would  consist 

of  only  the  spin  damping  coefficient  (C  ).  No  data  are  available  on  the 

n,ii 

actual  magnitude  and  angle  of  attack  dependence  of  this  coefficient  since  its 
effect  was  undetected  in  tests  so  far.  If  there  are  aerodynamic  driving 
torques,  these  would  be  expressed  as  C 

o 
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Spin  torques  and/or  jet  damping  due  to  the  flare  exhaust  would  require  special 
definitions  tailored  to  the  assumed  characteristics  (usually  related  to  thrust 
and/or  mass  flow  rates). 

The  inertial  and  dynamic  axes  used  for  analysis  of  the  flight  dynamics  will  be 
presented  in  detail  in  Appendix  B. 

In  this  report,  the  word  camber  will  be  used  to  refer  to  the  curvature  of  the 
centerline  of  a  cross  section  of  the  disk.  In  airfoil  theory,  the  characteristics 
of  a  basic  finite  thickness  symmetrical  airfoil  may  be  considered  to  be  super¬ 
imposed  on  the  effects  of  the  curvature  and  orientation  of  an  infinitesimally 
thin  (or  flat  plate)  airfoil.  Most  low-speed  aircraft  have  thick  airfoils  whose 
camber  is  at  least  sufficient  to  provide  flat  lower  surfaces,  while  in  the 
Wright  brothers  aircraft  the  wings  were  essentially  cambered  flat  plate  air¬ 
foils. 

4.  Trajectory  Simulation  Development 

On  the  basis  of  the  separation  of  the  high  frequency  dynamic  motions  from 
the  slow  precession  and  linear  velocities  defining  the  trajectory,  a  quasi¬ 
steady  math  model  and  digital  computer  program  were  developed  to  simulate 
flare  trajectories.  This  program  allows  for  gyroscopic  precession  under 
the  aerodynamic  moments  and  translation  under  the  aerodynamic  forces 
experienced  as  a  function  of  angle  of  attack,  spin,  and  angular  rates. 

The  forces  are  considered  radial  (axial)  and  normal  to  the  disk,  and  the 
moments  are  considered  to  be  either  overturning  (pitching)  in  the  aerodynamic 
plane  or  rolling  about  a  radial  axis  in  the  plane.  The  program  allows  for 
spin  decay  under  air  friction.  It  allows  for  linear  decrease  of  mass  with 
time,  but  this  linearity  could  easily  be  modified  for  a  nonlinear  rate.  The 
mass  is  assumed  to  be  lost  from  the  inner  surface  of  an  annulus  of  pyro¬ 
technic.  The  aerodynamic  forces  and  moment  coefficients  defined  in 


-21- 


section  I.  A.  3.  may  be  tabulated  as  a  function  of  the  cavity  radius  of  this 
annulus  for  four  values  of  cavity  radius,  or  as  one  function  representing  a 
solid  body  or  a  burning  body  with  negligible  change  in  aerodynamics. 

The  aerodynamic  forces  and  moments  are  assumed  to  by  symmetrical  with 
respect  to  the  aerodynamic  plane  containing  the  velocity  vector  and  the  spin 
axis,  except  for  the  Magnus  moment  which  is  the  rolling  moment  induced  by 
the  interaction  of  spin  and  the  airflow.  At  the  present  time  this  Magnus 
rolling  moment  coefficient  is  assumed  to  be  a  parametric  fraction  of  the 
pitching  moment  coefficient.  The  fraction  can  be  varied,  or,  in  effect, 
the  aerodynamic  plane  could  be  considered  to  be  rotated  slightly  by  the  spin. 

The  program  can  be  commanded  to  simulate  a  nonprecessing  constant  attitude 
trajectory.  This  permits  definition  of  desired  reference  trajectories,  and  in 
this  mode  the  aerodynamic  moments  do  not  have  to  be  known  or  assumed. 

A  truncated  version  of  the  nonprecessing  trajectory  simulation  has  also  been 
used  for  special  in-house  preliminary  studies  on  a  teletype  console  time 
shared  computer.  The  unpublished  special  program  in  the  PACTOLUS 
language  was  expected  to  be  useful  for  remote  operation  from  Hurricane  Mesa, 
but  the  phone  lines /party  lines  would  not  have  permitted  it.  Its  use  was 
inspired  by  the  apparent  negligible  precession  in  the  first  series  of  Hurricane 
Mesa  tests. 

A  sample  of  some  early  trajectories  used  in  examining  the  initial  Hurricane 
Mesa  tests  and  planning  later  tests  is  shown  in  Figure  6.  The  shaded  area 
indicates  possible  variability  in  the  Hurricane  Mesa  face  at  impact.  Agree¬ 
ment  was  considered  reasonable  at  that  stage  of  development  for  the  50  foot 
per  second  launcher.  The  figure  indicated  that  it  would  require  a  150  foot  per 
second  launch  velocity  to  utilize  the  maximum  vertical  drop  available  at 
Hurricane  Mesa. 
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PITCH  ATTITUDE'  =  0  DEGREES 


5.  Preliminary  Flight  Dynamics  Analyses 


Trajectory  simulations  and  analytical  solutions  were  used  to  study  the  flight 
dynamics  of  the  flare.  These  investigations  involved  primarily  planar  tra¬ 
jectories  for  flares  which  were  assumed  not  to  p recess,  and  which,  therefore, 
required  only  lift  and  drag  force  predictions,  and  studies  of  precessing  bodies 
with  available  aerodynamic  moment  data.  The  studies  were  initiated  with 
estimated  data,  and  measured  data  were  used  as  they  became  available. 

Some  interesting  and  worthwhile  effects  were  observed  in  these  early  planar 
analyses  and  exploited  or  observed  physically  in  subsequent  investigations. 

It  was  found  that  the  orientation  of  the  disk  (held  constant)  was  the  dominant 
influence  in  determining  the  terminal  flight  behavior,  rather  than  the  initial 
flight  path  direction.  There  was  little  difference  for  a  level  disk  whether 
tossed  level,  upward,  or  downward,  but  if  tipped,  there  was  a  significant 
difference.  If  tipped  edgedown,  it  would  travel  farther  than  the  level,  and, 
if  tipped  edgeup,  it  would  go  less  far  or  even  ccme  back.  This  can  be 
observed  for  frisbees  and  clay  pigeons  also.  There  appeared  to  be  some 
optimum  orientation  for  gliding  rather  than  falling,  and,  in  one  case,  it 
appeared  capable  of  keeping  up  with  the  launch  aircraft.  Disks  rolled  to  the 
right  would  eventually  reach  a  straight  glide  path  to  the  right. 

The  effect  of  the  different  actual  angle  of  attack  histories  on  a  freely  precess¬ 
ing  disk  is  another  problem.  The  response  of  disks  permitted  to  precess 
was  disappointing  in  these  early  studies, especially  with  preliminary  aero¬ 
dynamic  estimates  which  turned  out  to  correspond  to  unrealistically  thin 
disks.  The  disks  would  precess  in  roll,  eventually  reaching  a  vertical  orien¬ 
tation  which  eliminated  the  angle  of  attack.  While  quite  nonlinear,  the  rapidity 
of  precession  increased  with  launch  speed  and  decreased  with  spin  or  inertia. 

The  most  baffling  results  were  those  preliminary  results  obtained  with  pre¬ 
liminary  clay  pigeon  data,  since,  although  no  one  was  sure  that  a  disk  would 
fly  aside  from  the  limited  Hurricam*  TMi-ya  tests,  a  clay  pigeon  does  success¬ 
fully  fly. 
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An  attempt  was  made  to  match  the  significant  parameters  of  a  typical  trap 
range  trajectory  using  photo  data  on  velocity  and  spin,  but  the  initial  pre¬ 
cession  under  the  unsymmetrical  moment  and  the  precession  and  travel  at 
impact  at  first  did  not  appear  to  be  consistent  with  the  clay  pigeon's  reputa¬ 
tion  of  successful  flight.  Of  course,  the  most  experienced  observers  of  clay 
pigeon  flight  were  the  least  likely  to  have  observed  a  large  number  of  com¬ 
plete  trajectories,  and  the  launch  dispersion,  wind  effects,  and  lack  of  quan¬ 
titative  historical  data  marie  it  impossible  to  define  departure  from  the  launch 
plane  at  impact.  As  a  matter  of  fact,  calculated  path  departure  was  within 
the  realm  of  reasonable  dispersion.  Later  results  (in  section  IV.  B.  1) 
showed  a  reasonable  precession  history. 

As  a  means  of  checking  the  simulation  results,  examining  dynamic  stability, 
and  possibly  obtaining  design  criteria,  the  equations  of  motion  were  linearized 
and  a  solution  attempt  was  made  for  stability  roots  and/or  dynamic  response. 
These  early  and  limited  analytical  studies  validated  the  assumption  of  motion 
separation  used  in  the  trajectory  analysis.  They  also  appeared  to  validate 
the  predicted  roll  precession  of  right  circular  cylinders  using  estimated  and 
wind  tunnel  data  by  identifying  the  primary  forcing  function  of  gravity  but  not 
any  compensating  stabilizing  terms.  They  did  not  correlate  with  the  well- 
known  behavior  of  frisbees  and  clay  pigeons  and  the  apparent  behavior  of  the 
dummy  flares  at  Hurricane  Mesa.  The  latter,  however,  changed  direction 
and  velocity  far  outside  the  bounds  of  the  linear  analysis  and,  of  course,  were 
considerably  thicker  than  the  configurations  for  which  data  were  available. 

The  preliminary  analytical  solutions  indicated  that  the  Magnus  (rolling) 
aerodynamic  moment  due  to  spin  was  a  stabilizing  influence  so  that  perhaps 
this  was  much  greater  than  anticipated.  They  also  indicated  that  the  smaller 
the  pitching  moment,  the  better,  so  that  perhaps  the  pitching  moment  was 
much  lower  than  measured  without  spin.  The  preliminary  analytical  solutions 
also  did  not  appear  to  properly  account  for  the  effects  of  the  camber  of  the 
clay  pigeon  and  frisbee  shapes. 

These  preliminary  analytical  solutions  were  later  extended  to  obtain  criteria 
which  did  account  for  and  explain  these  camber  effects,  which  correlate 
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with  observed  behavior,  and  which  serve  as  guidelines  for  future  efforts 
on  the  Self-Suspended  Flare.  These  criteria  are  presented  in  section  IV.  B.  2. 

6.  Fuzing  Investigations 

As  an  aid  to  planning  both  aerodynamic  and  flare  and  fuze  design  activities/ 
requirements,  a  cursory  examination  was  made  of  potential  fuzing  concepts. 
The  goal  was  to  provide  ignition  at  the  optimum  altitude  with  a  preset  sealed 
device  that  did  not  require  setting  at  launch.  The  preferred  approach  at 
NAD/ Crane  at  the  outset  was  to  calibrate  spin  decay  with  altitude  loss  from 
launch.  However,  because  of  indications  that  this  might  not  be  a  valid  fuzing 
signature,  other  concepts  were  examined.  Concepts  examined  and  discussed 
with  the  Project  Engineer  included  the  following: 

.  Spin  decay  (settable)  -  arming  by  a  combination  of  spin  velocity  and  spin 
acceleration,  and  firing  by  spin  decay. 

.  Mechanical  pyrotechnic  delay  (settable  in  six  increments)  -  arming  by 
a  borerider  pull  pin  and  spin  velocity,  and  firing  by  pyrotechnic  delay. 

.  Barometric  pressure  (settable  or  preset)  -  arming  by  spin  and  differ¬ 
ential  dynamic  pressures,  and  firing  by  barometric  pressure. 

These  concepts  were  not  carried  beyond  the  conceptual  stage,  but  helped  to 
focus  mutual  attention  on  requirements,  constraints,  possible  approaches, 
and  the  aerodynamic  contractor's  capability  for  providing  or  developing  what¬ 
ever  type  of  fuzing  device  might  be  required  at  any  stage  of  Self- Suspended 
Flare  development,  from  experimental  to  production. 

7.  Aerodynamic  Investigations 

Relatively  unsuccessful  literature  searches  were  made  to  obtain  aerodynamic 
data  on  flying  disks.  The  only  one  obtained  from  a  DDC  search,  an  Israeli 
report, ,  concentrated  on  the  purely  rotational  stability  theory,  with  results 
corresponding  to  those  of  the  writer,  but  it  neither  included  the  translational 
trajectory/ gravity  effects  nor  offered  aerodynamic  data.  A  paper  was 
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discovered  through  word-of-mouth  on  the  lift  and  drag  of  the  discus  which 
included  another  application  of  the  optimum  initial  orientation  phenomenon. 

As  the  Aerodynamic  Analysis  of  the  Self- Suspended  Flare  progressed,  some 
college  professors  indicated  informally  that  their  students  had  made  un¬ 
published  wind  tunnel  tests  on  frisbees  with  conventional  external  balances. 
The  available  data  had  varying  amounts  of  correlation  with  that  obtained  in 
this  study.  One  case  showed  the  only  spin  effect  to  occur  at  a  very  high  spin 
rate,  with  a  sudden  spike  in  the  data  which  does  not  correspond  to  any  of  the 
usual  spin  scaling  parameters  except  apparently  the  critical  Reynolds  number 
for  that  specific  condition. 

A  few  references  were  found  on  thin  disks  and/or  cylinders.  From  the 
available  data  base,  surprisingly  good  estimates  were  made  in  order  to 
initiate  the  dynamic  analyses  and  to  make  the  wind  tunnel  tests,  which  would 
be  the  prime  source  of  quantitative  data,  more  meaningful. 

When  completed,  these  wind  tunnel  data  were  scrutinized  and  estimates 
updated  to  cover  configurations  of  interest  for  which  tests  were  not  avail¬ 
able.  These  scrutiny/estimates  picked  up  the  trends  of  the  effects  of  thick¬ 
ness  and  corner  radius  leading  to  synthesis  of  new  candidate  test  configura¬ 
tions  expected  to  bracket  the  desired  effects.  While  the  effect  of  detailed 
corner  radius  on  basically  cylindrical  bodies  was  not  yet  isolated,  there 
seemed  to  be  a  definite  effect  of  thickness  and  of  gross  corner  radius. 

A  brief  study  was  made  of  the  effects  of  flare  plume  on  the  aerodynamics. 

With  the  magnitudes  of  flame  pressure  and  temperature  supplied  by  NAD/ 
Crane,  this  analysis  indicated  that  the  flare  plume  would  have  a  negligible 
effect  on  the  aerodynamics.  The  aerodynamic  characteristics  would  then 
correspond  to  the  equivalent  hollow  cool  body.  It  was  considered  desirable 
to  conduct  a  more  detailed  study  at  some  future  date,  but  such  a  study  was 
beyond  the  scope  of  this  project. 

An  example  of  the  aerodynamic  data  obtained  or  extrapolated  is  shown  in 
Figure  7.  More  detailed  results  are  presented  and  discussed  in  section 
IV.  B.  3. 
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♦  OVER 
TURNING 


Figure  7.  Pitching  Moment  Coefficient  for  Typical  Self- Suspended 
Flare  Configurations 
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8.  Preliminary  Hurricane  Mesa  Tests 


'tear  the  outset  of  the  aerodynamics  analysis,  a  preliminary  investigation 
cf  the  use  of  Hurricane  Mesa  for  free-flight  tests  was  conducted  by  NAD / 
Crane  and  the  Denver  Research  Institute.  Although  remote,  this  site  offered, 
in  addition  to  sheer  drops,  the  availability  of  electrical  power,  shops,  dark 
room,  launch  and  camera  stations,  access  to  both  top  and  impact  areas,  and 
a  sled  track. 

Dummy  flares  representing  right  circular  cylinders  12  and  8  inches  in 
diameter  and  2  inches  thick  and  expected  flare  density  were  launched  at 
maximum  velocity  and  spin  and  primarily  in  a  level  orientation.  One  was 
launched  with  the  launcher  rolled  about  15  degrees.  The  trajectories  were 
quite  repeatable  and,  in  general,  correlated  well  with  computer  calculations 
for  nonrolling  flares  within  the  limits  of  the  site  survey  data  available. 
One-half  pound, 4-inch  disks  were  not  correlated  because  the  launches  were 
not  repeatable.  They  were  more  subject  to  the  winds  and  were  not  used 
in  later  tests. 

The  dummy  flares  did  not  appear  to  roll  any  significant  amount.  Computer 
predictions  with  estimated  data  based  on  thin  disks  and  later  with  measured 
data  for  an  8;  1  diameter-to-thickness  disk  (section  IV.  B.  3.  )  indication  a 
relatively  severe  amount  of  rotation,  as  did  the  available  linear  theory. 

This  apparent  lack  of  correlation  provided  a  strong  incentive  for  the  later 
detailed  free-flight  and  wind  tunnel  testing.  With  hindsight,  the  differences 
could  be  largely  attributed  to  the  differences  in  thickness  assumed  in  the 
computer  studies  and  that  actually  flown.  The  wind  tunnel  model  was  for 
an  8:1  diameter-to-thickness  ratio,  the  8-inch  dummy  flares  were  4:1,  the 
12-inch  dummies  were  6: 1,  and  the  first  coefficient  estimate  was  for  essen¬ 
tially  infinite  ratio. 
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9.  Nonspinning  Wind  Tunnel  Tests 


During  the  early  planning  of  the  analysis,  it  was  anticipated  that  the  flare 
could  possibly  be  exhausting  from  both  the  top  and  bottom  of  the  grain  cavity, 
so  that  the  aerodynamic  configuration  would  be  an  annulus  with  increasing 
hole  size  (doughnut),  rather  than  an  inverted  cup  with  increasing  bottom 
cavity  diameter  as  was  actually  the  preferred  configuration. 

Support  of  either  a  spinning  annulus  or  a  hollow  cup  on  a  wind  tunnel  sting 
was  considered  beyond  the  scope  of  the  project.  Best  estimates  were  that 
the  spin-dependent  aerodynamic  terms  were  not  large  enough  to  justify  the 
expense  involved  in  devising  fixtures  to  measure  them.  Thus,  the  decision 
was  made  to  perform  initial  studies  with  nonspinning  models  in  order  to 
maximize  the  number  of  configuration  variations  examined. 

The  other  consideration  in  wind  tunnel  fixture  design  was  the  necessity  to 
place  the  "center"  of  the  strain  gage  balance  as  close  as  possible  to  the 
center  of  the  model  to  maximize  sensitivity  to  the  forces  and  moments.  This 
could  not  be  done  with  an  internal  balance  on  hollow  models.  Therefore,  a 
side  mount  was  used  to  place  the  pitch  measurement  gages  in  line  with  the 
pitch  axis.  This  would  also  make  it  easier  to  cover  the  full  range  of  0-90 
degrees  in  angle  of  attack. 

A  typical  mount  is  shown  in  Figure  8.  In  this  case  the  complex  clay  pigeon 
configuration  was  obtained  by  plastic  molding.  The  other  models,  the  right 
circular  cylinder  and  the  frisbee  shape,  were  machined  from  aluminum. 

The  results  of  these  tests  correlated  fairly  well  with  each  other  and  with 
preliminary  estimates  and  with  some  available  discus  data.  The  most 
severe  disagreement  with  preliminary  estimates  was  in  a  favorable  direction. 
The  center  of  pressure  of  the  thin  disk  was  assumed  relatively  far  forward 
so  that  the  estimated  pitching  moments  were  considerably  larger  than  those 
measured.  The  more  airfoil-like  bodies  had  a  more  pronounced  stall-peak 
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on  both  normal  force  and  pitching  moment  than  did  the  cylinder,  but  all  had 
similar  trends  before  and  after  stall.  As  expected,  from  conventional  airfoil 
practice  and  theory,  the  cambered  (domed)  clay  pigeon  and  frisbee  had  a 
definite  shift  of  the  symmetry  of  the  moment  curves  so  that  the  aerodynamic 
pitching  moment  was  zero  at  a  finite  angle  of  attack  rather  than  at  zero.  The 
hollow  counterparts  had  more  pronounced  reductions  of  moment  and  greater 
shift  of  the  moment  curve  intercepts. 

A  spring  restrained  oscillatory  side  mount  fixture  was  used  to  examine 
aerodynamic  pitch  damping.  The  damping  coefficients  were  quite  small  and 
of  questionable  sign.  Therefore,  the  results  may  have  been  masked  by  the 
nonlinear  hysterisis  type  instability  often  displayed  by  thin  disks  in  free 
oscillation  in  response  to  separation  from  the  body  itself  or  to  flow  turbulence. 

While  the  side  mount  provided  measurement  sensitivity  advantages,  it  also 
introduced  the  problem  of  flow  interference,  both  external  and  internal  (for 
hollow  models).  The  internal  interference  could  have  caused  the  measured 
moments  to  be  too  high  and  therefore  conservative.  The  interference  also 
could  have  caused  drag  results  to  be  too  high,  since  the  lift-to-drag  ratio 
appears  to  be  15  -  20  percent  too  low  to  correspond  with  informal  flight 
observations. 

10.  Interim  Investigations 

As  a  result  of  the  apparent  inconsistencies  between  preliminary  free-flight 
tests,  analytical  solutions,  and  simulation  results  with  available  nonspinning 
wind  tunnel  data,  it  was  concluded  that  either  there  was  a  greater  spin 
dependence  of  the  aerodynamics  than  anticipated,  causing  a  drastic  reduction 
in  overturning  moment,  or  there  was  a  stabilizing  aerodynamic  mode  ignored 
in  the  simulation. 

A  short  exploration  of  this  problem  area  and  of  approaches  to  the  solution 
was  undertaken.  This  included  observation  of  a  smoke  tunnel  flow,  observation 
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of  and  consideration  of  instrumentation  of  three -degree-of-freedom  wind 
tunnel  motion,  investigations  of  fixture  designs  for  wind  tunnel  measurements 
of  spinning  models,  and  consideration  of  free-flight  approaches.?.  These 
brief  explorations  were  undertaken  at  Honeywell  and  through  consultation  at 
Notre  Dame.  Consideration  was  also  given  to  magnetic  suspension  at  MIT. 

The  Notre  Dame  smoke  tunnel  flow  over  a  spinning  clay  pigeon  (Figure  9) 
showed  a  considerable  degree  of  three  dimensional  effect,  and  indicated  some 
slight  variation  with  spin,  but  not  enough  to  represent  a  drastic  reduction  in 
pitching  moment  from  that  of  a  nonspinning  body.  Three-degree-of-freedom 
observations  of  symmetrical  disks  on  Notre  Dame's  jewel  bearing  and 
Honeywell's  gas  bearing  displayed  unexpected  levels  of  attitude  stability  with 
slow  coning  which  do  not  agree  with  linear  theory  and  could  be  attributed  to 
mass  unbalance  or  the  sought  for  neutral  overturning  moment. 

This  three-degree-of-freedom  behavior  indicated  that  it  might  be  feasible  to 
support  a  model  for  3  degree-of-freedom  tests,  but  the  computer  program 
required  to  extract  the  aerodynamic  coefficients  was  not  available  and  would 
require  time  and  expense  to  develop.  The  programs  available  at  that  time 
assumed  solutions  for  quasi-symmetrical  bodies  such  as  projectiles  and 
missiles  with  spin  axes  near  the  velocity  vector  and  which  have  symmetrical 
differential  equations.  Because  these  symmetrical  equations  assume  com¬ 
parable  aerodynamic  stiffness  for  both  angular  rotations  of  interest  (pitch 
and  yaw),  they  are  not  applicable  to  gliding  disks  (pitch  and  roll)  or  gliding 
flettners  (yaw  and  roll)  in  which  there  is  no  aerodynamic  moment  due  to  roll. 
They  would  be  applicable  to  flight  with  the  flare  broadside  to  the  flow,  how¬ 
ever. 

Concurrent  discussions  with  Notre  Dame  regarding  free-flight  techniques 
brought  out  some  of  the  difficulties  and  possible  approaches  to  motion  simu¬ 
lation  and  data  extraction.  The  value  of  light  weight  models  for  amplification 
of  dynamic  motions  to  assist  in  instrumentation  was  emphasized. 
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The  effect  of  the  support  on  three -degree -of-freedom  data  is  another  area 
of  concern.  Discussion  of  magnetic  wind  tunnel  support  systems  with  Dr. 
Covert  of  MIT  indicated  that  present  systems  would  necessitate  small 
(quarter/half  dollar  sizes)  models,  and  additional  development  of  test  data 
extraction  techniques  would  be  required  to  obtain  six-component  data  from 
spinning  models. 

From  these  investigations,  it  appeared  that  the  least  additional  time  and  cost 
would  be  required  for  a  plan  involving  free-flight  tests  and  aerodynamic  data 
extraction  at  Hurricane  Mesa,  because  a  basic  project  and  equipment  were 
already  well  underway,  and  a  six-component  wind  tunnel  test  of  a  spinning 
solid  cylinder  model,  because,  if  successful,  quantitative  data  would  be 
directly  measured  to  define  the  effect  of  spin  on  the  overturning  moment. 

The  technique  proposed  would  increase  the  sensitivity  of  pitching  moment 
measurements. 

Subsequent  adoption  of  this  interim  plan  resulted  in  project  redirection  to 
participate  in  the  Hurricane  Mesa  tests,  reduce  the  test  data,  build  and  test 
a  spinning  wind  tunnel  feasibility  model,  and  incorporate  the  results  in  the 
dynamic  analysis  of  the  flare  behavior. 

11.  Hurricane  Mesa  Tests 

In  conjunction  with  the  NAD/Crane  Project  Engineer,  it  was  concluded  that 
the  most  expedient  approach  within  the  available  time  and  funds  would  be  to 
attempt  to  extract  dynamic  motion  and/or  aerodynamic  data  from  the  Hurri¬ 
cane  Mesa  free-flight  tests  which  were  already  scheduled  for  live  flare  tests 
and  formal  recording  of  nominal  trajectories.  This  approach  would  involve 
addition  of  special  models  and  test  conditions  and  additional  instrumentation. 
The  aerodynamics  contractor  contributed  planning  knowledge,  data  require¬ 
ments,  equipment  and  personnel  to  these  augmented  experiments  and  instru¬ 
mentation. 
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The  tracking  station  cameraman  employed  by  the  aerodynamics  contractor 
had  several  years  of  experience  as  a  cameraman  at  Hurricane  Mesa.  The 
aerodynamics  principal  investigator  made  on-site  observations  of  all  tests 
and  provided  selection  and  modifications  of  the  special  configurations. 

These  on-site  observations  and  experiment  selections  based  on  early  behavior 
observations  were  more  useful  than  the  preplanned  aerodynamic  experiments 
and  quantitative  data,  as  often  turns  out  to  be  the  case. 

The  tests  were  conducted  in  November  1969.  The  special  models  and  test 
conditions  devised  by  the  contractor  included  removal  of  steel  ballast  and/or 
wood  layers  from  some  of  the  standard  dummy  flares  prepared  by  the  test 
contractor,  removal  of  upper  corners  to  simulate  frisbee/clay  pigeon  con¬ 
tours,  edge  painting,  and  cross  hairs,  and  tilting  the  launcher.  Figure  10 
presents  typical  dummy  flare  configurations  (unpainted). 

In  addition  to  the  launcher  at  the  launch  site,  there  was  a  tracking  camera 
with  zoom  lens  above  the  launcher  on  a  forklift  (not  usable  for  live  firings), 
a  test  number  board  visible  to  tracking  cameras,  a  launch  signal  flash  bulb, 
and  camera  synchronizing  and  timing  electronics. 

At  the  control  site,  a  safe  distance  from  boresite  on  the  edge  of  the  cliff, 
was  located  the  launch  control  box  which  controlled  launch  and  ignition  event 
pulses  to  the  synchronizing  circuits  and  the  flashbulb.  At  the  tracking  sight, 
there  was  a  "tracking"  camera  with  telephoto  lens  mounted  on  a  panhead  with 
pickoffs  for  aximuth  and  elevation,  a  data  camera  recording  box,  panoramic 
motion  picture  camera,  and  associated  timing  electronics.  These  were 
supplemented  by  a  panoramic  still  camera  for  the  live  flare  tests. (See  Fig¬ 
ures  11,  12,  13,  14  and  15). 

Below  the  launch  site  on  the  sloping  shoulder  of  Hurricane  Mesa  was  placed 
a  grid  marking  at  1  degree  intervals  of  line-of-sight  depression  from  the 
boresight  station.  The  markers  were  5  foot  horizontal  poles  with  alternating 
1  foot  red  and  white  stripes. 
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Figure  14.  Top  View  of  Flare  Launcher 
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For  the  particular  system  evolved,  the  maximum  possible  data  was  obtained. 
One  feature  contributing  to  maximizing  data  was  the  use  of  color  film  on  all 
the  tracking  cameras  to  provide  contrast  with  the  rugged  terrain. 

The  results  of  the  Hurricane  Mesa  tests  were  somewhat  spectacular  and 
unexpected.  Some  of  the  modified  dummy  flares  greatly  amplified  the  pre- 
cessional  motions,  which  were  not  perceptible  for  the  standard  dummy  flares 
in  these  or  earlier  preliminary  tests, instead  of  the  hoped  for  stable  coning 
modes.  The  live  flares  developed  sufficient  thrust  to  climb  above  the  mesa 
instead  of  falling  at  the  predicted  rate  along  the  basic  parabolic  path  of  the 
dummy  flares.  Lighter,  thinner  dummy  flares  rolled  to  the  left  on  the  order 
of  10-30  degrees,  thus  curving  the  path  significantly  from  the  nominal  plane. 
When  launched  at  a  downward  angle,  they  rolled  to  the  left  until  edge-on 
after  even  greater  trajectory  curvature.  Even  at  full  or  greater  weight, 
dummies  with  rounded  corners;  that  is  solid  frisbee/clay  pigeon  shapes, 
peeled  off  quite  significantly.  In  some  cases,  the  initial  rapid  roll  preces¬ 
sion  appeared  to  stop  or  slow  down  and  the  disks  attempted  to  glide  in  a 
manner  corresponding  to  this  new  orientation,  as  predicted  by  nonprecessing 
nonrolling  theory  in  computer  runs.  As  a  matter  of  fact,  one  case  was 
"gliding"  so  well  that  it  was  airborne  at  least  10  seconds  even  though  it 
impacted  near  the  base  of  the  cliff  instead  of  the  usual  6-7  second  drop  to  the 
nominal  impact  area  down  the  sloping  shoulder.  These  manifestations  of 
steady  precession  without  visible  wobble  tended  to  validate  the  precessional 
model  used  in  the  trajectory  computer  program,  especially  since  closer 
scrutiny  and/or  parametric  analysis  indicated  that  the  scaled  precession  of 
standard  dummy  flares  might  not  be  readily  perceptible  under  the  Hurricane 
Mesa  test  conditions. 

Continual  strong  updrafts  made  it  impossible  to  conduct  instrumented 
experiments  with  standard  toy  or  oversized  toy  frisbees.  The  updrafts  also 
tend  to  cloud  the  quantitative  data  reduction  for  the  dummy  flares. 
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Meaningful  observations  of  clay  pigeon  flight  which  tend  to  validate  the 
criteria  developed  were  obtained  in  spite  of  the  high  winds,  however.  It  was 
relatively  easy  to  find  a  stable  nose -down  orientation  which  would  glide 
"wings  level"  at  a  stable  equilibrium  velocity.  The  straight  line  flight  path 
would  undergo  changes  in  direction  as  wind  direction  shifted  with  altitude, 
but  few  cases  of  precession  or  curving  back  into  the  mesa  were  observed. 

It  was  guessed  that  some  had  traveled  1/4  mile  before  impact  far  down  on 
the  base  of  the  mesa.  Unfortunately,  no  quantitative  data  were  obtainable 
on  these  flights. 

In  all  analyses,  it  has  been  assumed  that  the  flare  would  generate  negligible 
thrust  on  the  basis  of  illumination  requirements  and  preliminary  pressure 
estimates  at  NAD/ Crane.  Due  to  a  combination  of  factors,  most  of  the  live 
flares  generated  sufficient  thrust  to  halt  the  fall  of  the  flare  and  lift  it  above 
the  top  of  the  mesa.  The  one  flare  which  did  not  rise,  of  course,  provided 
considerable  illumination  of  the  face  of  the  mesa  (Figure  16).  The  rising 
flares  are  typified  by  Figure  17.  The  flares  were  apparently  quite  stable 
during  both  descending  and  ascending  flights  although  there  was  an  impression 
of  tilting  toward  the  tracking  site  that  caused  a  lapse  of  tracking  on  one  flight. 

The  launcher  was  developed  by  the  flight  test  contractor  as  originally  speci¬ 
fied  by  NAD/Crane  for  carriage  in  a  light  aircraft  so  that  its  primary  pur¬ 
pose  was  to  impart  a  controllable  pre-spin  and  sufficient  velocity  for  safe 
separation.  It  uses  a  pneumatic  cylinder  to  drive  a  carriage  which  mounts 
a  turntable  on  which  the  flare  is  placed.  The  launcher,  consisting  of  a  rubber 
belt  wrapped  around  the  disk,  a  turntable,  and  an  idler,  imparts  1000  rpm 
to  a  12  inch,  12  pound  disk  while  the  carriage  imparts  a  50  foot  per  second 
velocity. 

This  velocity  would  seem  to  provide  a  safety  margin  for  early  air  launches 
and  did  provide  a  free-flight  spinning  environment  for  the  live  flares,  but 
it  was  not  sufficient  to  provide  a  ground  based  simulation  of  air  launch.  It 
did  not  use  the  range  capability  of  the  Hurricane  Mesa,  although  greater 
range  would  have  strained  the  instrumentation  capability. 
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Almost  a  mile  of  movie  film  plus  numerous  still  negatives  were  obtained 
from  the  test  series  (48  official  launches). 

The  trajectories  obtained  for  standard  dummy  flares  agreed  fairly  well  with 
predictions  for  the  launch  conditions,  and  the  observed  precessions  tended 
to  validate  the  precessional  computer  program.  There  was  reasonable 
correlation  with  the  analytical  criteria  solution,  but  the  rapid  change  in 
angle  of  attack  for  the  low  speed  launch  and  velocity  variations  makes  it 
difficult  to  correlate  such  linearized  analyses  closely  and  makes  it  impossi¬ 
ble  to  use  these  solutions  to  extract  the  desired  aerodynamic  coefficients. 

A  better  simulation  would  have  required  much  higher  launch  velocities  (per¬ 
haps  from  a  sled  or  aircraft). 

While  not  successful  in  producing  quantitative  data,  not  too  much  had  been 
expected,  and  these  tests  were  valuable  because  they  did  provide  qualitative 
validation  of  the  math  models  and  insight  into  the  dynamic  criteria,  revealed 
problems  in  flare  design,  gave  experience  in  field  testing  and  data  collection, 
and  led  to  experience  in  data  reduction. 

12.  Hurricane  Mesa  Data  Reduction 

Plans  for  reducing  the  Hurricane  Mesa  test  films  by  the  aerodynamics  con¬ 
tractor  were  predicted  on  the  availability  of  film  reading  equipment  with  an 
automatic  computer  input  interface  and  personnel  familiar  with  film  data 
extraction,  with  Hurricane  Mesa  itself,  and  with  the  math  models  for  which 
the  data  were  intended  and  'i.e  availability  of  large  and  small  computers  and 
their  associated  programmers.  The  film-reader  and  interface  could  record 
crosshair  position  and  angle  or  change  therein  and  frame  number  or  time. 

It  soon  became  apparent  that  it  was  not  feasible  to  write  a  master  program 
which  would  accept  inputs  from  all  cameras  and  print  out  position,  attitude, 
and  aerodynamic  data,  even  though  the  maximum  obtainable  data  and  generally 
high  quality  photography  were  generated.  Each  test  had  to  be  examined 
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separately,  and  the  most  appropriate  combination  of  film  data  correlation 
assumptions  and  computation  applied.  This  often  involved  multiple  readings 
of  the  same  camera  and/or  use  of  data  from  another  camera  to  interpret  or 
compute  a  set  of  parameters. 

It  was  recognized  fairly  early  that  no  quantitative  attitude  data  would  be 
obtained  because  no  operator(s)  could  consistently  locate  and  measure  the 
axes  of  the  elliptical  image  of  the  flare,  which  varied  from  a  line  to  a  near 
circle. 

As  could  be  expected  in  such  a  system,  there  were  lapses  in  synchronization 
or  timing  signals,  but  conditions  were  repeatable  enough  to  justify  attempts 
to  reconstruct  or  ignore  these  data  in  many  cases. 

It  was  found  that  the  most  convenient  approach  to  trajectory  position  data 
extraction  was  terrain  matching.  This  was  especially  true  for  the  bo  resight 
camera  because  the  curving  flight  path  and  narrow  fields  of  view  caused 
drastic  departure  of  the  flare  image  from  the  test  range  grid,  which  itself 
was  limited  by  terrain.  It  was  relatively  easy,  however,  to  recognize  the 
prominent  terrain  features  in  the  motion  picture  films  and  to  locate  them 
on  still  pictures  made  of  the  terrain  from  both  tracking  and  boresight  stations. 
Coordinates  of  the  images  on  these  stills  were  corrected  for  camera  orien¬ 
tation,  itself  extracted  from  image  position  of  reference  points.to  obtain 
measures  of  line  of  sight  "angles.  "  These  angles  for  boresight  and  tracking 
were  combined  to  determine  three-dimensional  flare  position  coordinates. 

The  boresight  calculations  required  zoom  lens  calibration  of  focal  length 
variation  with  time  from  event  and  image  size  measurement. 

The  planar  trajectories  obtained  for  standard  inert  flares  agreed  reasonably 
well  with  previous  computer  predictions,  but  attempts  to  extract  detailed 
motion  parameters  from  these,  and  especially  from  precessing  special 
flares,  were  not  particularly  successful.  This  disappointing  return  on  effort 
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expended  reflected  a  combination  of  factors  which  included  data  scatter, 
intermittent  gaps  in  data,  poor  resolution  near  the  launch,  strong  updrafts, 
unknown  atmospheric  density,  the  difficulties  in  instrumentation  and  calibra¬ 
tion  due  to  the  terrain,  impracticality  of  reading  an  adequate  number  of 
points  for  accurate  data  fitting,  and  the  multiple  differentiation  required 
and  the  unavailability  of  standard  software  for  these  special  cases.  All 
parties  were  aware  of  the  gamble  involved. 

13.  Spinning  Model  Wind  Tunnel  Tests 

A  spinning  model  was  made  and  tested  to  obtain  data  and  check  the  feasibility 
of  the  proposed  technique.  The  model  consisted  of  a  nonspinning  center  ring 
to  which  the  rear  sting  was  attached  and  upper  and  lower  spinning  disks 
connected  by  a  through  shaft  (see  Figure  18).  These  disks  were  spun  by  a 
cold -gas  turbine  mounted  on  the  shaft.  The  8  inch  diameter  model  tested 
had  a  4: 1  thickness  to  diameter  ratio  so  that  it  was  essentially  a  full  sized 
mockup  of  the  8  inch  diameter  live  flares  launched  at  Hurricane  Mesa.  The 
Hurricane  Mesa  tests  indicated  that  the  thicker  shape  would  be  more  stable, 
and,  as  expected,  the  pitching  moment  was  smaller  and  virtually  neutral  at 
lower  angles  of  attack.  The  spin  effect  itself  was  small.  The  internal 
mechanism  precluded  use  of  much  thinner  shapes  or  large  cavities,  but  a 
variety  of  corner  radii  and  bottom  cavity  contours  can  be  accommodated 
with  modifications  of  the  upper  and  lower  disks.  The  turbine  will  require 
modification  if  higher  spin  rates  are  desired. 

This  model  is  suitable  for  both  spinning  and  nonspinning  tests. 

This  model  reflected  our  experience  that  the  level  of  moments  being  con¬ 
sidered  cannot  be  accurately  measured  with  balances  that  are  displaced  from 
the  model  center  let  alone  exterior  to  the  model.  Thus,  the  accurate  meas¬ 
urement  of  moments  on  a  thin  hollow  body  is  still  a  problem  area. 
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Figure  18.  Spinning  Wind  Tunnel  Model  Installation 


To  maximize  sensitivity  in  pitch  in  order  to  detect  the  anticipated  spin  effect 
on  pitching  moment,  the  strain  gage  balance  was  placed  to  the  side  of  the  bearing 
assembly  so  that  its  pitch  center  was  aligned  with  the  spin  axis  (see  Figure  19). 
This  meant  that  the  roll  center  was  displaced  approximately  1  inch  and  pro¬ 
duced  strong  coupling  of  the  normal  force  into  the  roll  gage,  making  it 
extremely  difficult  to  detect  small  Magnus  moments.  Since  the  balance  could 
not  be  placed  at  the  center  because  of  the  bearing  i.nd  turbine,  solution  of 
this  sensitivity  conflict  with  standard  six-component  balances  would  require 
either  two  separate,  partially  redundant,  balances  or  two  passes  with  the 
same  balance  relocated  between  passes. 

Existing  so  called  Magnus  balances  are  designed  for  the  stiffness /sensitivity/ 
geometry  requirements  of  symmetrical  projectiles.  Consequently,  it  would 
require  a  specially  designed  and  built  yoke  and  gage  assembly  with  the  gages 
separated  to  provide  the  required  sensitivity  for  this  special  configuration. 

This  delicate  device  could  be  relatively  expensive. 

Indications  are  that  the  Magnus  moment  must  have  been  insignificant  on  the 
shapes  tested  so  far  but  that  it  could  become  more  important  for  the  thicker 
shapes  envisioned  as  a  desirable  direction  for  study.  Consequently,  con¬ 
sideration  should  be  given  to  construction  of  such  a  special  device  for  future 
investigations. 

14.  Stability/ Precession  Criteria 

On  the  basis  of  the  observations  of  Hurricane  Mesa  flight  of  the  standard 
and  special  dummy  flares  and  clay  pigeons,  the  preliminary  analytical 
solutions  were  updated  and  re-examined.  The  test  observations  revealed  no 
evidence  of  a  dynamic  mode  stable  about  the  vertical,  but  rather  steady 
precession  "proportional"  to  aerodynamics  and  gyroscopic  stiffness.  The 
re-examination  revealed  the  sought  for  effect  of  camber  which  both  provided 
nonprecessing  flight  for  the  actual  clay  pigeons  and  frisbees,  but  caused  the 
precession  of  the  solid  frisbee/clay  pigeon  shaped  dummy  flares,  recon¬ 
firmed  the  forcing  effect  of  gravity  drop,  and  reinforced  the  validity  of  the 
trajectory  math  model. 
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These  analysis  provided  an  explanation  of  the  flight  of  the  toy  frisbee  and 
tentative  criteria  for  application  of  the  principle  to  ordnance  delivery. 

These  preliminary  criteria  presented  in  section  IV.  B.  2.  represent  an  inter¬ 
action  of  configuration  spin  and  delivery  conditions  which  are  easily  satisfied 
for  the  frisbee  and  the  clay  pigeons  tossed  from  Hurricane  Mesa.  Satisfaction 
for  dense  flares  or  munitions  may  require  a  brute  force  approach  if  the 
frisbee  analogy  is  not  practical. 

The  criteria  may  be  summarized  as  follows. 

A  spinning  thin  disk  will  precess  (roll)  about  its  velocity  vector  if  there  is  an 
aerodynamic  moment  resulting  from  its  inclination  to  the  flow.  A  cupped  or 
domed  disk  will  have  a  neutral  aerodynamic  moment  at  a  particular  inclina¬ 
tion  (angle  of  attack)  at  which  angle  it  will  also  generate  finite  lift. 

If  the  lift  equals  the  weight  at  this  angle  of  attack,  the  disk  will  stay  at  that 
angle  of  attack  without  precessing. 

If  the  orientation  and  glide  angle  are  such  that  the  velocity  stays  constant, 
this  condition  could  continue  indefinitely. 

15.  Trajectory  Computations 

Subsequent  to  the  wind  tunnel  tests,  the  trajectory  program  was  exercised 
with  data  and  launch  conditions  derived  from  the  4: 1  cylinder  wind  tunnel 
tests,  the  Hurricane  Mesa  tests,  the  frisbee /pigeon  tests,  and  the  insights 
gained  from  scrutiny  of  the  linearized  precessional  criteria.  There  were 
varying  degrees  of  correlation  with  flight  observations  and  criteria  inter¬ 
pretations.  See  Section  IV.  B.  1.  for  typical  results. 

For  the  clay  pigeon/ frisbee,  results  were  obtained  for  simulated  trap  range 
or  downhill  tosses  that  appear  to  be  within  the  limits  of  visual  observation 
under  the  normal  terrain  and  wind  conditions  and  the  variations  in  measured 
lift  and  drag. 
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Clo3er  examination  of  the  roll  precession  for  simulated  standard  Hurricane 
Mesa  launches  of  standard  flares  led  to  the  observation  that  the  computed 
roll  angles  near  impact  could  very  well  not  be  detectable  at  the  distances 
involved  even  with  a  less  angular  background. 

Attempts  to  find  a  launch  angle /velocity/ spin  combination  that  would  give 
negligible  precession  by  gliding  at  a  neutral  point  for  the  4: 1  cylinder  were 
not  completely  successful  because  the  aerodynamic  moment  was  apparently 
still  too  large  during  the  transient  to  the  neutral  point.  It  is  hoped  a  thicker 
disk  will  give  the  desired  level. 

The  effect  of  launch  angle  was  not  correlated  with  the  Hurricane  Mesa  re¬ 
sults.  However,  the  roll  increase  could  be  attributed  to  the  fact  that  the 
level  launch  would  cause  the  disk  to  pass  through  a  stall  point  and  reach 
the  lower  moments  attendant  to  high  angles  of  attack,  while  the  inclined 
launches  would  stay  at  lower  angles  of  attack  and  high  moments  longer. 

16.  Recent  Efforts 

Results  were  not  yet  conclusive  concerning  feasibility  or  required  configura¬ 
tions.  On  the  basis  of  the  trends  established  by  the  wind  tunnel  data  and  the 
alternatives  offered  by  the  criteria,  an  outline  of  the  next  iterations  in  con- 

V 

figuration  was  drawn  up.  These  configurations  included  increases  in  thick¬ 
ness,  variations  in  corner  radii,  variations  in  depth  of  a  partial  cavity  which 
would  still  leave  pyrotechnic  volume,  and  combinations  thereof.  These 
studies  would  hopefully  isolate  the  separate  effects  which  were  combined  in 
the  typical  candidate  configurations  tested  so  far.  Funds  did  not  permit 
pursuit  of  this  or  of  any  other  attacks  on  solution  so  that  they  must  be 
reserved  for  future  efforts. 

Discussions  were  held  with  the  NAD/Crane  Project  Engineer  regarding  air 
drops.  It  was  concluded  that  the  most  expedient  approach  to  simulation, 
feasibilit>  evaluation,  and  test  instrumentation  development  would  be  preliminary 
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drops  of  inert  flares  at  NAD/Crane  with  available  aircraft  and  manual 
tracking  photography.  This  would  reveal  the  problems  to  be  anticipated  in 
more  formal  instrumented  drops  of  inert  and  live  flares. 

These  technical  efforts  have  been  suspended  pending  further  authorization 
and  flight  test  activity. 

It  should  be  noted  that  tracking  photographs  of  recent  gas  gun  tests  using 
spinning  circular  sabots  launched  edge -on  with  the  disk  vertical  showed  a 
tendency  to  roll-out  towards  "wings  level,  "  nose  downhill  flight  in  a  manner 
analogous  to  that  of  a  frisbee  tossed  on-edge.  Exploitation  of  this  prelimi¬ 
nary  impression  obtained  from  light-weight  sabots  to  dense  flares  was  not 
possible  since  it  was  obtained  after  suspension  of  technical  activity.  How¬ 
ever,  it  bears  consideration  in  future  studies  as  a  possible  means  of  auto¬ 
matically  compensating  for  the  launch  transient  that  would  otherwise  have 
to  be  calibrated. 
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B.  RESULTS 


Typical  results  of  the  efforts  described  in  section  IV.  A  are  presented  in  the 
following  subsections  collected  under  the  following  categories: 

.  Theoretical  trajectories. 

.  Stability/precession  criteria. 

.  Configuration  aerodynamics  characteristics. 

.  Free-flight  testing. 

.  Math  models. 

These  results  are  supplemented  by  the  appendices  to  this  final  report. 

1.  Theoretical  Trajectories 

During  the  aerodynamic  analysis,  computer  simulations  were  performed  on 
the  flight  trajectories  of  prototype  flare  concepts  with  a  range  of  constraints 
on  orientation  and  mass.  Most  of  the  prototype  flare  studies  were  made  with 
nonprecessing  modes  of  the  simulations  to  determine  idealized  reference 
planar  trajectories  while  searching  for  a  configuration  which  would  fly 
stably  and/or  for  the  criteria  it  must  satisfy.  Since  there  could  be  an  in¬ 
definite  altitude  loss  between  launch  and  desired  ignition  altitudes,  the 
trajectory  behavior  of  inert  flares  was  studied.  Since  it  was  initially 
assumed  that  burning  pressures  would  not  produce  thrust,  simulations  were 
run  with  mass  loss  corresponding  to  burning  with  and  without  thrust  and 
with  and  without  effects  of  this  volume  loss  on  the  aerodynamics.  Because 
their  supposedly  well  known  behavior  helped  inspire  the  concept  and  their 
contours  represented  easily  identifiable  candidate  flare  configuration 
variations,  simulations  of  the  flight  of  clay  pigeons  and  toy  frisbees  were 
utilized  to  check  the  computer  math  models  and  analytical  criteria  developed. 
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a.  Typical  Nonprecessinn  Inert  Flare  Trajectories  -  A  typical  paramet¬ 
ric  analysis  of  the  idealized  flight  of  the  flare  prior  to  ignition  is  shown  in 
Figure  20  as  a  function  of  launch  velocity  for  falls  from  launches  at  10,  000 
feet  to  ignition  at  3000  feet  or  above.  In  this  example,  a  12  pound,  12-inch- 
diameter,  2-inch-thick  flare  was  assumed  to  be  launched  in  level  flight  in  a 
level  orientation  and  held  gyroscopically  in  that  orientation  throughout  the 
flight.  Similar  curves  were  obtained  for  other  sizes. 

b.  Typical  Trajectories  for  Clay  Pigeons  and  Frisbees  -  In  order  to 
validate  the  precessional  computer  program  and  the  precession  criteria 
of  Section  IV. B.  2  relating  to  both  aerodynamic  characteristics  and  delivery 
conditions,  the  following  runs  were  made  using  the  measured  aerodynamic 
data. 

Figure  21  shows  the  results  of  the  simulation  of  a  toy  frisbee  toss  from 
a  50  foot  elevation,  with  initial  orientation  as  prescribed  by  the  theory, 
including  a  slight  droop  to  the  left  as  one  naturally  tosses  a  frisbee.  Its 
path  is  relatively  straight  and  level  (Figure  21).  The  peeloff  (Figure  22) 
which  does  occur  is  well  beyond  the  point  at  which  a  normal  flight  would 
strike  the  ground. 

It  is  felt  that  the  computed  glide  path  in  Figure  21  was  steeper  than  normally 
expected  for  a  frisbee  for  two  reasons.  The  wind  tunnel  model  mounting 
may  have  given  a  lower  than  actual  lift-to-drag  ratio.  Those  shallow 
recreational  tosses  which  do  not  appear  to  peeloff  probably  strike  the 
ground  first  anyhow,  and  many  do  peeloff.  This  peeloff  is  explained  by 
the  analytical  model  as  due  to  slowdown  because  the  flight  is  at  a  shallower 
angle  than  that  required  for  equilibrium  glide. 

Figures  23  and  24  illustrate  the  simulated  behavior  of  the  clay  pigeon  on  a 
trap  range.  Even  through  it  never  reaches  a  steady- state  glide  orientation 
nor  satisfies  the  available  steady- state  precession  criteria,  it  did  maintain 
a  relatively  straight  path  and  relatively  constant  orientation  in  the  time 
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Figure  20.  Typical  Idealized  Self- Suspended  Flare  Trajectories  —  Nonburning 
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Figure  21.  Typical  Frisbee  Trajectory  Simulation  Results  -  Position 
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Figure  22.  Typical  Frisbee  Trajectory  Simulation 
Results  -  Orientation  Angles 
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Figure  23.  Simulated  Trap  Range  Flight  of  Clay  Pigeon  -  Position 
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Figure  24.  Simulated  Trap  Range  Flight  of  Clay  Pigeon  - 
Orientation  Angles 
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before  impact.  The  launch  velocity  and  angle,  peak  altitude,  and  impact 
point  agree  fairly  well  with  preliminary  motion  picture  data  obtained  at  an 
actual  range,  and  the  roll  angle  and  cross  range  departure  are  considered 
below  the  level  of  observation  on  a  range  subject  to  winds,  gusts,  erratic 
launchers,  and  informally  instrumented  tests.  Again,  the  behavior  may 
have  been  exaggerated  slightly  by  interference  effects  on  the  measured 
aerodynamics,  but  the  computed  behavior  is  considered  to  validate  the  math 
model  itself. 

An  even  more  dramatic  validation  of  the  math  model  and  the  precession 
criteria  is  shown  in  Figures  25  and  26  for  the  simulated  toss  from  a  cliff 
of  a  clay  pigeon  at  the  prescribed  orientation  for  stable  glide  and  near  the 
prescribed  glide  velocity.  In  this  simulation,  quite  straight  and  "wings 
level"  flight  was  achieved  for  an  indefinite  period.  This  result  correlates 
with  informal  trials  during  the  Hurricane  Mesa  tests.  Clay  pigeons  tossed 
downhill  over  the  cliff  flew  stably  in  straight  line  segments  which  appeared 
to  correspond  to  wind  direction  changes  and  reached  the  base  of  the  mesa 
without  turning  into  the  cliff.  Unfortunately,  no  photographic  record 
of  this  is  available,  but  several  of  the  observing  personnel  watched  these 
experiments  while  awaiting  better  weather  for  testing  and/or  photography. 

c.  Simulated  Mesa  Toss  of  Inert  Flares  -  Figures  27  through  30  present 
the  results  of  the  application  of  the  trajectory  simulation  and  the  estimated 
aerodynamics  based  on  the  available  data  to  the  Hurricane  Mesa  launches 
of  the  12-inch,  12-pound  and  8-inch,  5-pound  dummy  flares,  respectively. 
As  indicated  by  the  impact  point  from  flight  test  35,  there  was  reasonable 
validation  of  the  model  and  the  data  for  the  12  inch,  12  pound  unit.  The 
behavior  of  the  8 -inch,  5-pound  unit,  similar  to  that  of  the  live  flares 
launched  at  Hurricane  Mesa,  corresponded  to  the  simulation,  and  if  any¬ 
thing,  appeared  slightly  more  stable  at  Hurricane  Mesa. 

The  abrupt  peeloff  shown  in  the  computed  results  occurred  at  times  and 
falls  beyond  that  required  for  impact  at  Hurricane  Mesa.  The  computed 
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Figure  25.  Typical  Simulation  Results  for  Cliff  Toss 
of  Clay  Pigeon  -  Position 


Figure  26.  Typical  Simulation  Results  for  Cliff  Toss  of  Clay  Pigeon 
Orientation  Angles 
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Figure  27.  Simulated  Hurricane  Mesa  Trajectory  for  12-inch,  12-Pound 
Dummy  Flare  -  Position 
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Figure  28.  Simulated  Hurricane  Mesa  Trajectory  for  12-Inch,  12-Pound 
Dummy  Flare  -  Orientation  Angles 
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Figure  29.  Simulated  Hurricane  Mesa  Trajectory  for  8-Inch, 
5- Pound  Dummy  Flare  -  Position 


ORIENTATION  ANGLES  -  DEGREES 


Figure  30.  Simulated  Hurricane  Mesa  Trajectory  for  8-Inch 
5- Pound  Dummy  Flare  -  Orientation  Angles 
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departure  from  a  planar  trajectory  prior  to  this  impact  time  was  considered 
below  the  threshold  of  visual  observation  against  the  rugged  Hurricane  Mesa 
terrain.  However,  this  prediction  of  extreme  trajectory  curvature  for 
greater  falls,  even  at  this  low  speed,  raises  questions  concerning  per¬ 
formance  at  actual  launch  speeds  and  altitudes  and/or  the  validity  of  the 
models  and  data  which  can  apparently  only  be  answered  by  more  realistic 
flight  simulation. 

d.  Simulated  Air  Drop  of  Inert  Flare  -  Figures  31  and  32  illustrate  the 
application  of  the  trajectory  program  and  the  associated  aerodynamic  data 
for  the  8-inch,  5-pound  dummy  flares  to  a  higher  launch  speed  (175  feet  per 
second),  corresponding  to  the  very  low  end  of  the  air  launch  spectrum. 

This  shows  a  200  foot  cross  range  excursion  of  the  flare  in  the  first  10 
seconds  and  a  30  degree  excursion  in  roll,  both  of  which  would  be  percep¬ 
tible  to  an  observer.  It  is  interesting  to  note  that  in  Figures  27  through  32, 
the  roll  angle  has  seemed  to  approach  a  steady- state  value  which  would  tend 
to  result  in  a  steadily  turning  flight  until  the  angle  of  attack  approaches 
equilibrium. 

If  this  simulation  result  is  valid,  the  flight  performance  would  probably  be 
unacceptable,  but  there  are  enough  inconsistencies  and  missing  aerodynamic 
data  to  indicate  the  need  for  more  free-flight  simulation. 

e.  Simulated  Live  Flare  Trajectories  -  As  displayed  vividly  in  the  live 
tests  at  Hurricane  Mesa,  generation  of  thrust  by  the  flare  could  provide  a 
truly  self- suspended  or  over- suspended  flare.  Most  of  the  trajectory  analy¬ 
ses  performed  were  based  on  the  pyrotechnic  design  intent  of  negligible 
thrust,  however. 

(1)  Effect  of  Mass  Loss  on  Trajectory  -  Figures  33  and  34  present 
some  typical  trajectories  which  were  computed  to  show  the  effect  of  mass 
loss  due  to  burning  on  the  trajectory  of  a  stable  flare. 


*. 
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Figure  31.  Simulated  Air  Launch  of  8-Inch,  5-Pound  Dummy  Flare 
Position 
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Figure  32.  Simulated  Air  Launch  of  8- Inch,  5- Pound  Dummy  Flare 
Orientation  Angles 
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Figure  33.  Typical  Effect  of  Flare  Mass  Loss  on  Trajectory 
of  Nonprecessing  12-Inch,  12-Pound  Flare 
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Figure  34.  Typical  Effect  of  Flare  Mass  Loss  on  Trajectory 
of  Nonprecessing  8-Inch,  5-Pound  Flare 
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These  figures  are  for  an  arbitrary  launch  velocity  of  200  feet  per  second, 
which  is  a  likely  speed  for  preliminary  live  air  drops  and  for  light  recon- 
naisance  aircraft.  They  cover  typical  diameter,  weight,  and  burning  time 
combinations.  On  the  figures,  overlapping  curves  are  shown  for  ignition 
delays  of  0,  2,  4,  6  and  10  seconds,  while  infinite  ignition  time  refers  to  the 
nonburning  case  for  comparison.  The  nonburning  cases  have  end  points 
corresponding  to  burning  time  for  the  burning  items. 

In  these  idealized  planar  trajectory  simulations,  the  thrust  was  assumed  to 
be  negligible,  precession  was  assumed  to  be  negligible  with  respect  to  a 
level  orientation  and  launch,  and  the  effect  of  burning  and  mass  loss  and 
cavity  on  the  aerodynamics  was  assumed  to  be  negligible  on  the  same  basis 
as  the  assumption  of  negligible  thrust.  Because  the  aerodynamics  did  not 
change,  the  net  effect  displayed  was  a  balanced  effect  on  accelerations  in  both 
lift  and  drag  directions.  Burning  just  reduced  the  arc  length  distance 
traveled  compared  to  the  distance  along  the  paih  of  an  equivalent  nonburn- 
ing  flare  for  a  given  amount  of  time. 

This  reduction  was  on  the  order  of  1500  feet  and  1000  feel  of  :.r< 
the  12-inch,  12-pound  and  8-inch,  5-pound  flares,  respect  ire] 

In  other  calculations,  the  inclusion  of  cavity  effect  on  lift  resulted  in 
some  deviation  above  the  nonburning  path.  This  effect  might  then  provide 
an  additional  20  percent  reduction  in  altitude  loss  over  and  above  that  due 
to  mass  loss.  This  effect  is  not  available  for  the  latest  aerodynamic- 
configurations,  however. 

An  example  of  the  combined  effects  of  thrust,  mass  loss  and  cavity  is 
shown  for  the  8:1  cylinder  for  which  cavity  data  are  available  (Figures  35, 

36,  and  37). 

In  Figure  35,  the  computed  trajectory  for  a  "live"  thrusting  8:1  cylindrical 
flare,  assumed  to  have  negligible  precession  for  15  seconds  of  burning  time. 
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Figure  35.  Typical  Effect  of  Thrust,  Mass  Loss,  and  Cavity 
Aerodynamics  on  Nonprecessing  Flare  Trajectory 
(8:1  Cylinder) 
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Figure  36.  Typical  Effect  of  Thrust,  Mass  Loss,  and  Cavity  Aerodynamics 
on  Precessing  Flare  Trajectory  (8:1  Cylinder)  -  Position 
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Figure  37.  Typical  Effect  of  Thrust,  Mass  Loss,  and  Cavity  Aerodynamics 
on  Precessing  Flare  Trajectory  (8:1)  Cylinder)  -  Orientation 
Angles  .78. 


-  y - 

is  compared  to  the  corresponding  inert  flare  trajectory.  The  combined 
effect  is  an  almost  a  two-thirds  reduction  in  altitude  loss  with  0.  5  pound 
of  thrust.  A  higher  level  of  thrust  would  be  required  to  produce  a  climb 
above  the  launch  point,  but  the  trajectory  has  almost  leveled  off  at  burnout 
for  this  particular  case. 

In  Figures  36  and  37,  the  trajectory  parameters  for  the  same  "live"  and  in¬ 
ert  flares  are  compared  when  they  are  permitted  to  precess  under  the 
aerodynamic  torques  experienced.  As  may  be  inferred  from  the  inert 
response,  the  8:1  cylinder  itself  is  considered  too  thin  to  provide  the 
desired  precessional  characteristics  from  a  level  launch  so  that  conclusions 
based  on  this  simulated  "live"  response  would  be  somewhat  academic  and 
premature. 

In  the  absence  of  corresponding  cavity  effect  data  for  thicker  configurations 
and  the  significant  pressures  associated  with  thrusting,  it  is  recommended 
that  future  studies  with  thrust  assume  that  the  aerodynamics  are  the  same 
as  those  for  a  solid  flare  of  the  desired  dimensions.  Typical  coefficients 
for  such  a  solid  flare  and  other  configurations  studies  are  shown  in  Appendix 
D. 


(2)  Constant  Altitude  -  While  the  ballistic  coefficients,  burning 
time,  and  launch  and  ignition  conditions  define  the  altitude  lost  for  a  flare 
delivered  in  a  level  orientation,  a  combination  of  these  factors  plus  the 
proper  orientation  could  be  utilized  to  reduce  the  average  altitude  excursion. 
This  is  illustrated  in  Figure  38,  where,  with  the  same  aerodynamic  data 
and  nonprecessing  simulation  as  for  the  Figures  33  and  34,  it  was  indicated 
that  a  stable  3  inch  'liameter,  0.  5  pound  -  10  second  flare  could  be  kept 
within  ±100  feet  of  launch  altitude  at  net  flare  launch  forward  speeds  above 
about  250  feet  per  second.  Net  forward  speed  is  the  vector  sum  of  aircraft 
and  launcher  velocities.  Launcher  velocity  vector  and  aircraft  were 
horizontal,  but  flare  pitch  attitude  and,  consequently,  launch  angle  of 
attack  were  forward  edge-up  at  the  angles  indicated.  These  cursory 


-79- 


FLARE  DOWN  RANGE  TRAVEL  FROM  LAUNCH  >  FEET 


133J  -  soruinv  HONnvn  wow  tivj  3v\hj 


80 


Figure  38.  Locus  of  Burnout  Points  for  Typical  Constant  Altitude  Self- Suspended  Flare 


results  indicate  that  some  further  exploration  might  be  worthwhile.  They 
certainly  bring  out  the  influence  of  flare  orientation. 


(3)  Effect  of  Mass  Loss  on  Terminal  Velocity  -  While  the  inter¬ 
action  at  gliding  flight  conditions  requires  numerical  simulation,  consider¬ 
able  insight  into  the  boundaries  can  be  obtained  from  examination  of 
terminal  vertical  fall.  This  idealized  case,  approached  for  slow  launch 
speeds  or  long  ignition  delays,  has  yielded  closed  form  solutions  which 
could  permit  analytical  tradeoffs  between  flare  ballistic  and  pyrotechnic 
parameters.  This  analytical  solution  is  presented  in  graphical  examples 
and  general  equation  form  in  Figure  39. 


The  simple  analytical  solution  is  that  the  average  terminal  fall  velocity  for 
a  disk  of  variable  mass  but  constant  drag  coefficient  and  area  is  two- thirds 
the  terminal  velocity  at  the  start  of  burning;  (Vt&ve  =2/3  V  ).  Thus,  the 
total  fall  for  a  given  size  and  pyrotechnic  weight  varies  directly  with  the 
burning  times.  A  16  inch  disk  with  15  pounds  of  pyrotechnic  (0.  5  pound 
case)  is  predicted  to  fall  about  3100  feet  in  60  seconds. 


With  regard  to  tradeoffs,  these  curves  show  that  the  thir  ner  disk  is  better 
for  altitude  loss  of  a  level  flare  of  a  given  weight,  whereas  precessional 
stability  studies  indicate  that  the  thicker  the  better.  Therefore,  additional 
overriding  tradeoffs  would  appear  to  be  in  order.  A  thin  disk  would  pre¬ 
sumably  fall  faster  if  precessed  on  edge  than  a  thick  disk  only  partially 
tipped. 


f.  Typical  Equilibrium  Glide  Conditions  -  In  line  with  the  observation 
that  precession  can  be  prevented  only  by  gliding  at  an  angle  of  attack  at 
which  the  overturning  pitching  moment  is  zero  (a  =  aDp)  because  of  camber 
or  thickness,  the  desired  equilibrium  conditions  may  easily  be  computed 
from  the  aerodynamic  coefficients. 
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Figure  39.  Normalized  Fall  Distance  for  Self- Suspended  Flare 


Given,  a^,  find  the  glide  values  of  and  CA  at  this  angle  of  attack,  then 
the  glide  pitch  angle  0q,  and  glide  velocity,  Vq, 


aG  " 

aop 

(glide  angle  of  attack) 

*G  = 

tan'1  <-Ca/Cn>g 

(glide  pitch  angle) 

yG  = 

dG  -  aG 

(glide  path  angle) 

r  2W  cos  0Q  -j1/2 

VG  * 

L  P s  J 

(glide  velocity) 

or  more  familiarly  using  and 

yG  *  "tan  1  / ( c  l/ cd  )q 

2  ,  r  -2  W  sin  VG  -I 

VG  L  OS  Cdg  J 


VG  +  aG 
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Some  typical  values  for  these  parameters  based  on  the  aerodynamic  data 
are  shown  in  Table  1. 


The  problem  is  insuring  that  the  transient  between  launch  and  this  equlibrium 
condition  will  not  cause  precession.  With  the  current  estimates  of  very 
small  Magnus  moments  the  pitch  attitude  change  during  the  transient  should 
be  small,  therefore  the  glide  pitch  attitude  is  expected  to  define  the  required 
launch  attitude: 


0Q  9  Qq  (initial  pitch  altitude) 

<xq  =  0q  -  yq  (initial  angle  of  attack) 
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respectively. 


g.  Typical  Trajectories  for  Other  Applications  -  An  obvious  extension 
of  the  flying  spinning  disk  is  to  stabilize  delivery  of  high  explosive  ordnance.. 
One  form  could  be  launch  in  a  short  range  attack  or  defensive  mode  from  the 
surface.  Some  samples  have  been  computed  using  the  disk  aerodynamics  in 
the  nonprecessing  time- share  simulation.  Typical  results  for  a  4  inch. 

1  pound  disk  are  shown  in  Figures  40  and  41  as  a  function  of  launch  angle 
and  launch  velocity.  Spin  rate  was  not  specified  since  this  trajectory 
program  option  assumes  that  constant  attitude  is  obtainable.  Figure  41 
illustrates  the  effect  observed  throughout  the  study, and  capitalized  on  in 
the  steady  glide  analysis, that  the  disk  attitude  has  a  greater  influence  on 
the  terminal  flight  conditions  than  the  launch  flight  path  direction.  The  tick 
marks  on  the  curves  are  intended  to  denote  the  assumed  constant  orienta¬ 
tion  of  the  disk  at  those  points. 


2.  Stability  and  Precessional  Criteria  for  Spinning  Disks 


a.  Introduction  -  Simplified  analytical  solutions  have  been  obtained  for 
the  linearized  stability  and  precessional  response  equations  of  frisbee-like 
objects.  Correlation  with  local  and  Hurricane  Mesa  observations  is  con¬ 
sidered  quite  good  with  respect  to  the  following: 


Hand-launched  toy  frisbees,  which  some  persons  can  fly  well  for  a 
while. 

Light-weight  solid  inert  flares,  which  precessed  inversely  with  weight. 

Beveled  edge  solid  inert  flares,  which  precessed  much  worse  at  any 
weight. 
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DIMETER  -  4  INCHES 
THICKNESS  •  1  INCH 


DIAMETER  -  4  INCHES 
HICKNESS  ■  1  INCH 


Figure  41.  Typical  Trajectory  for  Flying  Disks  -  45  Degree  Launch  Elevation 


.  Clay  pigeons  thrown  from  the  cliff,  which  glided  straight  for  long 
periods  of  time. 

Correlation  with  downward  launches  at  Hurricane  Mesa  and  clay  pigeons 
at  the  trap  range  appears  to  be  dependent  on  the  magnitudes  for  the  non* 
linearities  of  the  parameters.  Qualitative  correlation  with  the  trajectory 
computer  program  is  considered  good  with  respect  to  the  shape  of  the 
curves  and  the  validity  of  the  nutation  (wobble)  suppression  assumption. 

The  preliminary  conclusions  that  may  be  drawn  regarding  design  criteria 
for  negligible,  or  acceptable,  precession  of  ordnance  items  utilizing  the 
frisbee  principle  may  be  extremely  difficult  to  meet. 

The  required  configuration  is  not  yet  apparent.  It  will  require  additional 
aerodynamic  data  to  define  whether  a  desired  or  practical  configuration  for 
a  given  application  will  provide  sufficient  orientation  stability. 

The  present  equations  and  available  data  indicate,  however,  that  the 
configuration  should  be  a  thick  disk  (or  short  cylinder)  with  relatively 
blunt  symmetrical  corners  for  arbitrary  launch  conditions,  or  a  domed/ 
cupped  configuration  on  the  order  of  a  solid  or  half- solid  clay  pigeon. 

b.  Stability  and  Precessional  Criteria  Defined  -  As  usually  defined, 
dynamic  stability  requires  that  angular  rates  and  linear  velocities  will 
return  to  equilibrium  when  disturbed,  with  some  exponential  envelope. 

This  does  not  necessarily  (or  even  usually)  mean  that  attitude  or  position 
will  return  to  their  original  values.  This  definition  of  dynamic  stability  is 
that  used  in  this  study.  The  precessional  criteria  referred  to  in  this  present 
report  are  especially  defined  to  represent  constraints  to  prevent  changes 
in  attitude  and  rates  under  forcing  functions  which  act  to  change  the  equilib¬ 
rium  conditions  for  a  dynamically  stable  body.  As  indicated  in  section  IV.  A.  1, 
a  precession  is  herein  defined  to  be  the  steady  gyroscopic  rotation  of  a 
spinning  body  under  the  applied  torques. 
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It  must  be  noted  that  both  the  stability  criteria  and  the  precession  criteria 
must  be  satisfied  to  obtain  straight  flight.  Satisfaction  of  the  stability 
criteria  is  the  necessary,  but  not  wholly  sufficient,  condition  that  small 
disturbances  will  not  grow.  However,  the  problem  is  that  the  response  of 
a  disk  that  is  dynamically  stable,  as  usually  defined,  to  an  initial  angle  of 
attack  (produced  by  pitch  from  horizontal  or  by  the  trajectory  slope,  to  an 
overturning  moment  due  to  asymmetry  (camber),  or  to  trajectory  curvature 
under  gravity)  is  to  steadily  precess  its  spin  axis  and  lift  vector  out  of  the 
vertical  plane  and,  thus,  to  turn  the  trajectory. 

Thus,  one  must  strive  to  minimize  the  response  to  these  forcing-functions 
in  addition  to  providing  stability  in  the  usual  sense.  The  most  obvious 
trivial  solution  for  arbitrary  forcing  functions  is  to  have  the  center  of 
pressure  on  the  center  of  gravity,  the  next,  to  provide  overwhelming 
gyroscopic  momentum,  but  they  may  or  may  not  be  obtainable. 

c.  Stability  and  Precessional  Criteria  for  Level  Launch  -  The  flare 
equations  of  motion  in  Appendix  B,  transformed  to  velocity  axes  which  are 
parallel  or  perpendicular  to  the  velocity  vector,  were  linearized  with  respect 
to  a  steady- state  planar  trajectory.  They  were  solved  for  stability  roots  of 
the  characteristic  equation  and  for  response  of  the  roll  angle  of  the  lift 
vector  and  of  the  angle  of  attack  to  initial  conditions,  trajectory  curvature, 
and  aerodynamic  asymmetries  by  use  of  the  La  Place  transform  technique. 

The  stability  roots  may  be  easily  separated  into  two  modes,  the  high 
frequency  nutation  and  the  low-frequency  undulation  modes.  The  high 
frequency  nutation  mode,  corresponding  to  the  nutation  of  projectiles  and 
Magnus  rotors  is  always  oscillatory  if  gyroscopic  stabilization  has  been 
satisfied.  The  undulation  mode,  however,  usually  degenerates  into  two 
aperioidc  submodes,  one  of  which  could  still  be  called  the  undulation  mode 
and  the  other  the  precession  mode.  Undulation  is  dominant  in  magnitude 
of  the  root,  while  the  precessional  root  approaches  zero,  or  pure 
precession. 
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The  dynamic  stability  criteria  for  these  modes  represent  the  inequality  con¬ 
straints  that  must  be  met  to  insure  that  the  stability  roots  are  at  least 
neutral  (i.  e. ,  zero),  or  preferably  damped  (i.  e. ,  negative).  The  nutation 
damping  criterion,  the  nutation  frequency,  and  the  undulation  damping 
criterion  obtained  turn  out  to  be  identical  to  those  for  a  Magnus  rotor  in  a 
steady  glide  with  spin  axis  horizontal.  The  criterion  for  gyroscopic  stability 
is  identical  to  that  for  a  spinning  projectile. 

The  amplitude  of  the  nutational  mode  response  of  a  stable  prespun  disk  to 
the  forcing  functions  of  concern  was  negligible  and,  therefore,  was  suppressed 
in  the  results  and  in  the  computer  program.  The  damping  exponent  of  the 
undulation  submode  dominates  the  response. 

For  arbitrary  elevation  of  the  initial  velocity  vector  the  linearized 
equations  yielded  undulation-precession  roots  which  appeared  to  be  quite 
slow  convergences  or  divergences  depending  on  the  magnitudes  and  signs 
of  the  trajectory  parameters.  This  variation  makes  general  analytical 
solutions  difficult  to  obtain  and  interpret,  let  alone  to  describe  and/or  to 
present  criteria. 

The  troublesome  terms  are  zero  or  negligible  for  level  launch  flight  paths, 
however,  so  that  quite  simple  expressions  were  obtained  for  the  roll  rate  in 
a  region  near  enough  to  a  level  launch  flight  path  for  linearity  to  apply. 

These  level  launch  equations  were  simplified  into  the  combined  precession 
constraint  equations  for  initial  roll  and  pitch  and  steady-  state  angle  of 
attack  and  velocity  of  Table  2.  If  these  constraints  are  satisfied  simultane¬ 
ously,  steady- state  roll  rate  will  be  zero  and  the  roll  attitude  will  be  less 
than  some  arbitrary  maximum  value. 

Basically,  the  simultaneous  satisfactioa  which  could  constitute  a  juggling 
act,  requires  that  the  velocity  be  such  after  a  short  transient,  the  lift  will 
equal  the  weight  at  the  angle  of  attack  for  neutral  overturning  moment  and 
the  initial  pitch  and  roll  attitude  be  such  that  after  the  transient,  residual 
roll  is  some  acceptable  low  value  and  the  velocity  will  be  constant. 


Table  2.  Linear  Criteria  for  Negligible  Precession  from 
Level  Launch 
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These  criteria  were  derived  from  the  more  general  steady- state  roll-pitch 
equations  in  Table  3  which  are,  in  turn,  based  on  the  undulation/precession 
response  equations  and  lumped  coefficient  definitions  of  Table  4.  The 
glossary  for  these  tables  is  presented  in  Table  5. 

Table  6  shows  the  dynamic  stability  criteria  applicable  to  the  level  launch 
case.  These  criteria  for  undulation  and  nutation  mode  dynamic  stability  and 
the  necessary  gyroscopic  stability  also  apply  to  other  trajectory  conditions 
but  would  in  those  cases  be  supplemented  by  additional  trajectory-dependent 
criteria  which  are  not  available  for  presentation.  The  stability  an  i  response 
equations  should  be  extender  :o  these  other  conditions  at  the  earliest 
opportunity. 

Some  numerical  examples  of  equilbrium  glide  trajectory  parameters  based 
on  the  exact  equations  and  the  nonlinear  aerodynamic  data  have  been  shown 
in  Table  1  for  typical  configurations.  The  corresponding  initial  roll  orienta¬ 
tions  required  by  the  linear  criteria  of  Table  2  are  also  shown  in  Table  1 
under  the  simplifying  assumptions  that  launch  occurred  near  the  equilibrium 
glide  velocity  and  that  the  peripheral  speed  ratio  was  approximately  0.  8  for 
all  cases.  Naturally,  where  the  roll  angle  predicted  by  linear  theory  and  the 
velocity  change  exceeds  the  usual  bounds  of  linear  theory,  the  required 
launch  conditions  in  both  roll  and  pitch  would  require  an  iterative  numerical 
solution  of  the  nonlinear  equations. 

d.  Application  of  the  Criteria  to  Friabees  and  Clay  Pigeons  -  Satisfaction 
of  the  precessional  criteria  represents  simultaneous  satisfaction  of  a 
relatively  precarious  simultaneous  conditional  equilibrium  condition  or 
coming  close  enough  that  the  error  from  equilibrium  causes  an  acceptably 
slow  divergence  from  a  straight  flight.  This  is  apparently  where  athletic 
skill  comes  into  the  picture  for  frisbees.  However,  the  aerodynamic  and 
inertia  characteristics  of  the  frisbee  make  it  relatively  easy  to  satisfy 
these  simultaneous  requirements  for  a  frisbee  or  clay  pigeon  and  difficult 
for  a  heavy  disk. 
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Table  3.  Linear  Criteria  for  Steady-State  Precession  from  Level  Launch 
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Table  4.  Quasi-Steady  Response  Equations  for  Level  Launch 
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Table  5.  Glossary  for  Stability/Precession  Criteria 
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Table  6.  Linear  Criteria  for  Dynamic  Stability  for 
Level  Launch 
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A  key  feature,  as  could  be  expected,  is  the  curvature  or  camber  of  the  thin 
disk.  Thin  flat  disks  are  aerodynamically  unstable  edge- on  and  the  natural 
tendency  is  to  overturn  to  the  broadside  orientation.  The  curvature  causes 
the  moment  to  be  a  tuck-under  pitching  moment  at  low  angles  of  attack  and 
overturning  at  higher  angle  of  attack,  with  equilibrium  at  a  finite  angle  of 
attack.  This  equilibrium  angle  for  zero  pitching  moment  ( aQp  in  the 
equations)  is  on  the  order  of  10  degrees. 

The  disk  is  still  aerodynamically  unstable  about  this  shifted  equilibrium 
point  so  that  gyroscopic  momentum  due  to  spin  is  needed  for  stability. 
Available  data  do  not  show  significant  effect  of  the  spin  on  the  aerodynamics 
themselves.  Under  the  right-hand  rule,  the  gyroscopic  reaction  at  low 
angles  of  attack  would  be  a  roll  to  the  right  This  helps  explain  the  rollout 
of  a  frisbee  tossed  on-edge.  At  higher  angles,  the  gyroscopic  reaction 
would  be  a  roll  to  the  left.  This  shift  helps  permit  out- smarting  the 
gyroscopic  right  hand  rule. 

For  the  available  data,  this  roll  is  ±1  degree  if  the  initial  angle  of  attack 
were  ±4  degrees  from  the  neutral  value.  Because  of  the  light  weigh*  of  the 
frisbee,  it  is  possible  to  make  the  lift  equal  the  weight  at  the  angle  of  attack 
for  zero  pitching  moment.  If  the  lift  equals  the  weight,  the  trajectory 
curvature  will  be  zero  and  the  angle  of  attack  will  not  change  from  this 
equilibrium.  If  launched  at  a  lower  angle  of  attack,  the  frisbee  will  drop 
until  the  lift  equals  the  weight.  If  the  proper  initial  tilt  to  the  left  is  used, 
some  of  which  occurs  naturally  from  the  weight,  the  roll  to  the  right  during 
the  transient  will  bring  it  up  level.  If  thrown  at  a  higher  angle  of  attack, 
the  flight  path  will  climb  until  satisfied  and  it  would  roll  to  the  left.  If  the 
angle  of  attack  when  lift  equals  weight  is  the  rotational  equilibrium  angle, 
the  transition  right  or  left  roll  will  stop. 

The  velocity  at  which  the  lift  equals  the  weight  at  the  equilibrium  angle  of 
attack  is  between  20  -  25  feet  per  second  for  the  available  wind  tunnel  data 
and  corr  y^onds  to  noninstrumented  observations  of  actual  flights.  After 
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the  initial  transient  to  this  steady-state  condition  as  the  velocity  decays  due 
to  air  drag,  the  path  will  curve  downward,  the  angle  of  attack  will  depart 
from  the  equilibrium  point  and  peeloff  occurs.  It  should  be  apparent  that 
if  the  toss  is  such  that  the  precession  criteria  are  met  on  a  downward  glide-  * 
path  which  will  also  maintain  constant  velocity  (the  tangential  gravity 
component  equals  drag)  the  peeloff  will  be  delayed  until  spin  decay  or  other  f 
disturbances  (such  as  hitting  the  ground  or  a  wall)  intervene. 

For  long  term  straight  flight,  the  combination  of  initial  orientation  and 
velocity  has  to  be  such  that  when  the  lift/ weight  force  equilibrium  is  reached 
(as  quickly  as  possible)  the  angle  of  attack  is  that  for  rotational  equilibrium, 
the  roll  orientation  is  negligible,  and  the  velocity  is  in  equilibrium. 

For  the  frisbee  on  level  ground,  the  velocity  equilibrium  is  not  too  critical 
because  the  decay  is  slow  relative  to  the  desired  flight  time.  For  a  clay 
pigeon  on  a  trap  range,  the  time  is  also  apparently  too  short  for  significant 
observable  peeloff  to  occur  since  there  is  definite  trajectory  curvature. 

The  principle  was  validated  with  clay  pigeons  tossed  with  edges  slightly 
downward  off  the  cliff  at  Hurricane  Mesa.  They  quickly  reached  stable 
equilibrium  in  both  attitude  and  velocity  and  flew  straight  with  wings  level 
for  significant  times,  changing  direction  primarily  with  wind  shifts. 

e.  Application  of  Criteria  to  Self-Suspended  Flares  -  It  is  relatively 
easy  to  meet  the  precessional  criteria  with  a  fristee  because  the  light 
weight  makes  it  easy  to  rapidly  obtain  equilibrium  between  lift  and  weight 
and  the  flight  time  prevents  visible  buildup  of  near  misses  on  orientation. 

For  an  air-launched  heavy  disk  with  finite  flight  time,  achievement  of  the 
equilibrium  is  more  difficult  because  of  the  possible  wider  difference 
between  launch  and  equilibrium  and  a  possible  discrepancy  between  lift 
equilibrium  and  moment  equilibrium.  The  transition  phase  could  also 
cause  excessive  roll  excursion  even  if  a  satisfactory  equilibrium  condition 
exists. 
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The  more  solid  and/or  symmetrical,  the  smaller  the  angle  for  rotational 
equilibrium,  but  the  heavier,  the  larger  the  angle  required  for  lift 
equilibrium. 

As  may  be  shown,  the  precessional  criteria  involve  a  combination  of: 
pitching  moment  characteristics,  gyroscopic  stabilization  through  spin, 
and  flight  conditions.  Thus,  alternate  approaches  could  be  pursued. 

(1)  Brute  Force  Mode  -  If  gyroscopics  can  be  used  to  overwhelm 
aerodynamic  characteristics,  then  flight  conditions  would  not  be  critical. 
This  brute  force  approach  would  involve  development  of  a  shape  whose 
aerodynamic  moments  are  negligible  over  the  angle  of  attack  region 
anticipated  with  arbitrary  launch  attitude  compared  to  the  gyroscopic 
momentum  achievable  with  practical  spin  rates  and  mass  distributions. 
Thicker  cylinders  approaching  2:1  or  1:1  diameter- to- thickness  ratios 
appear  to  be  a  possible  means  for  achieving  neutral  pitching  moments, 

at  least  at  small  angles  of  attack.  Our  measurements  of  the  4:1  cylinder 
showed  a  trend  in  this  direction.  There  would  have  to  be  careful  juggling 
to  insure  that  the  Magnus  moments  that  are  apt  to  be  generated  by  the  edges 
of  the  thicker  disks  would  not  be  detrimental. 

(2)  Frisbee  Mode  -  The  other  alternate  would,  of  course,  be  to 
provide  a  configuration  and  delivery  scheme  completely  analogous  to  frisbee 
delivery.  Thus,  the  configuration  would  be  domed  and  preferably  cupped 
and  the  launch  orientation  would  be  near  that  nose  down  attitude  required 
for  equilibrium  glide  at  the  rotational  equilibrium  condition,  rolled  to 
compensate  for  the  transition  from  launch  (see  Table  1). 

Another  interesting  approach  represents  a  combination  of  the  proceding 
approaches  to  the  aerodynamic  and  launch  to  glide  transition  problem. 
Rather  than  a  shift  of  otherwise  linear  moments  through  camber  (doming/ 
cupping),  a  symmetrical  cylinder  which  has  a  neutral  pitching  moment  over 
a  significant  ±  range  sufficient  to  provide  lift  equilibrium  might  be  launched 
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in  the  glide  orientation  with  negligible  precession  during  transition  and  the 
angle  would  not  be  critical  as  long  as  it  was  within  the  neutral  zone. 

The  4:1  cylinder  data  approaches  this,  but  the  deviations  from  neutrality 
appear  to  still  be  too  great  and  the  angular  region  of  near  neutrality  is 
not  quite  great  enough  for  expected  launch  condition. 

3.  Configuration  Aerodynamic  Characteristics 

The  configuration  synthesis  aerodynamic  definition  process  constituted  an 
iterative  procedure  which  is  not  yet  complete.  These  studies  started  with 
identification  of  the  three  basic  configurations  based  on  inputs  from  the 
customer  and  the  origins  of  the  concept.  The  tools  used  included  literature 
searches  for  available  data,  estimates  based  on  available  data,  nonspinning 
and  spinning  wind  tunnel  tests,  flow  theory,  observation  of  flight  tests, 
computer  simulation  of  trajectories  and  all  the  other  "tricks"  of  the 
aerodynamicists  trade,  including  omphaloskepsis  and  engineering  judgment. 
From  this  evolved  prediction  of  aerodynamic  characteristics  for  configura¬ 
tions  tested  and  new  generations  to  be  tested. 

The  basic  configurations  used  were  an  8:1  diameter-to-thickness  ratio  right 
circular  cylinder  with  a  small  radius  on  the  top  corner,  representing  a 
16-inch- diameter,  2-inch- thick,  60- second  burning  flare  that  was  an  inter¬ 
im  design  goal,  with  four  levels  of  flare  cavity,  and  called  the  baseline 
flare;  a  frisbee-Iike  configuration  with  two  cavity  levels,  hollow  and  solid; 
and  a  solid  and  \\iollow  clay  pigeon.  To  this  was  later  added  a  solid  4:1 
cylinder,  with  a  small  radius  on  the  top  corner,  which  corresponds  to  the 
8  inch  diameter,  2-inch*thick  live  flare  launched  at  Hurricane  Mesa  The 
data  obtained  on  these  slopes  were  interpolated  for  8:1,  6:1,  and  4:1 
symmetrical  and  beveled  inert  flares  and  extrapolated  for  thicker  bodies. 

a.  Configuration  Aerodynamic  Data  -  The  results  of  the  wind  tunnel 
tests  of  the  hollow  models  of  the  three  basic  configurations  are  typified 
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by  the  pitching  moment  coefficient  plot  of  Figure  42.  These  domed  clay 
pigeon  and  frisbee  configurations  display  similar  linear  increases  until  a 
stall  point  or  flow  separation  is  reached  and  display  a  similar  shift  of  the 
point  of  neutral  aerodynamic  moment  to  12  -  13  degrees  of  angle  of  attack  due 
the  asymmetry.  This  effect  could  be  expected  from  the  theory  of  cambered 
thin  airfoils.  A  similar  rise  to  a  stall  point  is  displayed  by  the  normal 
force  coefficients,  C^j,  in  Figure  43.  The  baseline  cylindrical  flare 
virtually  eliminates  the  linear  region  and  sharp  stall  at  finite  angle  of 
attack  in  both  these  figures  because  the  bluff  contours  undoubtedly  cause 
earlier  flow  separation  and  early  stall.  The  corresponding  axial  force 
coefficient,  C^,  parallel  to  the  disk  is  shown  in  Figure  44.  These  coeffi¬ 
cients  were  defined  previously  in  section  IV.  A.  3. 

Similar  trends  are  shown  in  the  results  for  solid  shapes  with  varying  degrees 
of  contouring  or  thickness  corresponding  to  launch  configurations  of  the 
flare  prior  to  ignition.  The  pitching  moment  coefficient  is  shown  in  Figure 
45,  the  normal  force  coefficient  in  Figure  46,  and  the  axial  force  coefficient 
in  Figure  47. 

These  data  show  both  the  original  8:1  diameter-to-thickness  ratio  baseline 
cylinder  and  the  4:1  cylinder  tested  on  the  spin  fixture.  In  both  cases, 
the  cylinders  had  an  approximately  0.  1  inch  (full-scale)  upper  corner 
radius  and  square  lower  corner  similar  to  the  prototype  live  flares.  This 
introduces  a  slight  asymmetry  into  the  data  at  small  angles  of  attack,  which 
was  suppressed  for  computer  simulation  of  the  symmetrical  dummy  flares 
launched  at  Hurricane  Mesa. 

Figure  45  gives  a  combined  indication  of  the  effects  of  thickness  and  corner 
radius  on  the  pitching  moment  coefficient  which  appear  to  dominate  in  the 
region  below  45  degrees  angle  of  attack.  Aside  from  the  aforementioned 
suppression  of  stall,  there  is  a  much  less  pronounced  effect  on  normal  force 
(Figure  46).  As  a  matter  of  fact,  there  is  not  much  departure  from  the 
often  used  crude  sine  wave  first  approximation  corresponding  to  the  normal 
component  of  drag  due  to  uniform  pressure  on  an  inclined  plate. 
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PITCHING  MOHENT  COEFFICIENT 
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Figure  42.  Pitching  Moment  Coefficients  for  Typical  Hollow  Self- 
Suspended  Flare  Shapes  (Nonspinning) 
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NORMAL  FORCE  COEFFICIENT 
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Figure  43.  Normal  Force  Coefficients  for  Typical  Hollow  Self- Suspended 
Flare  Shapes  (Nonspinning) 
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AXIAL  FORCE  COEFFICIENT 
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Figure  44.  Axial  Force  Coefficients  for  Typical  Hollow  Self 
Suspended  Flare  Shapes  (Nonspinning) 
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pitching  moment  coefficient 


Figure  45.  Pitching  Moment  Coefficients  for  Typical  Solid  Self- Suspended 
Flare  Shapes  (Nonspinning) 
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Figure  46. 


Normal  Force  Coefficients  for  Typical  Solid  Self- Sura*  ad  *d 
Flare  Shapes  (Nonspinning) 
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b.  Thickness  Effects  -  Examination  of  the  solid  configuration  data  and 
the  available  literature  leads  to  the  impression  of  a  trend  towards  reduction 
of  the  overturning  moment  at  small  angles  of  attack  as  the  effective  leading 
edge  thickness  increases.  We  know  from  Magnus  rotor  work  that  longer 
cylinders,  with  a  diameter-to-thickness  ratio  of  approximately  1:2,  have 
relatively  strong  restoring  moments.  The  crossover  between  overturning 
and  restoring  is  near  1:1.  In  our  case,  not  only  is  the  slope  at  small 
angles  of  concern  but  the  average  level  over  a  ±  30-40  degree  region. 

In  order  to  obtain  gliding  flight,  if  that  mode  is  desired,  the  drag  should 
not  be  much  higher  than  the  lift.  None  of  the  thicknesses  considered  seem 
to  be  neutral  enough  over  a  large  enough  angle  of  attack  range  to  permit 
brute-force  delivery  in  a  horizontal  orientation.  The  thicker  the  disk,  the 
more  neutral  the  pitching  moment  with  angle  of  attack,  but  the  thicker  the 
disk,  the  steeper  the  trajectory  and  the  higher  the  angle  of  attack.  Thus, 
it  would  appear  that  a  3:1  or  possibly  2:1  diameter-to-thickness  ratio  is 
likely  to  provide  a  desirable  combination  of  stability  and  trajectory  charac¬ 
teristics  for  delivery  as  a  disk  when  coupled  with  the  ideal  delivery 
orientation,  leading  edge  down. 

c.  Corner  Radius  Effects  -  The  detailed  effect  of  small  corner  radii 
on  basically  cylindrical  disks  has  not  yet  been  completely  isolated  because 
the  radius  used  on  only  the  top  corner  of  the  bodies  tested  introduced  a 
slight  asymmetry  into  the  data.  This  was  arbitrarily  faired  through  ze^o 
for  the  symmetrical  dummy  flares  which  had  a  bevel  radius  top  and  bottom. 

The  effect  of\  gross  corner  radius  depends  on  whether  it  is  symmetrical 
or  on  the  top  surface  only.  On  the  basis  of  the  frisbee,  clay  pigeoa  discus, 
(and  lentricular  body  data  obtained  since  completion  of  the  work  and  analysis,  ) 
it  would  appear  to  provide  a  region  where  the  normal  force  and  pitching 
moment  are  linear  with  angle  of  attack.  The  extent  of  this  region,  or  the 
stall  point  where  flow  separates  from  the  top  surface,  thus  appears  to  be  a 
function  of  this  radius  and  the  effective  leading  edge  thickness. 
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The  frisbee  and  clay  pigeon  had  similar  stall  points  but  different  basic 
thicknesses  and  different  apparent  corner  radii,  but  similar  leading  edge 
thicknesses.  Both  of  these  were  approximately  half  that  of  the  8:1  cylinder. 
While  the  slopes  of  the  curves  were  all  similar,  and  the  stall  point  and 
depth  of  drop  after  the  stall,  if  any,  varied  with  corner  radius,  the  dominant 
effect  observed  was  the  shift  of  the  point  of  zero  moment  for  the  non- 
symmetrical  bodies.  An  increase  in  effective  camber,  by  imposing  an 
ever  increasing  nose-down  couple  on  the  linearly  increasing  overturning 
moment  due  to  the  lift  force,  increases  the  angle  of  attack  at  which  the 
moments  balance.  This  could  be  true  only  up  to  a  point  since  the  clay 
pigeon  has  more  camber  than  the  frisbee  and  a  higher  equilibrium  angle, 
but  the  cylinder's  abrupt  change  in  camber  gives  a  different,  nonlinear 
effect  from  this  idealized  airfoil  model. 

The  optimum  corner  radius  has  yet  to  be  identified,  and  as  this  discussion 
should  have  implied,  would  depend  on  whether  one  wants  to  obtain  a  maximum 
nonlinear  region  of  neutral  moment  with  a  symmetrical  body  or  desires  to 
exploit  the  more  precarious  discrete  neutral  moment  point  with  its  associated 
finite  moment  at  zero  angle  of  attack,  of  the  unsymmetrical  configurations. 

A  proposed  future  wind  tunnel  study  to  explore  the  effects  of  corner  radii 
and  their  interaction  with  thickness  is  described  in  section  IV.  B.  3.  j. 

d.  Spin  Effects  -  The  spinning  model  wind  tunnel  tests  tended  to  confirm 
the  initial  impression  that  the  pitching  moment  and  lift  would  not  themselves 
be  affected  by  the  spin,  but  that  the  predominant  effect  would  be  the  pro¬ 
vision  of  gyroscopic  stiffness,  with  small  or  insignificant  Magnus  rolling 
moment. 

The  test  data  are  cross  plotted  in  Figure  48  for  the  pitching  moment  coeffi¬ 
cient  and  rolling  moment  versus  the  dimensionless  spin  rate  scaling  param¬ 
eter,  tud/2V,  (the  peripheral  speed  ratio  Vr/V)  at  typical  values  of  angle  of 
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Figure  48.  Effect  of  Peripheral  Speed  Ratio  on  the  Aerodynamic 
Characteristics  of  a  Right  Circular  Cylinder 
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attack.  The  variations  in  the  data  up  to  the  maximum  peripheral  speed 
ratio  attainable  with  the  present  turbine  were  within  the  scatter  on  the 
data  itself.  It  was  concluded  that  the  effect  was  negligible. 

This  Vr/V  of  0.  6  corresponds  to  spin  rate  of  1800  rpm  at  100  feet  per 
second.  3600  rpm  at  200  feet  per  second,  and  9000  rpm  at  500  feet  per 
second  for  an  8  inch  flare.  Larger  flares  or  slower  speeds  would  reach 
this  0.  6  value  at  lower  velocities. 

While  the  rolling  moment,  C^,  at  first  glance  appears  to  show  some  effect 
of  spin,  no  quantitative  use  can  be  made  of  the  data  for  several  reasons. 

It  is  well  within  the  scatter  for  the  balance  used  on  the  expanded  scale 
plotted  in  Figure  48  and  is  completely  overwhelmed  by  the  transfer  of  the 
normal  force  acting  near  the  center  to  the  laterally  offset  roll  balance 
strain  gage  location.  This,  in  turn,  necessitated  a  use  of  a  balance  with  a 
high  full  scale  capacity  for  the  roll  gage,  thus  reducing  its  sensitivity  to 
the  minute  roll  torques. 

While  not  providing  usable  plots  of  magnus  moment  as  a  function  of  angle  of 
attack  and  spin  rate,  the  results  should  bracket  it  as  being  less  than  the 
tolerance  on  the  roll  gage.  1/4  percent  of  full  scale  of  ±  25  inch-pounds. 

In  terms  of  the  computer  program  which  currently  computes  the  Magnus 
moment  coefficient  as  a  constant  fraction  of  the  pitching  moment,  this  is 
less  than  1/8  percent  of  full  scale  pitching  moment  of  ±  50  inch-pounds. 

Ratios  of  under  10  percent  at  peripheral  speed  ratios  of  1.  0,  (C.  =  -.  1  C^) 

did  not  seem  to  significantly  effect  the  computed  trajectories.  Ten  percent 
might  be  visualized  physically  as  rotating  the  normal  force  center  of 
pressure  5.  73  degrees  to  the  right  of  the  aerodynamic  plane.  This  requires 
a  considerably  more  significant  effect  on  the  force  picture  than  normally 
observed  on  projectiles  or  predicted  theoretically  during  initial  stages  of 
this  project.  Therefore,  the  results  appear  to  have  been  adequate  for  the 
work  to  date. 


-Ill- 


It  is  also  recognized  from  our  Magnus  rotor  work  that  finite  length  cylinders 
(1:1,  1:2)  have  finite  Magnus  moments  with  Scale  factors  on  the  order  of  50 
percent  of  the  restoring  (pitching)  moments  for  the  1:2,  and  different  non- 
linearities.  For  one  with  neutral  pitching  moment,  the  scale  factor  would 
approach  infinity.  Therefore,  the  Magnus  moment  should  receive  more 
serious  attention  if  thicker  disks  are  studied  in  the  future. 

The  balance  had  to  be  offset  in  order  to  put  the  pitch  measurements  center 
on  the  pitch  axis  to  provide  the  sensitivity  needed  to  measure  any  differences 
in  pitching  moment  due  to  spin,  while  still  allowing  room  for  the  bearings 
and  shaft  and  spin  turbine  (see  Figure  19). 

The  only  solutions  apparent  would  be  a  separate  balance  on  the  centered 
roll  axis  with  drastically  reduced  scale  or  a  specially  built  Z  shaped 
balance  with  the  pitch  strain  gages  at  the  pitch  axis  as  now  and  the  roll 
strain  gages  on  the  centered  roll  axis.  The  existing  and  available  standard 
balances  did  not  have  the  required  ranges. 

Existing  Magnus  balances  at  other  facilities  were  designed  for  projectiles 
with  a  different  orientation  of  sensitive  axes.  The  problem  in  all  such 
balances  would  be  even  worse  here  because  of  the  large  normal  forces  on 
such  a  disk.  That  common  problem  is  providing  sufficient  stiffness  in  one 
axis  to  withstand  the  high  loads  and  sufficient  weakness  in  the  other  to 
permit  sensing  the  Magnus  loads  while  preventing  cross- coupling  of  the 
high  loads  into  the  more  sensitive  sensor  and  structure.  It  should  be 
obvious  that  such  a  special  balance  would  be  expensive,  but  it  should  be 
tried  someday  if  sufficient  funds  are  available. 

Development  of  hardware  and  software  to  extract  the  Magnus  data  from 
3-degree-of-freedom  dynamic  test  data  could  also  be  complex  and 
expensive,  since  this  body  has  a  different  axis  orientation,  and  inter¬ 
ference  effects  from  the  support  could  be  significant. 
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e.  Combustion  Effects  -  A  preliminary  analysis  of  the  effect  of  the 
combustion  wake  on  the  aerodynamic  characteristics  was  performed,  based 
on  preliminary  estimates  of  burning  pressures  and  temperatures  obtained 
from  NAD/Crane  prior  to  the  manifestation  of  sufficient  pressure  to 
generate  finite  thrust 

This  simplified  analysis  is  presented  in  Appendix  A. 

The  conclusion  was  reached  that  the  effect  of  the  flare  plume  on  the  aero¬ 
dynamic  characteristics  will  be  negligible  if  thrust  is  negligible.  Extension 
of  this  cursory  analysis  to  more  adequately  account  for  the  actual  combus¬ 
tion  behavior  would  in  itself  appear  to  be  a  major  research  task  for  a 
combustion  research  laboratory,  and  outside  the  scope  of  anticipated  near- 
term  developmental  efforts. 

If  the  actual  combustion  effects  do  place  the  pressures /forces  outside  the 
state  of  the  art  of  analytical  prediction,  then  the  only  recourse  is  to 
relatively  arbitrary  assumptions  for  analytical  flight  mechanics  studies 
and  to  free-flight  testing  for  definition  or  verification  of  the  required 
assumptions.  As  indicated  in  the  next  subsection  on  cavity  effects,  a 
recommended  interim  arbitrary  assumption  would  be  that  the  aerodynamic 
forces  and  moments  are  the  same  as  the  corresponding  solid  body,  with 
thrust  added  as  indicated  by  the  burning  characteristics  or  laboratory 
measurements. 

f.  Cavity  Effects  -*  The  measured  effects  of  the  variable  cavity  in  the 
8:1  baseline  cylindrical  flare  wind  tunnel  model  are  shown  in  Figure  49 
for  the  pitching  moment.  The  variation  here  is  very  erratic  and  did  not 
lend  itself  to  a  specific  or  generalized  model  for  interpolation.  It  was 
noted,  of  course,  that  very  little  change  occurred  until  the  cavity  radius 
exceeded  65  percent  of  the  body  radius.  Then  the  most  significant 
variation  occurred  between  this  condition  and  the  hollow  condition.  This 
apparent  insensitivity  to  cavity  until  virtually  hollow  was  even  more  pro¬ 
nounced  on  the  normal  force  and  axial  force  coefficients  (Figures  50  and  51). 


-113- 


ANGLE  OF  ATTACK, 0*  -  DEGREES 


Figure  49.  Effect  of  Cavity  on  Pitching  Moment  Coefficient  of  Baseline 
Self- Suspended  Flare  Shape  (Nonspinning) 
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Figure  50.  Effect  of  Cavity  on  Normal  Force  Coefficient  of  Baseline 
Self- Suspended  Flare  Shape  (Nonspinning) 
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Figure  51.  Effect  of  Cavity  on  Axial  Force  Coefficient  of  Baseline 
Self- Suspended  Flare  Shape  (Nonspinning) 


Even  though  the  overturning  moment  for  this  8:1  cylinder,  and  the  anticipated 
precessioa  was  reduced  compared  to  the  clay  pigeon  and  frisbee  it  still 
appeared  to  be  excessive  in  computer  simulations. 

The  testing  and  study  has  more  recently  concentrated  on  thicker  solid 
configurations  in  a  search  for  a  configuration  which  would  be  stable  during 
the  indefinite  period  between  launch  and  ignition,  before  extending  the  study 
of  cavity  burning  effects.  There  is  not  yet  verification  that  such  a  con¬ 
figuration  and  associated  delivery  technique  has  been  achieved,  and  the  burn¬ 
ing  pressures  are  presently  higher  than  anticipated  and  desired.  Therefore, 
no  attempt  has  been  made  to  estimate  the  effects  of  the  burning  cavity  on  the 
4:1  cylinder  data  or  on  the  6:1  dt.ta  interpolated  between  the  solid  8:1  and 
4:1  data  nor  to  extrapolate  for  thicker  configurations. 

The  available  data  on  cavity  effects  did  not  reveal  discernible  trends  upon 
which  to  base  a  standardized  interpolition/prediction  model.  Therefore, 
the  computer  program  is  configured  to  accept  the  primary  aero  coefficients 
as  measured  or  estimated  curves  for  four  different  cavity  radii  or  to 
ignore  the  cavity  change  completely  by  accepting  only  one  curve.  In  view 
of  the  present  and  anticipated  state  of  knowledge  of  the  high  pressure  plume 
effects  and  the  fact  that  the  8:1  cylindrical  configuration  for  which  cavity 
data  are  available  is  not  presently  recommended  for  tactical  application, 
the  recommended  interim  approach  would  be  to  use  one  curve  and  ignore 
the  changes  during  burning. 

As  could  be  observed  by  comparing  Figures  42  through  44  and  45  through 
47,  the  basic  form  of  the  coefficients  for  the  clay  pigeon  and  frisbee  is 
similar  whether  in  the  normal  hollow  or  the  flare /ordnance  oriented  solid 
condition.  In  both  cases,  filling  in  the  bottom  raises  the  level  and  slope 
of  the  overturning  moment  and  reduces  the  angles  for  neutral  overturning 
moment.  It  slightly  increases  the  normal  force  level  of  the  clay  pigeon 
while  decreasing  that  of  the  frisbee.  While  no  data  were  obtained  for  inter¬ 
mediate  cavities,  the  implication  is  that  a  partially  hollow  base  may  be 
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beneficial  in  tailoring  the  pitching  moment  asymmetry.  Such  a  lip,  if 
shallow  enough,  could  be  designed  for  nesting  so  as  to  maximize  the 
packing  density  while  filling  the  center  dome  with  payload  material.  A 
plan  for  exploring  partially  cupped  configurations  in  the  wind  tunnel  is 
discussed  in  section  IV.  B.  3.  j. 

g.  Aerodynamic  Damping  -  The  aerodynamic  pitch  damping  measured 
in  the  wind  tunnel  for  the  solid  and  hollow  frisbee  and  8:1  cylinder  con¬ 
figuration  is  shown  in  Figure  52.  The  tests  were  conducted  with  non¬ 
spinning  models  on  a  spring  restrained  side-mount.  The  interpretation 
made  by  the  wind  tunnel  personnel  was  that  the  damping  coefficients  were 
very  small  based  on  their  experience  with  conventional  and  unconventional 
bodies.  The  presence  of  some  positive  values  which  were  also  observed 
visually  as  amplitude  buildups  may  be  attributed  to  hysterisis-like  flow 
separation,  which  has  been  observed  for  thin  and/or  sharp  edged  disks 
and  is  the  cause  of  autorotation  of  thin  nonspinning  disks  and  plates.  This 
effect  [which  would  have  undamped  the  high  frequency  nutations  (wobble)] 
was  not  preceptible  on  any  of  the  flight  tests  conducted  to  date.  As  a 
matter  of  fact,  the  small  wobble  induced  by  the  launches  appeared  to  damp 
out  fairly  rapidly  even  though  quantitative  data  on  decay  were  not  obtained. 
Thus,  the  spin  either  prevents  this  phenomenon  from  occurring  aero- 
dynamically  or  resists  it  gyroscopically. 

The  results  of  the  trajectory  program  were  insensitive  to  the  magnitude  of 
the  aerodynamic  damping. 

h.  Idealized  Aerodynamic  Requirements /Alternate  Configurations  -  The 
computer  simulations  and  linearized  analysis  indicate  that  delivery  in  a 
manner  somewhat  analogous  to  that  of  the  frisbee  could  require  launching 
at  a  spin  stabilized  orientation  corresponding  to  that  which  it  would  have  in 
steady  glide,  front  edge  15-20  degrees  downward.  The  glide  angle  of  attack 
would  be  on  the  order  of  20-30  degrees  positive  required  to  support  the 
weight.  At  this  glide  angle,  the  pitching  moment  would  have  to  be  zero  to 
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PITCH  DAMPING  DERIVATIVES 


MODEL  ANGLE  OF  AT’/*"  ",  DEGREES 
Figure  52.  Effect  of  Angle  of  Attack  on  Pitch  Damping 
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prevent  precession.  As  previously  discussed,  the  frisbce  possesses  a  point 
of  neutral  moment  at  an  angle  sufficient  to  support  the  weight. 

But  if  launched  from  an  airplane  in  level  flight,  this  edge  downward  orientation 
of  the  disk  would  correspond  to  a  large  negative  initial  angle  of  attack.  With 
a  displaced  linear  moment  characteristic  or  the  finite  nonlinear  curves 
observed  so  far,  the  area  under  the  curve  between  launch  and  glide  angle  of 
attack  would  seem  to  provide  an  excessive  precession  transient.  One  brute* 
force  approach  to  this  would  be  to  specify  that  the  pitching  moment  be 
negligible  for  all±  angles  of  attack  encompassing  steady  glide  and  launch. 
Naturally,  it  would  seem  to  be  desirable  to  have  it  negligible  at  all  angles  of 
attack.  Such  an  idealized  moment  curve  with  a  large  neutral  zone  is  sketched 
in  Figure  53. 

The  observed  effects  of  thickness  indicated  in  Figure  45a  trend  in  this  direc¬ 
tion  for  the  cylindrical  flare  as  the  diameter- to- thickness  ratio  decreased 
from  8:1  to  4:1.  Some  published  data  have  indicated  that  the  slope  of  the 
pitching  moment  versus  angle  of  attack  may  be  zero,  at  least  at  zero 
angle  of  attack,  for  ratio  of  1:1,  but  such  a  shape  would  not  stay  near 
zero  long,  and  the  Magnus  rolling  moment  and  side  force  would  undoubtedly 
become  significant. 

It  is  hopeful  that  an  acceptable  combination  of  pitching  and  magnus  charac¬ 
teristics  will  be  obtainable  in  the  3:1  to  2:1  category. 

If  not,  and  thicker  disks  (becoming  cylinders)  are  indicated,  one  might  as 
well  deliberately  design  it  as  a  Magnus  ("flettner")  rotor  with  the  spin  axis 
horizontal  and  flame  out  both  ends.  The  glide  descent  angle  would  be  quite 
steep  for  a  1:1  rotor.  The  glide  velocity  would  be  less  than  150  feet  per 
second  for  a  5  pound  flare  with  the  same  volume  as  the  8  inch  diameter, 

2  inch  thick  prototype  tested  at  Hurricane  Mesa,  decreasing  as  mass  was 
lost  in  a  manner  similar  to  that  described  for  the  disks.  A  1:2  Magnus 
rotor  would  have  a  glide  velocity  less  than  160  feet  per  second  and  a  possibly 
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shallower  glide.  Prespin  would  help  insure  the  desired  orientatioa  while 
unspun  cylinders  with  spinup  ribs  would  have  a  random  direction  and  orienta¬ 
tion  transient.  Prespin  with  the  bottom  edge  retreating  would  tend  to  make 
it  fly  rearward  and  the  opposite,  forward,  depending  on  the  tactical  mission. 

Another  approach  to  generation  of  the  desired  neutral  overturning  moment 
characteristics  was  inspired  by  another  frisbee  variation  toy  which  has 
appeared  since  the  start  of  the  project  and  since  it  was  indicated  that  the 
flare  would  be  a  bottom-burning  cup  instead  of  a  top  and  bottom  burning 
annulus.  This  ring-shaped  toy  appears  to  be  more  stable  when  thrown  in 
the  frisbee  mode  than  the  more  conventional  frisbee-like  devices.  This 
could  be  attributed  to  a  more  neutral  overturning  moment  for  this  "tandem 
airfoil"  arrangement.  The  rear  surface  would  be  generating  a  lift  counter¬ 
acting  the  overturning  moment  of  the  lift  of  the  forward  surface,  whereas 
on  a  disk  or  conventional  frisbee-type  the  net  lift  force  acts  forward  of  the 
center.  Such  an  effect  might  be  possible  if  the  burning  pressures  permit 
airflow  through  the  cavity  to  reach  the  rear  portion  of  the  annulus. 

As  a  matter  of  fact,  the  measured  slope  of  the  frisbee  moment  with  angle  of 
attack  corresponds  to  a  center  of  pressure  location  about  40  percent  of  the 
diameter  aft  of  the  leading  edge  or  20  percent  of  the  radius  forward  of  the 
center  of  gravity,  while  the  standard  figure  for  a  thin  two-dimensional  air¬ 
foil  is  25  percent  aft  of  the  leading  edge  and  was  reflected  in  our  early 
estimates  of  this  disk. 

If  the  rollup  to  glide  that  can  be  obtained  for  a  frisbee  tossed  on  edge, 
tilted  to  the  left  in  a  near  vertical  plane,  and  tentatively  observed  for  light 
sabots  in  another  project  can  be  extrapolated  to  this  dense  body,  then  the 
desired  configuration  would  appear  to  be  a  domed  and  possibly  cupped  body. 

It  would  resemble  a  smooth,  almost  solid  clay  pigeon  or  a  thick,  almost 
solid  frisbee.  This  is  because  the  observed  frisbee  roll  up  is  attributed 
to  the  nose  down  pitching  moment  at  the  small  initial  angle  of  attack,  below 
the  equilibrium  angle,  which  gyroscopically  induces  a  right  hand  tilt  pre¬ 
cession  toward  wings  level. 
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i.  Wind  Tunnel  Model  Configurations  -  Drawings  of  the  wind  tunnel 

models  used  are  shown  in  Figures  54  -  57.  The  frisbee  and  clay  pigeon 

models  were  tested  either  hollow  at  normal  wall- thickness  or  solid  by  the 

use  of  plugs.  The  8:1  cylinder  used  nesting  plugs  to  provide  the  cavity  radii, 

covering  the  range  fpom  solid,  r  /r  =  0,  intermediate  stages,  r  /r  =  0.  25 

and  r  Jr  =  0.  65,  to  hollow,  r  / r  =  0.  95.  In  the  interest  of  simulating  the 
0  0 

averaging  of  the  star  shaped  grain  which  would  come  in  spinning  flight,  the 
cavities  were  circular  and  assumed  to  represent  the  average  annulus 
corresponding  to  the  star  shaped  grain. 

The  4:1  model  tested  on  the  spinning  rig  was  tested  only  in  the  solid  condition. 
The  model  design  support  system  does  lend  itself  to  moderate  cupping  of  the 
bottom  face  and  to  moderate  radiusing  of  the  upper  and  lower  corners.  A 
proposed  series  of  such  tests  to  explore  parametrically  the  effects  of  cor¬ 
ner  radii  and  cavities  within  the  bounds  permitted  by  the  model  is  shown  in 
the  following  subsection. 

A  scheme  for  accurate  measurements  on  a  spinning  hollow  model  has  not 
yet  been  devised.  As  a  matter  of  fact,  the  presence  of  the  side- sting 
mounting  block  could  very  likely  have  introduced  interference  effects  into 
the  nonspinning  hollow  model  data  obtained. 

j.  Recommended  Extension  of  Present  Wind  Tunnel  Test  Data  -  While 
it  is  recognized  that  Magnus  effects  due  to  spin  could  become  more  sig¬ 
nificant  for  the  thicker  shapes  under  consideration,  their  measurement 
would  require  development  of  special  equipment  and/or  data  reduction 
techniques.  It  was  therefore  concluded  that  the  most  expeditious  approach 
to  both  screening  gross  configuration  modifications  for  further  study  and 
isolating  some  of  the  finer  interacting  effects  of  present  configurations 
would  be  to  use  the  present  4:1  spinning  model  in  a  predominantly  non¬ 
spinning  mode  to  explore  these  effects  systematically. 

A  preliminary  plan  for  such  an  exploration  is  presented  in  Table  7,  cate¬ 
gorized  by  purpose  and  priority. 
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Figure  56.  Wind  Tunnel  Model,  8:1  Right  Circular  Cylinder 


I 


Table  7.  Interim  Wind  Tunnel  Test  Plan  for  Follow-On  Studies 


A.  CORNEA  RADIUS  AND  SOTTOM  LIP  VARIATION  (BASIC  4:1  MODEL) 

1.  TOP  CORNER 
>.  PRESENT 

b.  SAME  R  AS  8:1  IF  >  PRESENT 

c.  1/2  INCH  RAOIUS 

2.  TOP  AND  BOTTOM 

•  .  SAME  AS  PRESENT 

b.  SAME  R  AS  8:1  IF  >  PRESENT 

c.  1/2  INCH  RADIUS 

3.  LIP  ON  BOTTOM 

a.  REMOVE  BOTTOM  SECTION  AND  ADD  LIP  SO  THAT  4.1  THICKNESS 
RATIO  IS  MAINTAINED. 

b.  LEAVE  BOTTOM  SECTION  ON  BASIC  MODEL  AND  ADD  LIP  SO  THAT 
THE  ASSEMBLV  HAS  A  2:1  THICKNESS  RATIO. 

B.  THICKNESS  RATIO  VARIATION  (TWO  CORNER  RAIHI  DETERMINED  BV  A) 

1.  6:1 

2.  4:1  (FROM  A) 

3.  2:1 

C.  VARIABLE  CAVITV  MODELS  -  ARBITRARY  SHAPES  AND  SIZES. 


1.  TOP  CORNER  RADIUS  AND  THICKNESS 

a.  TWO  BEST  DETERMINED  BY  A  AND  B. 

b.  LARGER  RADIUS  ON  TOP  (APPROACHING  CLAY  PIGEON  SHAPE). 

2.  CAVITY  DEPTH  AND  RADII  TWO  CAVITY  DEPTHS  (TWO  DIFFERENT  tig) 
WITH  ONE  PLUG  FOR  EACH  BOTTOM. 

TEST  AT  M  =  0.2  = 


0  RANGE 

NUMBER  OF  CONFIGURATIONS 

A. 

•10  TO  40  DEGREES 

7 

B. 

-10  TO  40  DEGREES 

4 

C. 

-10  TO  90  OEGREES 

12 

4. 


Free- Flight  Testing 


Formal  free-flight  testing  as  such  was  not  the  direct  responsibility  of  the 
aerodynamics  contractor  in  the  Self- Suspended  Flare  project  but  the 
contractor  participated  in  planning,  instrumenting,  observing  and  inter¬ 
preting  the  tests  at  Hurricane  Mesa  in  November  1969,  with  respect  to 
flight  behavior  and  aerodynamic  characteristics,  and  was  made  responsible 
for  data  reduction.  This  participation  was  expanded  with  respect  to  original 
plans  in  an  attempt  to  maximize  the  aerodynamic  data  obtained  from  the  test 
series.  The  tests  were  originally  planned  for  inert  flare  trajectory  and  live 
flare  trajectory  and  illumination  tests. 

In  addition  to  these  formal  instrumented  tests  at  Hurricane  Mesa,  limited 
informal  flight  testing  of  clay  pigeons  and  frisbees  was  also  conducted  in 
conjunction  with  the  Hurricane  Mesa  tests. 

Free-flight  testing  for  a  completly  unrelated  project  at  the  site  at  Hurricane 
Mesa  and  in  Minnesota  also  gave  some  interesting  insight  into  the  problem 
because  of  the  observed  apparent  gliding  flight  of  some  of  the  light  weight 
disks  used  for  sabots  in  that  project. 

The  Hurricane  Mesa  tests  tended  to  validate  the  math  models  developed  and 
gave  valuable  insight  into  problem  areas,  but  in  themselves  neither  estab¬ 
lished  feasibility  of  the  concept  nor  yielded  quantitative  aerodynamic  data. 

a.  Hurricane  Mesa  Tests  -  A  log  of  the  configurations  tested  at  Hurricane 
Mesa  in  November  1969  for  the  48  official  launches  is  shown  in  Table  8. 
Almost  5000  feet  of  film  plus  numerous  still  negatives  were  obtained.  This 
table  indicates  the  purpose  of  each  test  number  and  gives  brief  remarks 
on  the  observed  flight  behavior.  The  remark  "OK"  means  that  the  orienta¬ 
tion  appeared  to  be  level  and  the  trajectory  as  expected. 
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Table  8.  Self- Suspended  Flare  Test  Configuration  Log 


TEST 

FLARE 

DIAMETER 

(INCHES) 

WEIGHT 

(POUNDS) 

THICKNESS 

(INCHES) 

LAUNCH 

PRESSURE 

(PSD 

PITCH 

DOWN 

(DECREES) 

PURPOSE 

REMARKS 

1 

-  5 

N/A 

12 

12 

200 

0 

BfcW  CHECKS 

6 

N/A 

12 

12 

? 

200 

0 

STANDARD  INERTS 

FELL  SHORT 

7 

N/A 

12 

12 

2 

200 

0 

STANDARO  INERTS 

FELL  SHORT 

• 

N/A 

12 

12 

2 

200 

0 

STANDARD  INERTS 

OK 

9 

N/A 

12 

12 

2 

200 

0 

STANDARD  INERTS 

OK 

10 

N/A 

5 

2 

200 

0 

STANDARD  INERTS 

OK 

11 

N/A 

5 

2 

200 

0 

STANDARD  INERTS 

OK 

12 

N/A 

5 

2 

200 

0 

STANDARD  INERTS 

OK 

13 

N/A 

5 

2 

200 

0 

STANDARO  INERTS 

OK 

14 

H*l 

12 

3.5 

1.312 

160 

0 

SOCIAL  -  TWO  LAYERS  WOOD 

PEELED  OFF  TO  LEFT. 

IS 

M-4 

12 

7.25 

1.39 

180 

0 

SPECIAL  -  TWO  LAYERS  WOOD 

SLOWER  PEEL  OFF-STOPPED 

ROLL  AND  CLIOED. 

16 

M-3 

12 

3.5 

1.2 

160 

16.5 

SPECIAL  -  TILT  DOWN 

ROLLED  ON  EOCE. 

17 

M*2 

12 

3.5 

1.187 

160 

16.5 

SPECIAL  -  REPEAT  OF  TEST  16 

ROLLED  ON  EDGE. 

18 

STANDARD 

12 

12 

2 

200 

17.5 

COMPARE  STANDARD  WITH 

TESTS  16  AND  17 

STRAIGHT 

19 

LIVE 

5 

2 

200 

17.5 

LIVE 

LEVELED  OFF. 

20 

LIVE 

5 

2 

200 

0 

LIVE 

WENT  STRAIGHT  UP. 

21 

H*  17 

2.25 

1.75 

200 

0 

LIGHT  8HNCH  SPECIAL 

ROLLED  OFF  NEAR  END. 

22 

H-18 

::  n 

3.7 

1-15/16 

200 

0 

MID  8-INCH  SPECIAL 

CLOSE  TO  8-INCH  STANDARD. 

23 

LIVE 

5 

2 

200 

0 

LIVE 

DOWN  AND  STRAIGHT  UP. 

24 

LIVE 

5 

2 

200 

0 

LIVE 

DUD. 

25 

-  33 

LIVE 

5 

2 

200 

0 

LIVE 

DUD  OR  CLIMBED. 

33 

•  36 

N/A 

12 

12 

2 

200 

0 

REPEATS  OF  12-INCH 

STANDARD  INERTS 

TEST  35  IS  SHORT  OF  OTHER. 

37 

H-10 

12 

7 

11-1/4© 

4-1/2© 

1-14/32 

200 

0 

FRISBEE  SPECIAL 

PELLED  OFF. 

38 

H-15 

12 

2 

200 

0 

FRISBEE  SPECIAL 

SOMC  ROLL.  | 

39 

H-19 

12 

2 

180 

0 

FRISBEE  SPECIAL 

ROLLED  ON  EDGE. 

40 

M-ll 

12 

9-5/8 

2 

160 

0 

FRISBEE  SPECIAL 

CONTINUOUS  PRECESS  AND  CLIDE. 

41 

H-21 

8 

4-1/2 

2 

180 

0 

FRISBEE  SPECIAL 

CONTINUOUS  PRECESS  - 

SOMC  PAUSE .  ; 

42 

H-7 

• 

3-1/2 

1.9 

160 

0 

FRISBEE  SPECIAL 

CONTINUOUS  PRECESS. 

43 

H-6 

8 

2 

1-8 

160 

0 

FRISBEE  SPECIAL 

CONTINUOUS  PRECES5. 

44 

H-22 

8 

5 

2 

160 

16 

STANDARD  FLARE  PITCHED 

DOWN 

NIL  PRECESS 

45 

H-24 

8 

2-1/4 

2 

140 

16 

LICHT  CYLINDER  SPECIAL 

SLOW  PRECESS  AND  GLIDE. 

46 

H-9 

8 

2 

1-3/4 

140 

16 

FRISBEE  SPECIAL 

PRECESS  ANO  CONSIDERABLE 
TRAVEL  TO  LEFT. 

47 

M-23 

8 

2-1/4 

140 

0 

LICHT  CYLINDER  SPC(  IAL 

SOME  PRECESS 

48 

N/A 

8 

5 

2 

200 

0-16  DECREE 
ROLL 

EFFECT  OF  ROLL  ON  STANDARD 
FURE 

NIL  PRECESS.  100  FEET 

TRAVEL  TO  RIGHT. 

Q  HONIVWCIL  A*D  Dm  «EC0I0S  ADC  DCVCAJCD  ON  THC  MIGHT  Or  TISTS  3«  AND  39. 
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b.  Hurricane  Mesa  Test  Instrumentation  -  The  photographic  instrumenta¬ 
tion  arrangement  at  Hurricane  Mesa  is  illustrated  in  Figure  58.  The  equip¬ 
ment  was  divided  between  the  launch  site,  the  nearby  control  site,  and  the 
tracking  site  which  was  placed  on  another  point  of  the  mesa  which  gave  full 
view  of  the  launch  and  impact  area.  The  tracking  site  was  730  feet  from 
the  launcher  and  the  launcher  was  aimed  perpendicular  to  the  line  joining 
the  launcher  and  tracking  camera.  The  tracking  camera  and  launcher  were 
within  inches  of  a  level  plane.  The  separate  sites  were  illustrated  in 
Figures  11  through  15. 

In  addition  to  the  launcher  at  the  launch  site,  (also  called  boresight)  there 
was  the  boresight  LoCam  tracking  camera  with  a  12-240  mm  power  zoom 
lens  above  the  launcher  on  a  forklift,  a  test  number  board  visible  to  the 
tracking  site  cameras,  and  a  launch  signal  fla$h  bulb,  plus  camera  synchro¬ 
nization  and  timing  electronics.  The  power  zoom  motor  was  preprogrammed 
and  was  synchronized  with  the  launch  signal. 

Launch  velocity  was  controlled  by  the  pressurization  level  which  had  been 
previously  calibrated  for  the  standard  flares  (50  feet  per  second  and  1000 
rpm  for  a  12  inch  diameter,  12  pound  inert  flare  and  75  feet  per  second 
and  1500  rpm  for  an  8  inch,  5  pound  inert  flare  at  200  psi).  The  launch 
condition  had  been  found  to  be  quite  repeatable.  Calibration  data  were  not 
available  for  the  field  modified  special  units,  but  it  was  hoped  to  obtain 
velocity  from  the  tracking  data.  The  launcher  was  released  by  a  solenoid 
on  electrical  signal  from  the  launch  control  box,  or  could  be  fired  manually. 

Below  the  launch  site  on  the  sloping  shoulder  of  the  mesa  was  placed  a  grid 
of  10  5-foot  long  poles  with  alternating  1-foot  red  and  white  stripes  laid 
horizontally  at  1  degree  intervals  of  line  of  sight  depression  from  the  bore¬ 
sight  station.  These  are  barely  discernible  in  the  still  photos,  but  were 
visible  to  the  color  tracking  cameras  when  they  were  in  the  field  of  view. 

This  grid  was  expected  to  provide  adequate  definition  of  flare  position  rela¬ 
tive  to  the  boresight  camera,  and  the  rugged  terrain  discouraged  more 
extensive  arrangements.  All  earlier  flights  had  shown  virtually  planar 
trajectories  (see  Figure  59). 
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The  impact  grid  survey  conducted  from  the  boresight  forklift  was  supplented 
by  Hasselblad  still  photos  of  the  impact  area  and  surrounding  terrain 
traversed  by  the  flare. 

At  the  control  site,  a  safe  distance  from  boresite  on  the  edge  of  the  cliff, 
was  located  the  launch  control  box  which  electronically  controlled  launch  of 
all  flares  and  ignition  of  the  live  flares  with  programmed  delays  and  provided 
event  pulses  to  the  camera  synchronization  circuit  and  the  flashbulb.  The 
launch  countdown  was  transmitted  by  radio. 

At  the  tracking  site  there  was  a  LoCam  tracking  camera  with  a  300  mm 
telephoto  lens  mounted  on  a  panhead  with  synchro  pickoffs  for  azimuth  and 
elevation,  a  recording  box  which  contained  the  azimuth  and  elevation  read¬ 
outs  and  the  recording  camera,  a  fixed  aim  Milliken  panoramic  motion 
picture  camera,  and  associated  timing  electronics.  These  were  supple¬ 
mented  by  a  Hasselblad  panoramic  still  camera  for  the  live  flare  tests 
and  the  terrain  mapping.  The  focal  length  of  the  Milliken  camera  was 
changed  from  17  mm  for  the  day  shots  to  10  mm  for  the  night  shots. 

Figure  58  indicates  the  relationship  of  the  respective  camera  lines  of  sight 
and  of  the  fixed  field  of  view  of  the  Milliken  camera.  On  this  figure, 
nominal  trajectory  refers  to  the  projection  of  the  apparent  trajectory  as 
seen  from  the  tracking  site  onto  the  nominal  ballistic  vertical  plane  through 
the  launcher  aim  axis.  If  there  were  no  lateral  deviation,  the  actual  and 
nominal  trajectories  would  coincide.  Because  the  cameras  were  aimed  down¬ 
ward  and  consequently  inclined  to  this  nominal  plane,  there  was  some  slight 
distortion  of  the  trajectory  images  that  was  corrected  for  in  the  data  reduc¬ 
tion. 

In  addition  to  transit  surveys  of  the  site  locations  and  the  trajectory  grid, 
calibration  procedures  included  checking  the  operation  of  the  boresight 
power  zoom  lens  against  dummy  flares  at  measured  distances.  This  was 
found  to  be  sufficiently  repeatable  such  that  one  standard  zoom  lens 
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calibration  curve  against  time  was  used.  This  lens  calibration  curve  of 
focal  length  versus  time  which  was  used  in  conjunction  with  boresight 
image  size  versus  time  (Figure  60)  to  determine  slant  range  from  bore- 
sight  to  the  flare,  which,  in  turn,  was  required  to  convert  image  apparent 
position  to  cross  range  position. 

The  flashbulb  at  the  launcher  was  intended  to  provide  a  launch  signal  in  the 
field  of  view  of  the  cameras  to  backup  the  synch  signals  on  the  film  margin 
transmitted  by  cable.  It  was  not  always  in  the  field  of  view  of  the  zoom  lens, 
however.  The  repeatability  of  the  framing  rate  of  the  tracking  cameras, 
which  were  allowed  to  build  up  to  speed  during  the  countdown,  not  only  pro¬ 
vided  a  backup  to  the  timing  light  signals  on  the  film  margin  but  greatly 
simplified  the  subsequent  data  reduction.  The  semi-automatic  film  reading 
equipment  to  be  used  could  both  count  frames  traversed  and  advance  a 
predetermined  frame  interval. 

c.  Typical  Hurricane  Mesa  Results  for  Standard  Inert  Flares  —  The 
planar  trajectories  obtained  agreed  fairly  well  with  predictions  for  the 
launch  conditions,  and  the  observed  precessions  tended  to  validate  the 
precessional  computer  programs.  There  was  reasonable  correlation  with 
the  analytical  criteria  solutions,  but  the  unrealistically  rapid  change  in 
angle  of  attack  plus  the  velocity  variations  makes  it  difficult  to  correlate 
such  linearized  analyses  closely  and  makes  it  impossible  to  use  these 
linear  solutions  to  extract  the  desired  aerodynamic  coefficients.  A  better 
flight  simulation  would  have  required  much  higher  launch  velocities  from  a 
sled  or  aircraft  to  keep  the  angles  small. 

While  not  successful  in  producing  quantitative  data,  not  too  much  was 
expected,  and  these  tests  were  quite  valuable  because  they  did  provide 
qualitative  validation  of  the  math  models  and  insight  into  the  dynamic 
criteria,  revealed  problems  in  flare  design,  gave  experience  in  field  test¬ 
ing  and  data  collection,  and  led  to  experience  in  data  reduction. 
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Figure  60.  Typical  Boresight  Camera  Flare  Image  Size 
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A  typical  trajectory  extracted  from  the  flight  test  photographic  data  is 
plotted  in  Figure  61  for  a  standard  12  inch  diameter,  12  pound  dummy  flare 
(Shot  36).  This  figure  indicates  the  errors  present  in  the  cross  range  data 
extraction  and  reduction,  but  shows  a  generally  planar  motion  as  expected 
and  observed  visually.  The  results  agreed  fairly  well  with  computer  pre¬ 
dictions  both  before  and  after  the  tests  (see  Figure  62). 

The  dotted  nominal  line  on  Figure  61  represents  the  corrected  projection 
of  the  trajectory  observed  from  the  tracking  site  onto  the  nominal  vertical 
ballistic  plane.  Since  there  was  little  cross-range  motion  there  is  little 
difference  between  it  and  the  fitted  actual  computed  with  appropriate  curve 
fitting.  For  that  matter,  there  is  little  difference  from  the  uncorrected 
raw  data  obtained  from  terrain  matching  in  Figure  63.  The  primary 
difference  is  the  curve  fit  smoothing  of  the  undulations  in  the  raw  data. 

Since  these  undulations  were  observed  in  other  cases,  they  may  very  likely 
have  been  real  rather  than  data  errors.  If  real,  they  could  be  attributed 
either  to  an  energy  interchange  oscillation  similar  to  the  "phugoid"  mode 
of  aircraft  or  merely  to  the  gustiness  of  the  winds. 

For  Shot  36,  good  correlation  was  also  obtained  between  terrain  matching 
of  Milliken  and/or  tracking  camera  observations  from  the  tracking  site  and 
positions  in  the  nominal  plane  computed  from  the  azimuth  and  elevation 
data  recorded  by  the  data  camera  from  the  tracking  camera  synchro 
pickoffs  at  the  tracking  site.  These  redundant  data  were  not  always 
available  and  still  required  cross-range  matching  and  correction  for  non- 
planar  trajectories.  It  did  serve  to  help  validate  the  method  eventually 
adopted  for  all  shots  that  were  reduced, however. 

d.  Typical  Hurricane  Mesa  Results  for  Special  Inert  Flares  -  A  typical 
extremely  nonplanar  trajectory  extracted  (and  partially  extrapolated)  from 
the  flight  test  photographic  data  is  plotted  in  Figure  64  for  Shot  45,  which 
was  a  special  light  weight  (2  1/4  pounds),  8  inch  diameter  cylinder  launched 
16  degrees  below  the  horizontal.  The  precessional  motion  demonstrated 
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Figure  61.  Typical  Flight  Test  Trajectories  for  a 
Standard  Inert  Flare 
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32.  Comparison  of  Typical  Computer  Trajectory  Simulation 
with  Flight  Test  Data 
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Figure  64.  Typical  Flight  Test  Trajectory  for  a  Special  Light 
Weight,  Inert  Flare 


quickly  took  the  boresight  field  of  view  out  of  view  of  the  range  grid,  and 
actually  took  the  flare  out  of  view  of  the  tracking  camera  before  impact. 

The  terrain  matching  profile  from  the  tracking  site  was  similar  to  that  for 
the  standard  flare.  The  fitted  actual  trajectories  required  some  extrapolation 
of  the  nominal  plane  trajectory  to  match  the  time  to  impact  of  the  cross 
range  trajectory.  Thus,  the  solid  line  actual  in  the  figure  departs  from 
the  dotted  line  projection  on  the  nominal  plane.  The  plotted  cross  range 
trajectory  corresponds  roughly  to  the  qualitative  observation  of  Table  8 
for  Shot  45-that  the  flare  first  precessed  considerably  to  the  left  then 
appeared  to  adopt  a  relatively  steady  glide  in  that  new  direction. 

Attempts  to  extract  meaningful  values  for  the  other  trajectory  parameters 
for  this  case  or  to  find  a  corresponding  computer  simulation  were  fruitless 
within  the  time  and  data  constraints. 

e.  Typical  Flight  Test  Results  for  Live  Flares  -  The  trajectory  of  a 
typical  thrusting  live  flare  is  shown  in  Figure  65.  This  corresponds  to  the 
still  picture  Figure  17.  This  trajectory  was  calculated  from  the  data 
camera  readouts  of  the  tracking  camera  line  of  sight  from  the  tracking  site, 
assuming  that  the  flare  stayed  in  the  nominal  vertical  plane,  since 
obviously  there  were  no  boresight  data  available.  The  trajectory  of  Figure 
65  shows  a  rise  higher  than  the  photo  of  Figure  17  because  the  still  camera 
field  of  view  cutoff  at  about  350  feet  above  the  mesa,  while  the  Milliken 
motion  pictures  cutoff  just  slightly  over  1 50  feet  above  the  Mesa.  The 
fields  of  view  had  been  selected  to  cover  the  intended  downward  trajectories. 
Unfortunately  for  some  reason,  no  quantative  data  were  retrievable  on  the 
tests  of  Figure  16,  the  only  one  for  which  the  illumination  data  would  have 
been  meaningful. 

Figure  66  illustrates  the  available  data  on  similar  live  shots  (27  and  29) 
obtained  from  the  Milliken  camera  data,  and  compared  to  the  corresponding 
portion  of  Shot  32. 
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Figure  65.  Typical  Live  Self- Suspended  Flare  Trajectory,  Shot  32 
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Figure  66.  Typical  Live  Self- Suspended  Flare  Trajectories, 
Shots  27,  29,  and  32 
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The  data  extraction  program  being  used  was  even  less  successful  in  extract¬ 
ing  valid  trajectory  parameters  from  these  data  because  of  the  sharp  corner, 
the  limited  number  of  points,  and  the  unknown  weight  and  thrust.  The  shapes 
do  qualitatively  match  the  results  of  constant  thrust  simulations  performed 
independently  at  the  University  of  Denver  and  on  the  trajectory  simulation 
described  in  this  report  in  section  IV.  B.  1 . 

f.  Flight  Test  Results  on  Spin  Decay  -  The  spin  rate  of  the  flare  was 
monitored  in  the  inert  flare  tests  by  the  boresite  camera,  and  extracted  by 
counting  revolutions.  For  the  standard  case  (Shot  36),  the  spin  rate  was 
very  close  to  the  calibrated  value  of  1000  rpm  with  scatter  between  923 
and  1000  rpm.  During  the  available  flight  time  at  low  velocity  and  high 
angle  of  attack,  there  was  no  noticable  decay  upon  which  to  base  an  estimate 
of  spin  damping  coefficient  to  use  in  predicting  decay  at  full-scale  velocities. 
This  lack  of  decay  under  these  conditions  does  provide  an  indication  that 
actual  spin  decay  will  probably  not  be  sufficient  for  a  reliable  SAF  signature. 

While  the  swirling  plumes  could  very  likely  have  provided  spin  rate  data 
on  the  live  flares,  which  appeared  to  stay  high,  the  generation  of  such 
high  thrust  and  the  observed  burned  through  spots  most  likely  also  con¬ 
tribute  spin  torques  which  would  have  distorted  the  interpretation.  Mass 
flow  from  an  effective  thrust  nozzle  would  contribute  to  damping  if  dis¬ 
placed  from  the  centerline. 

g.  Correlation  of  Flight  Test  and  Predicted  Results  -  Attempts  to 
extract  aerodynamic  coefficients  from  the  photographic  data  were  singularly 
unsuccessful  since  it  was  difficult  to  obtain  even  position  data  for  those 
flares  which  departed  significantly  from  the  nominal  ballistic  plane.  For 
the  planar  case  example  (Shot  36)  with  a  12  pound,  12  inch  dummy  flare, 
some  approximate  extraction  was  performed  and  compared  to  the  most 
representative  computer  simulation  of  the  same  flight  condition. 
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In  Figure  27  the  position  correlation  appears  to  be  quite  good  at  the  6.  4 
second  impact  point.  It  is  somewhat  disconcerting  to  note,  however,  that 
the  computed  trajectory  departs  considerably  from  the  planar  for  flight 
times  much  greater  than  that  attainable  at  Hurricane  Mesa. 

As  the  number  of  numerical  differentiations  required  to  obtain  the  parameters 
increased,  the  errors  increased.  Position  resolution  errors  against  the  sky 
background  or  launcher  pressurization  errors  could  have  contributed  to 
discrepancies  in  velocity  near  launch. 

h.  Auxiliary  Observations  from  Other  Flight  Tests  -  Some  preliminary 
indications  of  a  possible  approach  solution  to  the  flare  configuration/delivery 
problem  was  revealed  by  other  tests  at  Hurricane  Mesa  (horizontal  barrel) 
and  in  Minnesota  (elevated  barrel)  for  an  unrelated  project.  The  tracking 
camera  was  beside  the  gun.  Sabots,  each  consisting  of  two  6  inch  diameter, 

2  inch  thick  disks  (3:1  cylinders)  of  lightweight  foam,  were  used  in  rectan¬ 
gular  barreled  gas  guns  to  impart  a  controlled  forward  velocity  and  angular 
rate  to  a  munition  held  between  them.  The  axis  of  the  angular  rate  was 
horizontal  and  perpendicular  to  the  disks. 

In  some  tracking  photos  from  the  Hurricane  Mesa,  where  maximum  velocities 
were  600-800  feet  per  second  with  level  barrels,  there  were  some  observa¬ 
tions  that  the  sabots  appeared  to  roll  out  and  glide  with  the  forward  edge 
downhill  in  a  relatively  stable  orientation.  For  the  Minnesota  tests  at 
velocities  up  to  1300  feet  per  second  with  elevated  barrels,  there  were 
visual  impressions  that  the  sabots  rolled  out  into  a  stable  orientation, 
stopped,  and  then  glided  back  toward  the  gun.  The  tracking  cameraman 
was  the  same  individual  who  operated  the  tracking  camera  for  the  November 
1969  Self- Suspended  Flare  tests  at  Hurricane  Mesa. 

It  has  also  been  noted  that  toy  frisbees  tossed  on  edge  will,  in  some  cases, 
roll  out  to  the  desired  glide  orientation.  It  is  speculated  that  this 
phenomenon  could  be  exploited  to  eliminate  the  necessity  for  complex 
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tilting  of  the  launcher  and  or  brute- force  neutralization  of  the  aerodynamic 
moments,  by  letting  the  flare  seek  its  own  preferred  glide  orientation. 

5.  Math  Model/ Computer  Program 

a.  Basic  -  The  equations  used  to  develop  the  trajectory  computer  pro¬ 
gram  are  presented  in  Appendix  B.  These  equations  are  based  on  the 
aerodynamic  force  and  moment  model  of  Figure  5.  They  used  an  inter¬ 
mediate  set  of  nonspinning  aeroballistic  axes  for  the  flare  to  perform  the 
basic  calculations  with  appropriate  transformations  between  these  flare 
axes  and  the  aerodynamic  axes.  The  flare  axes  were  arbitrarily  chosen 
to  simplify  the  inertial/gyroscopic  and  kinematic  calculations  and  to 
minimize  the  possiblity  of  singularities  in  the  differential  equations. 

While  a  standard  conventional  sequence  of  Euler  angles  in  the  yaw,  pitch, 
and  roll  sequence  was  used  to  describe  the  instantaneous  orientation  of  any 
arbitrary  set  of  axes  in  the  flare,  the  axis  transformation  matrices  were 
computed  using  the  Euler  symmetrical  parameter  (sometimes  called 
quaternion)  approach  which  uses  four  parameters  instead  of  three,  but 
avoids  the  singularities  and  is  actually  simpler  to  program  than  other 
methods. 

The  fundamental  assumption  made  in  these  equations,  in  order  to  suppress 
the  high  frequency  nutational  oscillations  but  allow  quasisteady  rotational 
precession  and  full  translational  freedom,  was  that  the  derivatives  of  the 
pitch  and  roll  angular  rates  were  negligible.  Inverting  the  resulting 
equations  yielded  the  equations  for  gyroscopic  precession  under  the  applied 
aerodynamic  moments  resulting  from  angle  of  attack,  spin,  configuration 
shape,  and  angular  rates. 

Once  the  orientation  of  the  arbitrary  flare  axes  were  all  obtained  in  terms 
of  its  transformation  matrix,  the  orientations  of  several  other  arbitrary 
axes  were  readily  obtainable  and  were  output  by  the  program  by  making  the 
appropriate  auxiliary  axis  transformation  required  and  extracting  the 
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orientation  Euler  angles  by  manipulation  of  the  matrix  elements.  No  one 
set  of  axes  or  angles  seems  to  best  describe  the  orientation  of  such  a  disk 
type  body,  which,  if  and  '"hen  gliding  could  be  described  by  either  aero¬ 
dynamic  or  velocity  axes  but  unlike  an  airplane  the  wing  tips  are  not  visible. 
Consequently,  several  choices  of  output  axes  are  availaMe. 

b.  Linearized  Analysis  -  The  equations  used  for  the  stability /precession 
criteria  analytical  solutions  were  essentially  those  of  the  computer  program 
but  linearized  with  respect  to  the  velocity  axis  description  of  the  flare 
orientation  and  forces  (see  the  equations  of  Appendix  C).  In  this  axis  system, 
the  roll  orientation  of  the  lift  vector  determines  the  rate  of  turn  of  the 
trajectory  in  gliding  flight.  This  system  eliminated  the  effect  of  linear 
acceleration  and  lent  itself  to  varying  levels  of  degrees  of  freedom  for 
study  of  the  precessioa  with  and  without  translation  and  with  and  without 
wobble. 

The  linearization  was  performed  with  respect  to  an  ideal  planar  trajectory, 
which  would  assume  no  precessioa  and  the  equations  presented  were  solved 
for  the  simplest  case  of  level  launch,  as  presented  in  section  IV.  B.  2.  The 
solution  should  be  extended  for  the  more  general  and  more  proper  steady- 
state  glide  angle  conditioa  if  identifiable,  for  arbitrary  launch  conditions. 

c.  Typical  Aerodynamic  Coefficients  and  Inertial  Data  for  the  Computer 
Simulation  -  Typical  values  of  the  geometric  parameters  used  to  define 
the  flares  in  the  computer  runs  are  tabulated  in  Appendix  D. 

Static  aerodynamic  coefficients  CN>  CA,  and  CM  used  for  computation  with 
the  trajectory  program,  are  tabulated  in  Appendix  D  as  functions  of  the 
angle  of  attack  for  several  typical  configurations.  All  runs  used  the  same 
set  of  damping  coefficients  and  also  assumed  that  the  Magnus  moment  was 
negligible  (KM  =  0)  since  the  results  were  insensitive  to  these  parameters 
at  the  anticipated  levels.  The  configuration  numbers  as  tabulated  in  Appendix 
D  were  arbitrarily  assigned.  For  the  8:1  cylinder  with  variable  cavity  the 
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static  coefficients  are  tabulated  as  functions  of  angle  of  attack  for  four  values 
of  the  cavity  radius  ratio  to  outside  radius  (identified  as  Configuration  8C). 

d.  Sample  Simulation  Results  -  Samples  of  computer  inputs  and  outputs 
are  shown  in  Tables  9  and  10  for  a  case  with  cavity  dependent  aerodynamics 
input  as  four  sets  of  curves  and  with  constant  thrust,  and  in  Tables  11  and 
12  for  an  inert  case  with  only  one  set  of  aerodynamic  curves. 

e.  Possible  Math  Model/Computer  Program  Modifications  -  The  math 
model  presently  assumes  that  the  side  force  due  to  spin  (Magnus  force) 

is  negligible  and  that  the  pyrotechnic  burning  rate  and  thrust  are  constant. 

As  the  disk  thickness  increases,  this  Magnus  force  could  become  significant 
and  actually  provide  a  lifting  force  for  units  whose  spin  axis  is  nearer 
horizontal  than  vertical.  This  term  could  easily  be  added  in  the  same 
form  as  the  Magnus  moment  in  the  math  model. 

The  equations  for  mass  loss  and  thrust  would  merely  have  to  be  changed  to 
integrals  of  the  desired  arbitrary  burning  rate  and  the  equation  relating 
thrust  to  this  burning  rate  supplied. 

Associated  computer  program  modifications  would  be  to  introduce  the 
Magnus  moment,  the  Magnus  force  coefficients  (if  added),  and  the  spin 
damping  (when  available)  as  functions  of  angle  of  attack  and,  if  possible, 
cavity  radius,  and  to  introduce  the  burning  rate  functional  inputs  and 
add  them  to  the  differential  equation  solution. 

To  provide  finer  resolution  of  some  of  the  more  nonlinear  curves,  primarily 
for  cases  where  large  thrust  values  may  drive  the  angle  of  attack  to  -90 
degrees  in  vertical  flight,  the  data  input  should  be  modified  to  provide 
either  a  larger  number  of  data  points  per  variable  than  the  present  10  or, 
preferably,  a  selectable  number  for  each  coefficient  which  can  be  tailored 
to  the  available  data. 
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Table  9.  Original  Input  for  a  Series  of  Simulation  Runs 

with  Variable  Cavity  Aerodynamics,  Thrust,  etc. 

RUN  11-10  CONFIG  8C  -  INERT  -  NO  PRECESSION  -  NO  THRUST 


V  200, 

GA  0.0  -  -  .... 

HV  0.0 

THETA  .  0.0 

PHI  0.0  . . 

OMEGA  1000, 

H  .  .  5000, . .  .  . —  .... 

VA  0.0 

HP  2! 

HP  - 5.  _  ...  _  .  _  _  ... 

RHOP  2.86 

Oh i  0.0 

Tl  100. 

TF  15. 

HF  -1. 

KP  0.  .  -  . -  ... 

KM  0. 

THRUST  -  .  0... -  — - 

TINT!  .2 

TINT*  1.0 

TSTOP  1. 

CM0  -  . 

•10.  0.  10.  20.  30.  40.  50.  A0.  70,  90. 

-15,  -.15 . -.15  .  -.15  ... .  -.15 _ -.15 _ -.15 _ -.15.  .  -.15  -.15 

CLP 

•10,  0.  10,  20.  30.  43,  50.  6Q.  70.  90. 

-15,  -.15  -.15  -.15  -.15  -.15  -.15  -.15  -.15  -.15 

NCURVE  4. 

0.0 

•10.  -5.  .  .0.0 _  5.0  10. _ 15,  _  .  ..  30 ,  ...40,.  60.  90. 

-.46  -.26  -.04  ,17  .38  ,52  .86  1.05  1.74  1.2 

-.03  -.045  -.02  .035  .045  ,033.  .066  .076  .0(2  0.0 

,135  .135  .136  .1365  .138  ,165  .155  .111  .057  0.U 

,25  •• 

-10,  -5.  0.0  5.  10.  15.  30.  40.  60,  90. 

-.46  _ -.25  _ _ -.02— .  .2  _ .42 _ ,‘.»4 _ .87 _ 1.08  .  1.3  1.21 

-.03  -.045  -.02  .035  .045  .(133  .068  .078  .062  0.0 
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Table  9.  Original  Input  for  a  Series  of  Simulation  Runs 

with  Variable  Cavity  Aerodynamics,  Thrust,  etc. 
(Concluded) 


.125 

.63 

.1265 

,129 

,132 

_ .135 _  ,16 

.165 

.102 

.049 

0.0 

-10, 

•5  • 

0,0  . 

5. 

.  10.  16. 

30. 

45, 

70. 

90. 

•  .46 

-.25 

-.02 

,2 

.42  ,56 

.85 

1.13 

1.22 

1.19 

-.03  ... 

-.02. 

_  -.003.. 

.  ,023 

_ .042 _ ,04 _ 

_  .074 _ 

.09  _ 

,064 

0.0 

.15 

1.0 

.15 

.149 

,15 

.151  ,175 

.165 

.12 

,07 

0.0 

•10. 

•5. 

0.0 

5. 

12.  19. 

30. 

45. 

60. 

90. 

-.46 

-.2 

.08 

,36 

,78  1,18 

1.37 

1.46 

1.43 

1.26 

-.024 

.005 

.022 

.038 

0.0  -.01 

.007 

,024 

,022 

0.0 

.135 

,145  . 

-..155  — 

,17 

_ .202 _ ,24 - 

— ,17. 

,105.. 

.  .059 

0.0 

RUN  11- 
KP 

20  CONFIG  8C  - 
1. 

INERT 

•  PRECESSING  -  NO 

THRUST 

• 

HUN  11-30  CONFIG  8C  •  SUNNING  •  VARIABLE  AERO  -  NO  PRECESSION  •  .5  LB  THRUST 

_  THRUST  .» _ _ _ _ _  .  .  . 

TI  0. 

■  KP  0, .  . 

-  RUN  11-40  CONFIG  8C  -  BURNING  -  VARIABLE  AERO  -  -PRECESSING  -  ,5  LB  THRUST 
KP  1.0 


r 

■ 


Table  10.  Sample  Output  for  a  Typical  Computer  Run  with 
Cavity  Aerodynamics,  Precessing,  with  Fixed 
Format 


AL**fA 

rtL/H  «  u.U'JU 
•A.UUwMUA 

-3,  woi.fc^uw 

U. 

3.UlJ0h*»U 

l.UlHk.Ul 

1.9l'Uh»Ul 

0 ,  UUUb*0 A 

4,wnuF*0X 

6,g0gE»ri 

g,gngt«01 

Lia  — 

_  _-4*6UUc*0i 

-^*ULbfc-Ul 

•4.kUlib“tU 

1,/bilfc-ul 

UiOUgt-Vl 

O«6UUh«01 

4,U*UF«0ti 

;,44gE»r,o 

i.aoufc.og 

La 

•4 . t  uue-iu 

4v»UUtr«Urf 

4.9UUt-U, 

J..»l'Ufc"U4 

6,6l'Ub«U^ 

/,»l<UE*n< 

*, 449f.fi 

Ut 

LA 

A,  *>SUt*Ui 

1, obbK-Ui 

1. 40UI--U1 

A.  J03b«M 

i.0«Uk*Ul 

I.eUUfc.Ul 

1 .99Ub-UA 

X.ilUE-01 

>,/(i9E-r* 

u. 

AlHma 

HL/N  a  .JbU 
•1. J juI  ♦ll 

O.yl'lfUu 

J. 

9.UUUt*''U 

l.b«Uk*bl 

1.9UUtM'l 

4.  uUUbeMA 

4i,0UE*01 

B.gngf.r.i 

»,g0ut.0l 

t* 

■ u i 

•  < ,  9i'b»-«U* 

-X.ClIUflU 

.UUUK-lH 

... Uul.OX 

9.41'Uh-Ul 

o • 7UUb«UA 

i,g«OE*og 

i,440E.r« 

X.aiub.og 

La  - 

- 

•*1.0LkL*UL 

"<altVbt~V< 

*i  3UUfc-0* 

«.9Ubt*«* 

..JUUh-u* 

«/,buyb-U4 

/,anvE*o< 

*,4i)gE.r.a 

Ut 

l« 

i«<9ul  ••11 

i.*63k*Ui 

ii<HkallA 

i  .Ak'UK-l  l 

i.i9Uk*Ul 

1.0liUfc-Ul 

1  •  6*5  Uk  •O  A 

Xil)?UE*0l 

4.V99ET2 

Ut 

ALHh* 

Kc/it  i  ,odu 
•A.k'Uk^Ul 

•  5 ,  i»Ml»t-*UU 

0. 

>, U0uh*b J 

1. wU.k.Ul 

l,til>Uh«Ul 

4. UMUbeOA 

<I»PUE*01 

/.gngf.fi 

y.gftufc.oi 

CL 

•  4,CUuL*«Ul 

-4 .yutL-ui 

-A,  kUlit-lU 

* .UlUb*tA 

4.,Uwk>Ul 

9.#lgfc-Ul 

b  •  9Uub«UA 

X,13UF*0^ 

l,4>Uf»tU 

ltlVUb.OU 

L* 

•  )  i  b  (-  k  F  •  t  4 

•fc  ,  U<tlh- U* 

•J.U'iUb*U4 

«.<ligh.g. 

4,uruh-o«: 

/•4Mvb«U4 

,,U0UF-0, 

6t44gfc.f  i 

Ut 

LA 

-  -AH^i»4L»ua. 

x.  JLuL-Ui. 

A.lVUfc-UA 

1.3Utb*lA 

l.»luF*01 

l./9gh*Ui 

A.b^ub-UA 

X,<0ufc-f* 

1 1 949F»f a 

u. 

SU*M 

Hw/rt  «  1.UIU 
•laUi'Jb  +  liA 

•9. lUUL^UU 

J. 

UiUOOt'  +  l  u 

l.u"gh»ui 

A  •  tf  1‘Ub^OA 

4,9PgF*0i 

»,g4UE*n 

g.uoufc.oi 

La 

•4,o.  w**»l  * 

-r  .UUul-01 

0. gMUfc- U< 

4.6l)u^*l  1 

/.oggh-Hl 

l.»hgh»ug 

A  **>/Ub*Uii 

*,4f gfc.g j 

l,44UF»f u 

X.a6gt*og 

La 

3<(JlUt«lld 

4.  4 

Ut 

-It,'  ,►•1)4 

/  .yuub«0A 

4(4ngE»n, 

4i4ngF*ra 

g. 

L  A 

1  •  j*g»«Ux 

.  .HSgfUl 

1.33Ub-»)A 

L%  A 

,.u4uh.0i 

l. <tl‘Uh. til 

1 . 7<IUb«0A 

i.usup-g* 

»,»4gE.',a 

g. 

AU'aA 

u« 

A.uUUfc*UA 

4  «  U l It  ♦tl 

9it.guh.Ul 

4.g"gh«n 

•tVObb.Pl 

/,li!lgF»i  1 

v,unub*ni 

C.TJ 

-1 .31  u*  »ni 

-i.yiibt—  Ul 

•4.3U«b*UA 

•  1  •  &0U4--I  1 

-i.9l)utr*Ul 

*1 ,9LUt»l'i 

•X.»UU(r*0A 

-li»nuK-n. 

•l,»4UE.rl 

•l.»oufc*oi 

-A  •  u.’yb^Ui, 

w  , 

i.uUoh+UA 

*.unut^  A 

J. uu.i »ui 

4.l'>'Ub»bl 

9*  00Ub«UA 

«,gngE*01 

/.gogE.'i 

y.gogt.oi 

LLH 

•A  » 31  uH*0x 

-i.OMuk-Ui 

•A.9'iub-ui 

•A.3IJU*  *i  A 

-l.9Utth.lil 

-1.9tUb*IU 

•l«9l(Ub*0A 

•i,»nuF.ni 

•l,94„F.fl 

VO  UK^I 
2*QUuut*  j<: 

0  M*.) 

1 ,  "  jnu1  ♦’! s* 

liM^AAu  k  I'eti  J 

J. 

nJ  irll 

V.UU'iv*  ♦  U*> 

HVU  11 >fi<o 

U. 

va  tFrb* 

0  . 

u 

n  il'fcUJ 

t 

b  i IN) 

■  a«i'  i'i-*n4 

L-  »t'enj 

HP  *  *  N  » 
^.uuug(.ng 

cm  tf>*G) 
91 

xpfi  <LE 1 
4t4g9uh.au 

-  - 

NHCl'  \  1  L  »  1 
A-*-batuLt  'i.4 

HI  *1  1  1 
•  1  •  U  U 1 '  «•  •  j 

3*  J  J  J  k  L*U1 

iJmI  tLH| 

J  • 

1 

Tb  iStbl 
.9Jllt,L«0i 

U 

•J , 

l«tC) 

Th  (CfcC) 
»i4u9gr»ei 

hP»  i.j'lu  4l«  iJj'.U  IHMSl*  .9'  0 
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Table  10 


Sample  Output  for  a  Typical  Computer  Run  with 
Cavity  Aerodynamics,  Precessing,  with  Fixed 
Format  (Continued) 


r I nt  (att> 
u , 

v!X  (Mi 
'J. 

'J  r  |  H  ) 
j  • 

UH  (Ml 

•  a 

HUT) 

>.UUul)fe»U4 

*  (►  T) 

u, 

*T  IH'il 
<.U1u«k«0? 

-  .  . 

H  l-lflil 

3 » 

•ja**4A 

U  • 

•<a  (dcu 
l«VU0Ufc*l'4 

HH|a  (ULu> 

‘1 . 

A]  inefSl 

u, 

M  IH-LKI 

♦•a  Ufeti) 

■J. 

w  *  •  J  fc  u  ) 

J  • 

rhl*«  Ut(J) 

’*  a 

1)  (HPn) 
AabUUOfc^OJ 

UT  (HA'i/SfeD 
•  1  aOo'OH-Ol 

L“  <Uj/»FC> 
Mi 

-  - 

*<»J  |frH*J 

2  •  j  J 

fl* 

u. 

til  UHh) 
Ja 

u  (►na) 

c.uQUPt^d* 

Y  (fPS) 

V. 

m  nrj) 

M. 

CftT 

•  1 . 

JU 

-A#  U/U‘»fc*ul 

a 

Ml) 

J.ub42e*ui 

P  (WAVSfcC) 

1  ,hoOfe-01 

<)  (xtJ/3Pt> 
•J,3?8«t*0« 

-  — - - 

i: « 

•2a.il.Lj  Jr-U2 

wA 

l*2y4'jf  ji 

V  (Llll 
•.‘fa  4<v4H-m 

A  (LB) 

iig4ut)c*yn 

LDf  <LH) 

D«*u  <  LM  > 
i,uiunk»oo 

.♦»  ( Li  > 

♦yj 

n  1  /  n 

4. 3ho  >c*ul 

ib 

ta'M^Vf.j* 

I  1 

y .cHiuc-ui 

ui  T 

•  fr  ,iv>4fc.n* 

u  1  a 

•/, 

^o!*  («<ciiJ 

4. 

wl. 

duu  Jii-d2 

CM 

1  a  2v4ot  - JA 

HHiV  (uH»j) 

U  a 

ha 

-i.5?o0t-U4 

<4 

\  t 

hV 

i  •  od  v*>c*U  1 

w  v 

•1  a*»yl4h-'Jl 

Ai.pHAl) 

A.OH/tJfJl 

XV 

n. 

I^MiVu 

— 

i  V  u 

C* 

Ui Art  4  A  J 
-A.ftVA4^-Ul 

1  l.lfc  (SPV,) 
a.JlJtt-Jl 

j*  CM  t 

S.  ) /-i /•-♦y* 

JT  IM  ) 
3.944?b«iM 

JH  If  1) 
•’>  a  //dOh«ilA 

HUT) 

♦  a  Wv4fc*u4 

*  (Ml 
4*v;-*7HAni 

*T  (>0a) 

l  lfe*i> 

4. Ii&jfc»u4 

•iA«.1A  (jfc*) 

Ha  CJc) 
•laU*VfJi 

HH I  a  ( utu ) 
AaUU07c«UU 

A  (  (  OfeG  / 

1 afeltOHA^U 

«T  OT.l.il 
•l  |3«t»*e-y? 

-*»  \  >C«i  > 

*alJJ*’»«  Jrf 

cA  (J*u> 

•4a } 

f'rtl*  UKi) 
t  a'M^/t^Jd 

•J  (HPM) 
i,jnu0c*y-5 

0\  (HM^/SfeT) 

•1 «  Ay 

(,>  (ntu/akt) 
U. 

( u  1 ► rb  1 
I'l/UfU 

fiJ  UHa> 
AiiUU4e»U? 

t  ►  Hj  ) 

'J.yMjH*  JU 

u  (fHa> 

i , 

y  (►P^) 

1 ,  U0u0H-H 

V  (Ira) 

9 , ♦Mt' 

C  '.T 

•  1  .  *0'/  y-»U  J 

J'J 

•1.  y 7o/**d1 

Vw 

*a  5/40K- Ji 

Mi 

C  ,  J'J*4C»U1 

P  (Hu'l/SfeC) 

1  l  4  U  Xpt-T,; 

13  <"»u/»  =  c> 

r  v 

3.  J  rf  J 

l  * 

Aa  41?'Jb*U1 

Uh) 

'a  30H*M 

A  (LB) 

1  a  02b  rife  ♦UO 

L 1 t  T  ( 

A,*54“OH-ni 

nx»u  (LB) 
l.ouifc.uo 

*  1  ( ) 
3.4,J43f*llJ 

rt !  f*\ 

B  •  0  <J  0  *  fc  •  'J  ] 

IS 

1  aU4H*f  '* 

1 1 

3 . 66/Vt-i)4 

11*1 

Mia 

*V,u?/»fc-0^ 

-oIm  i  ••BUI 
i  .  l34*fc-lU 

cL 

4  a  /1>1B*U<T 

l;i 

A  a  .Ifiolr-UA 

^  H | V  (Utg) 

A  aUU/lfe«UU 

ha 

1 » 4 u  !yfe-n^ 

»> 

n  * 

4  •  ?  3 1  >  t  •  y  j 

H  v 

A • 4 J  « 1»»UP 

w  V 

•1  a  4<h**-*l.H 

AlHHAU 

A  .  APft/fUl 

«V 

KMIVO 

"Vu 

^  .  4  4^0*- -IN 

uAdMa<. 

“A  .  <»>Vr  fc»Ul 

) 
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r 


line  (att) 
«.yUule-Ll 


ll.it  Ibtu) 
b.uUblfbl 


Table  10.  Sample  Output  for  a  Typical  Computer  Run  with 


Cavity  Aerodynamics, 
Format  (Continued) 


Precessing,  with  Fixed 


U*  IM) 

lit  |>  u 
9i4bk.'k«g4 

fct*  l»  1 » 
•J!.19l<U»UU 

MOD 

«.yy)«e»u.i 

I  lit) 
/ivmMni 

yt  ilkk) 
XiV)yfk*w? 

►  V  lUfcUl 

<i*mma  lytu) 
•i.»>blk*Ub 

Hk  irui) 
■1.9i,8kt*Ul 

KKIb  tutu) 
».k/14f(ll 

ki  (Utai 
i),Vb«gt»ng 

Ht  (fi-un 

l,l?4lt>ll 

hA  XUtbl 

L  *  <||ttl> 

•llMbi'kaiil 

KHU  ll'tti) 
4.4/14K-01 

0  (NPM) 

;,gngot»gj 

ul  (map/SkC) 
l.inA/t-m 

If M  IMAU/kK) 
U.‘ 

Xu  UHbl 

ru  <)K»I 
J./liUkrgZ 

/U  lllb) 
l.UUBkt'Ui 

U  IfKb) 
A.V*>t9c»llZ 

V  OKS) 
l.iyJAfc-01 

•*  l)Pb) 

or 

dll 

•*(Udj4»»Ul 

Vu 

•  7,/y.Ui-Jt 

■  0 

4ti)Zk0e*U) 

k  ta«o/s»Ci 

•l.UZA/6-01 

0  (M»g/b*t) 
1.4»)?fc.||4 

C.'i 

la 

A. 

N  (lU  1 
1.4yil9ca'jg 

*  I  La  i 
A»klu9c*il) 

LlM  »U*4 » 

I, Jy44t*ng 

nuti.  at) 
X.kA47t*0D 

mT  (lH) 
iob/»:*0d 

Kl/M 

A, 

is 

l,t|«64i-*')Z 

11 

A'k»A4e>g.l 

U|T 

•  /,kz»gl--09 

b|» 

-I,g4kle»g4 

P»IA  Uic'H 

wL 

l.UV/Sfbl 

to 

1. Jo7Vc-l)i 

KHJV  (utltl 

4.k7/kK*Ul 

KA 

•l.WA»/k*01 

•  A 

A.iiTAjfc.Q* 

n  A 

•St 9*hvtm JS 

l*V 

WV 

•1 .ZiUke-lIl 

ALKHAD 

l.414«e-gi 

My 

l,J4«/M-04 

KM  l  Vl) 

•i,i)»/»f(ll 

nVu 

4**A¥*k- J4 

uAhkai.' 

•l.Altkk'gl 

J»  (Ml 

1  ,  1  Ol)0**IJ« 

IT  |UI 

i,M  If  T  » 

•A.bbikt’ug 

M  (1  T) 
4,y»»4e»gA 

»  lit) 
l.lePkt»i)it 

yt  iFPk) 

x,y4u«e.o? 

H»  (UctiJ 
.4*17  VOfUJ 

li.l'n*  <t|td) 

•  y.gse(!kayu 

H»  (Ulj) 

1. 7ykot*l)« 

KM IS  lutul 
i.<A9?t»VU 

AT  (lltP) 
4,b»»yfc*nu 

mt 

Z,g7y7t-ui 

HA  (UfcbJ 

•  A. 

c»  tttu) 
v. j«ve k*U« 

Pm  1  a  (bell) 
•l.i’A9<rt»Uu 

1)  IMPMI 
lfbUgUc'g.* 

OT  ts.IVStC) 

L“  (M»i)/Sft) 

Xu  llPbl 

T l  IkKjl 

k.UT/lfU^ 

in  (I  Kb  l 
1.4/Zot-Ui 

n  ()K») 

i.yjplfbi 

y  iri><?> 

•*,  )9ws*n 

i  DPb) 
i.A7A»t*(l 

C.1T 

Ul 

«u 

•A.ZVI)!,)  »l>ll 

■II 

1 ,6Jobe»bl 

K  (KA"7btr) 
•i,MgJkK-ni 

g  (MAu/kFti 
k,»t/kt-u4 

r.> 

.  l.ifle^fcjLux 

L  r 

*.4*eS*»Ul 

>  igHl 

K.lil*K*Ug 

•  (la> 
i./y»7t»ui 

tllT  ILM) 

1 ,V/Okfc.0g 

Data  (i  al 
I,y4g7e«00 

•a  1  (  lh  J 

•* .  duuuh*ug 

M 1  /  K 

lb 

1 .1’kkgt  *'U 

1 1 

9.kJ/4e»g.» 

HIT 

•B,A4bkK-n9 

bib 

-I,j4yit*ii4 

^b|M  UJtlii 

bl 

i.-C<uV«-yl 

UJ 

pH  I  V  Igfci,) 
*i.4AkAt»OI 

•  i  ,ymot-ni 

■  A 

4,k.1A7t»U4 

•  1.  t  *i  r*  4  bm  >14 

MV 

•*.y  ivls-gl 

gv 

*i ,gj«5a-ui 

«LKM« J 
l.uiBJe- Jl 

■  y 

• J , y /*yt*n j 

KMIVU 

-i,yi»u*UT 

*  *  w 


nAiV'nl' 

■a.  JdMfc-Ul 


I 
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Table  10 


Sample  Output  for  a  Typical  Computer  Run  with 
Cavity  Aerodynamics,  Precessing,  with  Fixed 
Format  (Concluded) 


1 1  1-  (j  J 
o  >  J  J  J  U  -  1 

J «  lr  1 1 
l  .  ♦  J< 

lH  (II) 

yN  I r  T 1 
•  / .  *Mo**Uy 

him 

«#y«<4t*u4 

r  f>  T  | 

j  ,pel>«»fr*r* 

»»  f»r,» 

•  *  1 CC 1 1 
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Table  11.  Typical  Input  for  a  Simulation  Run  with  Fixed 
Aerodynamics 


*W  "“V,'""  C,L  •  •  PStCESSIKO  -  l  LB  THHUST 


HP 

2. 

HP 

3. 

hhop 

2.06 

OMfcUA 

1500. 

V 

50. 

H 

3000. 

UA 

0. 

theta 

0  . 

PHI 

0. 

HF 

-1.0 

HV 

U  .0 

VA 

0.0 

PSI 

u .  0 

HC 

■j  .  0 

UH I 

o.o 

T  d 

15. 

T  1 

u .  0 

TF 

15. 

innusi 

1. 

KP 

1. 

KM 

u . 

TINT1 

.2 

TINT2 

1. 

TSTOP 

2. 

NCURVE 

1. 

-10. 

-3.  0. 

-0.40 

”0.24  -0,02 

-10. 

-5.  0 . 

-.008 

r.on  -.004 

-10. 

-5.  o. 

0.20 

0.17  0.16 

CMQ 

-10. 

0*  10, 

-15. 

"15.  -15. 

CLP 

-10. 

0.  10, 

-15. 

”15.  -15. 

5, 

0.24 

5. 

.014 

5. 

0.17 

20. 

-15. 

20. 

•15. 


10* 

0.44 

10* 

-.109 

10* 

0.21 

15. 

0.48 

15. 

-.008 

15. 

0.28 

20. 

0.60 

20. 

.012 

20. 

0.30 

40. 

1.20 

40. 

.076 

40. 

0.23 

60. 

1.37 

60. 

.062 

60. 

0.11 

90. 

1.34 

90. 

0. 

90. 

0.0 

* 

30. 

-15. 

40. 

-15. 

50. 

-15. 

60. 

•15. 

70. 

-15. 

90. 

-15. 

30* 

-11. 

40. 

-15. 

50. 

”15. 

60. 

-15. 

70. 

-15. 

90. 

-15. 
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Table  12 


Sample  Output  for  a  Simulation  Run  with 
Fixed  Aerodynamics 


•  •  INPUT  DATA  •  • 


NUN  1<A*1  CONI  IQ  1A  -  4 

10  1  CVL  -  F 1 460  AERO  - 

PRECISE  1 NC  • 

1  Id  THRUST 

ALPHA 

-1. J006*ul 

•5.0006*00 

J. 

5.0006*00 

1 . 000E ♦ 01 

1.5006*01 

2.0006*01 

4 .0006*01 

6.0006*01 

9.0006*01 

CN 

*4 .  uOOE-ul 

-2.400E-U1 

*4 . 000E*tl2 

2.4006-01 

4.400E-01 

4.00QE-01 

O.OOOE-Ol 

1.2006*00 

1.3706*00 

1 . 340E- or 

AlHhm 

-i.  Juoe*ui 

•5. 0006*00 

u . 

5.0006*00 

1.0006*01 

1.5006*01 

2.0006*01 

4 . 000E*01 

4.000E*01 

9.0006*01 

CM 

•».  iOOE-gJ 

-1.100E-U2 

-4.0006-03 

1.J00E-02 

•9.0006-03 

■t . 000E*03 

1.2006-02 

7 .40116-02 

6.2006-02 

0. 

AlHhA 

•1 . JU0E»ul 

-5,000E*U0 

J  • 

6.0006*00 

1.0006*01 

1.5006*01 

2.0006*01 

4 . UOOE*01 

6.0006*01 

9.0006*0] 

CA 

i. oUOE-ul 

1.700E-U1 

1 .*006-01 

1.7006-01 

2.100E-01 

2.1006-01 

3.000E-01 

2.300E-01 

1.1006-rl 

0. 

ALPHA 

•1. JOOE’Jl 

0. 

1 .0006*01 

2.0006*01 

3.000E*01 

4 • OOOE*01 

5 .0006*01 

4.0006*01 

7.0006*01 

9.0006*01 

CMU 

-1.5U0E*J1 

-i.5oot*ui 

*1.0006*01 

•1.500E*01 

-1.5006*01 

•1 . 50UE*01 

•1.500E*01 

•1.5006*01 

•1,500F*01 

•1.5006*01 

ALPHA 

-1.u00E*j1 

0. 

1.0006*01 

2.0006*01 

3.000E*01 

4 .0006*01 

5.0006*01 

4 .0006*01 

7,0006*01 

9.0006*01 

CLP 

-1.5U06*.l 

•1.5006*01 

-1.5006*01 

*1.500E*01 

-1.5006*01 

•1.500E»01 

-1.500E«01 

•1.5006*01  -1.5006*01 

•1.5006*01 

VO  ( FPS  > 

QAHMAU 

OfcU) 

HVO  (U6GI 

H  (066) 

6 

(DEO) 

PHI  (PEG) 

:  ,00uu6*0l 

0. 

0. 

0 

• 

0. 

0. 

0  (PPM) 

HO 

(1  T) 

VA  (FPS) 

0  (IN) 

HP 

(IN) 

PPO  (1.6*  > 

t .50Ug6*0J 

5.00106*03 

0. 

8 

.00006*00 

2.0000E*00 

5.00006*00 

«HOP  l DEN  1 

PC 

(LSI 

DPI  (LB) 

Td  (SEC) 

T 1 

(SEC) 

Tl  (SEC) 

<.06UU6*0j 

0. 

0. 

1 

•5000E*D1 

0. 

1.50006*01 

HI  (FT  1 

-1.00UG6*0  ' 


*P<  1....0  KM>  O.OUu  IHMUSTi  1.000 


) 
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Table  12.  Sample  Output  for  a  Simulation  Run  with 
Fixed  Aerodynamics  (Continued) 


nub  <sk) 

0. 

UX  (M) 

0. 

Dt  (FI) 

0. 

OH  (FT) 

-0. 

H  (FT) 
».0000k*U3 

X  (FT) 

0. 

VT  ()»S) 
9.00002.03 

•* 

MV  (DEG) 

0. 

GAMMA  (BEG) 

0. 

MS  (DEG) 
l.««00E*02 

PHIS  (DEG) 

0. 

AT  (DEG) 

0. 

MT  (FT-LB) 
-2.3970E-03 

HA  (l)2G> 

0. 

EA  (BEG) 

0. 

PH  I  A  (DEGJ 

0. 

0  (KPH) 

1 .  5000b  +  U3 

01  (RAD/SEC) 
-3.A309E-03 

DM  (R.U/SFC) 

0. 

XD  (FPS) 

9 . O00DE*01 

TO  (FPS) 

0. 

ZU  (FPS) 

0. 

U  (FPS) 
9.0000E»U1 

V  (FPS) 

0. 

M  (FPS) 

0. 

CMT 

-4.0000E-03 

UD 

-9.102*1-01 

VD 

0. 

HD 

2.Hm*01 

P  (RAD/SEC) 
1.A309E-03 

0  (RAU/SEC) 

•* • A429E-09 

CM 

-2 . OOOOE-02 

CA 

3. *0002-03 

N  ( LB  ) 

*1 . 7A04F-02 

A  (LB) 

3 . 4147E-U1 

LIFT  (LB) 
-1.7»*4E-02 

DRAG  (lB) 
3.4147E-01 

MT  (LB) 
9.0*0*E«00 

Ml/R 

2.S6G4E-01 

IS 

9.203AE-03 

IT 

4 .96392-03 

DIT 

-A.1993E-09 

DIS 

-7 . 993*E-09 

PS | A  IDES) 

0. 

CL 

-2. 00**2-02 

CD 

3.4000E-03 

PM | V  (DEU) 

0. 

PA 

1.6309E-03 

at 

-A.A429E-09 

HA 

0. 

PV 

1. *3091-03 

ov 

-6.4973E-03 

ALPMAD 
* . 49*6E-U1 

RV 

0. 

pmivd 

3.A309E-03 

MVS 

0. 

GAMMAD 
-  A . 4973E-01 

• 

TIME  <$Et> 
2.Uo01E-U3 

DX  (FT) 
9.9*l*E-00 

DY  (FT) 
-4.9322I-0A 

DM  (FT) 
-9.04f*E-03 

H  (FT) 

4 . 99 9 9E*03 

X  (FT) 

9.9*3*E*00 

VT  (FPS) 
9,00992*03 

MV  (DEC) 
-1.9091E-04 

GAMMA  (BEG) 
•9.72*9E*00 

MS  (DEG) 
3.77ttE*02 

PHIS  (DEU) 
2.9379E-02 

AT  (DEG) 
9.7220E*00 

MT  (FT-lB) 

S. *0422-03 

HA  (DEG) 
-2.73342-03 

EA  (DEB) 

1 • 02S3E-03 

PHIA  (DEG) 
-2.9390E-02 

0  (HPM) 
3.9000E-03 

OT  (RAD/SEC) 
4.0223E-03 

OM  (HAU/SEC) 

0. 

XU  (FPS) 
4,«Bias«oi 

Y0  (FPS) 

•1. *9*21-04 

20  (FPS) 

4.9903E*00 

U  (FPS) 
4.9S3»E*01 

V  (FPS) 
-2.19A9E-03 

M  (FPS) 
4.9909*. 00 

1 

CMT 

9. *2332-03 

UD 

-3 . 0374t*00 

VD 

-3.4393E-02 

MO 

2 . 43092*03 

P  (RAD/SEC) 

•4. *2231-03 

a  (Rao/sec) 

1 «»44*E-94 

I 

CM 

2.IIME-01 

CA 

1.797*1-03 

N  IlB) 
2.3S30E-01 

A  (LB) 
3.9979E-01 

LIFT  (LB) 
2.239SE-01 

DRAG  (lB) 
3.7*772-03 

1 

MT  (LB) 
4.933*E»00 

HI/A 

2 • S003E-03 

IS 

*.3**91-03 

IT 

4.9404E-03 

DIT 

-».*323E-(9 

DIS 

•*. 027*2-09 

| 

PS | A  (DEO) 
-2.71112-03 

CL 

2.90011-03 

CD 

2.0373E-03 

PM I V  (DEG) 
-2.9204E-02 

PA 

-4.02231-03 

Ga 

3, *42*2-04 

f 

HA 

-6,A6A«-04 

PV 

•4. G709E-03 

SV 

•* . 10A2E-01 

ALPMAD 

A.3079E-01 

RV 

-2.0220E-04 

PmIVU 

-4.0*922-03 

I 

l 

MVO 

-1.2*022-09 

GAMMAD 
•» . 30»2E-01 

j 
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Table  12.  Sample  Output  for  a  Simulation  Run  with 
Fixed  Aerodynamics  (Continued) 

* 

Tint  (SM) 

4.  ille-  1 

OX  (Ml 

1. 99226*1)1 

UT  (III 
-1.24276-04 

UH  (ft) 
-1.97566*00 

H  (FI) 

4 . 99606*03 

X  (IT) 
1.9V22E*01 

VT  UPS) 

5.05126*01 

HU  (Dili) 
-1.2Q616-U1 

GAMMA  (  UkG ) 

•1 . 10436*01 

HS  (UtU) 
1.77616*02 

PHIS  (UEU  > 
1.92466-02 

*1  (DtG) 

1 . 10446*01 

HT  (fT-LBl 
-5.2.916-03 

H  A  ( DIG) 
-4.959/6-U1 

e*  IU60) 

7 . 77566-04 

Ph|*  (DeS) 
-1.92326-02 

0  <NPM> 

1 . 50006*u3 

OT  (PiP/SEC) 
-3.67296-01 

OH  ( NiU/SFC ) 

«• 

XO  (EPS) 
4.997/6*1)1 

fD  UPS) 
-1.04306-03 

21)  (fPS) 
9.67566*00 

U  UPS) 
4.95776*01 

V  (TPS) 
-4.29146-01 

*  (IPS) 

9.67636*00 

CUT 

-6  •  79126  -  01 

uo 

-1.33906*00 

VI) 

2.46576-02 

nD 

2.20796*01 

P  ( 8AD/S6C ) 
3.67296-01 

U  (  NiO/SFC ) 

-1.52026-84 

CM 

4.48156-01 

CA 

2.24626-01 

N  <LH) 
4.04616-01 

*  (Lb) 
2.0270E-01 

Lilt  (LB) 
3.5I29E-01 

DHil  (LB> 

2.76466-01 

«T  (LH) 
4.866/6*00 

Hl/H 

3.01466-01 

15 

9.16746-01 

II 

4.93396-03 

BIT 

-7 .62906-05 

0 1  a 

-1.04616-04 

PSI4  (D6G) 
-4.99966-01 

CL 

1.97026-01 

CD 

3.06346-01 

PH  1 V  ( U6G ) 

•1 .95956-02 

PA 

3.0/296-03 

0* 

-1.51706-04 

HA 

4,99056-1)4 

PV 

3. /0056-03 

UV 

-5.76246-01 

AlPHAD 

5. 79096-Ul 

BV 

-2.13796-04 

Pm  I9U 

1.75366-03 

MVU 

-1.03396-05 

(iAHHAC 

-b.;e246-oi 

tint  (Si..) 
t.  016-  1 

OX  (IT) 
2.96076*01 

Ct  (Ml 
-1.49626-04 

DH  (M  ) 
*4.36546*00 

H  (f  1  ) 

4.99566*03 

X  (M  I 
2.98076*01 

vt  urs> 

5 . 1 ?09t  *  0 1 

M  V  (06UI 
-1 ,19706-UJ 

Gamma  iu6o) 
-1.60646*01 

MS  (DtG) 
-2.4525M0U 

PHIS  IU6u) 
1.76276-02 

tl  (DtG) 
1.60816*01 

hi  dt-Lb) 

-2.2725C-03 

' 

Hi  (  DeU  I 
1.92066-01 

6*  (1160) 
-7.95906-04 

PHI*  (1)6G> 
1.76116-02 

0  (HPH) 

1 . 50006*1)3 

01  (HAD/S6C ) 
-1.58196-93 

CM  (MAD/3PC) 

0. 

X  U  CAPS) 
4.9202b*Dl 

r D  UP*) 
-V.940Ofc-Q4 

2D  <FPS> 
1.42(146*01 

11  (IPs) 
4.92626*01 

V  (fPS) 
1.1/39E-01 

9  (IPS) 

1.42036.01 

CHT 

-i . 6662b ” 0 3 

Ul) 

-1.  77Q56*I)0 

VI) 

3.22366-02 

«D 

2.23136*01 

P  (HAO/SfeC) 
1.50196-03 

3  (P.U/SFC) 

-6.63696-05 

O 

5.00006-Ul 

VIA 

2.04336-01 

N  <Lfi) 
4.70466-01 

A  do) 
2.04366-01 

LHT  <L«> 

3 . 7080F -n 1 

Dm»u  (|  6  ) 

3.04156-01 

hT  (LH) 
4.80006*00 

h|/H 

3.2144fc-01 

IS 

9.14516-03 

1 1 

4 . 91 79t-03 

HIT 

-8.34586-03 

CIS 

-1.18956-04 

Pal*  M)6G> 

3. 92076-01 

CL 

4, 07436 -Ul 

CO 

4. 13396*01 

PH  1 V  (OEu) 
1.03206-02 

PA 

1. 50196-03 

U  A 

-6.0477fc.flS 

1 

HA 

6.5684fc-U4 

PV 

1. 70206-03 

UV 

-5.33446-01 

ALPHAD 

5.53376-01 

HV 

1.92866-04 

PmJVC 

1,09/46-03 

1 

I 

H  V  U 

1  .  A4P7P-fl,> 

uammali 

-5 . 5344fc -  01 

- 

1 

1 

■ 

1 
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Table  12 


Sample  Output  for  a  Simulation  Run  with 
Fixed  Aerodynamics  (Concluded) 


I  I  He  (Sti.) 

a.  oit-ii 


l|HE  (SE«) 
1 . <«00fc<uU 


UX  (Ft) 

3<9<22E<01 

DT  (FT) 

*4.9«7k<04 

DH  O  T) 
<7<04S9E<00 

H  (F  1  ) 
4<V<24E<03 

X  DTI 
3<9<22F<01 

VT  (TPS) 
9 . 2?44e<  01 

HV  (DkG) 
-7.4so;e-o4 

GAHHA  (DEG) 

•2  <  O700k<0 1 

HS  (DkG) 
1<77S<E<02 

PHIS  (DEG) 
1<190<E<02 

A)  (DkG) 
2<079<E<01 

HT  (FT-lB) 
9.3A22E-03 

HA  (DkG) 
<9<2<29E<U3 

fc*  (DEG) 

9  <  0993k*04 

PHIA  (DkG) 
<1<1<9<E<02 

0  (RPH) 
1<900<E<03 

OT  (HAD/SEC' 
<<9493E<03 

DH  (RtU/SFC) 
0. 

XU  (FPS) 
4<<887E<U1 

YD  OPS) 

•O  <  39  72k  -  04 

20  (FPS) 
1<A9<<E<01 

U  OPS) 
4<8<<7E<01 

V  (FPS) 
-4.4904F-03 

W  (TPS) 
1<B9<7E<01 

CRT 
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PA 
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OA 
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H* 
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PV 
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ALHHAD 
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RV 
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PMIVD 
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DY  (FT) 
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DH  (FT) 
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X  (FT) 
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VT  (TPS) 
9.3943k<0l 
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-3.9329E-03 
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HS  (DEG) 
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PHIS  (UEG) 
1<2<1BE<01 
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CN 
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A  (LB) 
2<»24fc<01 
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IS 
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CL 
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CD 

9<74<9E<01 
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PA 
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HA 
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PV 
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QV 
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RV 
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1.90Mk-04 

GAHHAD 
•  4 . 7  U03k <01 
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Another  useful  modification  of  the  program  to  permit  simplification  of  the 
input  data  for  nonburning  applications  (inert  or  ordnance  items)  would  be  to 
permit  input  of  total  weight  and  inertias  and  bypass  the  separate  case, 
pyrotechnic  parameters,  ignitioa  and  burning  times  which  have  to  be  care* 
fully  "faked  out"  to  simulate  a  solid  body  and  avoid  singularities  in  the 
automatic  calculations  of  inertias.  If  it  is  intended  to  alternate  between 
inert  and  burning  configurations,  this  could  be  handled  by  input  logic  which 
demand  only  the  inputs  needed  at  the  time. 

If  it  is  intended  to  use  the  program  for  nonburning  ordnance  items  then  it 
would  be  preferable  to  create  a  separate  deck  with  the  unneeded  variable 
mass  and  aerodynamic  terms  deleted.  This  latter  form  would  be  directly 
applicable  to  flettner  rotor  munitions  for  which  presently  available  pro¬ 
grams  are  either  nonprecessing  3  degrees  of  freedom,  using  planar 
aerodynamics,  or  full  6  degrees  of  freedom  requiring  the  full  set  of  six- 
component  static  and  six-component  dynamic  aerodynamics. 
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APPENDIX  A 

PRELIMINARY  COMBUSTION  WAKE  EFFECT  ANALYSIS 

The  following  analysis  of  the  aerodynamics  of  the  combustion  w*ke  was  a 
preliminary  study,  based  on  preliminary  estimates  at  NAD/Crane  of  com¬ 
bustion  pressures  and  temperatures.  The  predicted  pressures  were  not 
sufficient  to  produce  thrust  so  that  the  results  of  this  analysis  are  not 
applicable  to  thrust  producing  flares. 

1.  INTRODUCTION 

The  effect  of  the  combustion  wake  on  the  aerodynamic  characteristics  of  a 
spinning  disk-shaped  flare  will  be  determined  by  examining  the  combination 
product-air  mixture  in  the  flare  plume.  The  flare  considered  in  this 
analysis  is  assumed  to  burn  only  on  cavity  surfaces  (see  Figure  67 ).  Then 
by  symmetry  arguments  and  also  since  a  flare  is  usually  designed  so  that 
the  pyrotechnic  burns  at  atmospheric  pressure  (i.  e. ,  no  pressure  jump  in 
the  flare's  aerodynamic  behavior  is  caused  by  a  change  in  density  of  the 
medium  near  the  flare),  any  finite  combination  product  expulsion  velocity 
will  not  produce  a  thrust  force.  Thi*  change  in  density  results  from  the 
introduction  of  combustion  products  into  the  airstream  and  from  the  increase 
in  temperature  of  the  air  surrounding  the  flare.  Using  this  change  in  density 
approach,  upper  bounds  will  be  obtained  which  give  the  maximum  possible 
change  in  the  flare's  aerodynamic  coefficients.  Then  the  effect  of  the  flare 
wake  will  be  negligible  throughout  the  flare's  flight. 

2.  THEORETICAL 

The  mass  flux  of  combustion  products  into  the  flare  field  is  constant  if  a 
grain  configuration  of  constant  burning  area  is  used  and  if  the  pyrotechnic's 
burning  rate  is  insensitive  to  the  flare's  orientation  and  velocity.  Therefore, 
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Figure  67.  Flare.  Preliminary  Combustion  Wake 
Effect  Analysis 
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where  m^  is  the  combustion  product  mass  flux  into  the  flow  field  and  mp 
is  the  pyrotechnic  mass  consumed  in  tp  seconds. 

If  a  mass  balance  for  the  control  volume  (see  Figure  67,  exclusive  of  flare 
volume)  is  considered,  the  density  of  the  combustion  products  in  the  vicinity 
of  the  flare  can  be  determined.  For  the  control  volume  V. 


m  =  mJN  -  mQUT 


m„ 


B 


S' 


Op  V  *  ds. 


(2) 


where  v  is  the  velocity  of  the  particles  as  they  pass  through  S',  the  boundary 
of  V.  However,  at  any  given  time,  the  mass  of  the  combustion  products  in 
V  is 


Ppdv.  (3) 

V 

since  pp  may  not  be  constant  in  V. 

Therefore, 


m  =  Ht  °Pd  V)' 

so  that 


(4) 
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As  a  first  approximation  to  the  density  of  air  in  V,  the  ideal  gas  law  will  be 
used.  Therefore, 


P 


.  pamc 

A‘X*c“* 


(5) 


where  P^,  Mc#  and  Tc  may  be  functions  of  position  in  V  and  time.  Explicit 
determination  of  the  density  distribution  in  V  is  beyond  the  scope  of  this 
report.  However,  upper  and  lower  bounds  can  be  placed  on  p^  since  T^,  and 
Mc  are  bounded. 

3.  NUMERICAL  SOLUTION 

The  maximum  change  in  the  normal  force  coefficient  will  occur  when  the 
flare  is  at  90  degree  angle  of  attack,  while  the  axial  force  coefficient  and 
center  of  pressure  location  will  change  most  when  the  flare  is  at  0  degrees. 
Therefore,  combining  these  maximum  changes  will  give  the  greatest  effect 
of  the  combination  product  wake  on  the  aerodynamic  behavior  of  the  flare. 

If  pp  and  p^  are  to  be  obtained  from  equations  4  and  5,  an  appropriate 
control  volume  must  be  determined. 

The  pyrotechnic  considered  in  this  problem  burns  at  atmospheric  pressure 
so  that  the  change  in  pressure  force  is  small  compared  to  the  buoyancy 
force  or  change  in  skin  friction  coefficient.  Then,  for  maximizing  purposes, 
the  control  volume  of  equation  4  can  be  taken  as  the  flare  cavity  at  burnout 
so  that  S'  =  S  where  S  is  the  flare's  projected  area. 

Then,  with  V  corresponding  to  the  flare  cavity,  v  will  be  the  velocity  of  the 
combustion  products  leaving  the  cavity  with  component  vn  normal  to  S. 
Although  Pp  is  a  function  of  position  in  V,  it  will  be  assumed  for  simplicity 
that  Pp  is  a  function  of  time  only.  Then  equation  4,  which  reduces  to 
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(6) 


p_  v 
KP  n 


S, 


gives  an  equation  for  the  average  density  in  V.  Solving  for 
Op  at  t  =  tg  with  Op  =  0  at  t  =  0, 


°P 


mP/tB 

ITS- 
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1  -  exp 


"vn  hi 
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(7) 


4.  UPPER  BOUNDS  ON  FLARE  WAKE  EFFECTS 

a.  Buoyancy  Force 

An  increase  in  temperature  of  the  air  in  the  flare  cavity  adds  buoyancy  to 
the  flare.  The  maximum  buoyancy  force,  occurring  when  the  flare  is  at 
zero  degree  attitude,  can  be  determined  by  placing  limits  on  the  average 
temperature  of  the  air  in  the  cavity.  These  limits,  given  by  equation  8, 
although  simple  in  form,  give  an  upper  bound  on  the  buoyancy  force  for  all 
flare  orientations, 

ta  4  Tc  ‘  tb-  <8> 


where  T^  is  free  stream  static  temperature  and  Tg  is  the  pyrotechnic  flame 
temperature.  For  a  stoichiometric  mixture  of  pyrotechnic  and  oxidizer 
(02),  assuming  all  the  0^  in  the  cavity  is  used  during  reaction,  the  molecular 
weight  of  the  air  in  the  cavity  will  be  the  molecular  weight  of  the  inert  gases 
in  air.  This  condition,  combined  with  equation  8,  gives 


PAMa'  s  .  ,  s.  AMA 

“  A  tttj- 


TTT 


(9) 
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where  =  Molecular  weight  of  air 
*  Molecular  weight  of  Ng 
o  ^  =  Density  of  air  in  the  cavity. 

Therefore,  the  maximum  buoyancy  force  on  the  flare  will  be 


A' Is  h 

B  J 


» 


(10) 


which  assumes  that  the  thermal  conductivity  of  the  flare  case  is  zero,  i.  e. , 
the  flare  case  is  a  perfect  insulator. 


Bf  = 


"PA  MA 


PA  MA 


]s  h  g 


=  (2.38  x  10"3  -  .17  x  10“3)  Shg 
S  =  nd2  =  -J- 

g  =  32  S  h  g  =  4.  19 

h  =  2/12  ' 

=  2.21  x  10"3  (4.  19) 


Bp  =  9.  25  x  10"3  pounds  «  Weight 
Bp  <  <  weight  of  flare. 
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b.  Skin  Friction  Coefficient 


The  skin  friction  force  for  a  nonburning  flare  is  defined  by 


=  cF  qs  =  cF  d/2  p  os. 

NB  fnb  fnb  nb 


For  a  burning  flare,  the  effective  density  will  not  be  the  same  as  free 
stream.  Therefore, 


FB  =  CFB  pb(1/2v  1  s 


~  CF, 


[^r](1/2  “nb  v2>  s- 


since  Cj,  »  CFNB  *  Cp  and  defining  Cp'  =  c FH [  ,  ]  CF  [  n,R  ] 


l!s.  I 

r  °B  i 

L°NB  J 

MIN.  *  CF'  *  Cp  J 

MAX. 

However,  P  B  =  °p  +  ’ 


where  Pp  is  the  pyrotechnic  combustion  product  in  the  control  volume 
considered  in  o^'  is  atmospheric  density  in  this  volume.  If  6p  =  0  in  the 
cavity,  then 


rn  T  =  m 
IN  OUT, 


where 


?  =  /S°PvndS 


=  o'  V  S, 
P  n 
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where  Op'  is  the  density  in  this  volume  when  Pp  is  zero 
in  V').  Let  the  time  equal  t  ',  then 
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Mp/  tg 
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Substituting  p 
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which  is  only  true  when  t'  =  Therefore,  when  pp  =  0, 
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=  Const. 


(equilibrium  density 


Now,  since 
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For  this  problem, 

Ma'  =  22. 1 

Ma  =  28.  9 

T.  -  273 
A 

Tfi  =  2973. 


Therefore,  for  vn  «  10  feet  per  second, 

C  ' 

0.070  *  «  1.67. 

However,  the  skin  friction  coefficient  is  Cp  <0.  005 >  which  is  substantially 
less  than  the  axial  or  normal  force  coefficients. 

5.  CONCLUSIONS 

If  thrust  is  not  generated  by  the  burning  pyrotechnic,  the  effect  of  the  flare 
plume  on  the  aerodynamic  characteristics  will  be  negligible. 
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APPENDIX  B 


EQUATIONS  FOR  TRAJECTORY  SIMULATION 

The  following  equations  for  quasisteady  precession  and  translation  of  the 
flare  under  the  assumed  aerodynamic  torques  and  forces  were  programmed 
in  FORTRAN  IV  for  the  CDC  6600  computer  and  used  to  compute  most  of 
the  trajectories  presented  in  this  report.  This  program  permits  selection 
of  either  free  precession  or  suppression  of  precession.  Some  of  the  non- 
precessing  trajectories  shown  in  the  report  were  computed  with  a  truncated 
fixed  orientation  version  of  the  program  on  a  teletype  computer  in  a  special 
language  and  are  not  presented  here. 

The  order  of  presentation  is  as  follows : 

.  Axes 

.  Translation  equations.  * 

.  Aerodynamic  angles.  * 

.  Aerodynamic  forces  and  moments.  * 

.  Inertial  and  thrust  terms.  * 

.  Angular  rates.  * 

.  Axis  transformations  using  euler  parameters.  * 

.  Auxiliary  parameter  transformations.  * 

.  Initial  conditions  /  input  constants. 

.  Nomenclature.  * 


*  Reprinted  from  original  Honeywell  Inc.  Scientific  Programming 
Department  report. 
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1. 


AXES 


The  axes  may  be  defined  as  follows  -  all  orthogonal  triads  with  origin  at 
the  flare  center  of  gravity  except  for  inertial  axes  which  would  have  an 
arbitrary  initial  position,  (see  Figures  68-71). 

a.  Inertial  Axes 


The  Xi  axis  is  horizontal  parallel  to  the  arbitrary  reference  trajectory  plane 
(possibly  aircraft  direction  at  launch,  etc. ),  the  Z.  axis  is  vertical  down¬ 
ward,  and  the  Y.  axis  is  horizontal,  forming  the  other  axis  of  a  right  hand 
triad.  These  are  all  with  respect  to  a  flat  nonrotating  earth. 

b.  Flare  Axes 

The  Zp  axis  is  coincident  with  the  positive  spin  axis,  Xp  and  Yp  completing 
the  right  hand  triad  in  the  disk  plane.  The  component  of  angular  rate  of 
the  axes  themselves  along  the  spin  axis,  is  constrained  to  be  zero  (yaw 
rate,  R=0).  The  component  of  the  angular  rate  of  the  body  along  the  spin 
axis  is  then  the  spin  rate  ui .  For  zero  orientation  Euler  angles,  the 
flare  axes  would  be  parallel  to  the  inertial  axes. 

c.  Aerodynamic  Axes 

The  Z^  axis  is  coincident  with  the  positive  spin  axis,  X^  and  Y^  are  in  the  disk 
plane,  with  the  X^  axis  at  the  intersection  of  the  disk  plane  and  the  aero¬ 
dynamic  plane,  which  is  through  the  spin  axis  and  the  velocity  vector  per¬ 
pendicular  to  the  disk  plane.  The  normal  force  N  is  thus  by  definition 
parallel  to  and  opposite  the  Z^  axis,  and  the  axial  force  A  is  parallel  to  but 
opposite  the  X^  axis.  The  aerodynamic  axes  are  related  to  the  flare  axes 
by  a  yaw  rotation  about  the  Zp,  Z  axis. 
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Figure  68.  Aerodynamic  Axes,  Forces,  and  Moments 


MOVING  AXES 


(  )F  FLARE 

(  >A  AERODYNAMIC 

(  )v  VELOCITY 


Figure  70.  Euler  Angles  Between  Inertial  Axes 
and  Moving  Axes 
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Aerodynamic  Axes 


d. 


Velocity  Axes 


is  coincident  with  the  aerodynamic  Y  axis,  Y^>  or  pitch  axis,  X  ,  Z  ,,  in 
the  aerodynamic  plane,  with  Xy  aligned  with  the  velocity  vector  and  Zv 
forming  the  third  axis  of  the  orthogoral  triad.  The  lift  force,  L,  is  thus  by 
definition  parallel  to  but  opposite  the  Zv  axis,  and  the  drag,  D,  is  parallel 
to  but  opposite  the  Xv  axis.  The  velocity  axes  are  related  to  the  aerodynamic 
axes  by  a  pitch  rotation  through  the  angle  of  attack,  a^.,  about  the  Yft,  Yv 
axis. 

Each  of  the  axes  moving  with  the  flare  is  related  to  the  inertial  axes  by  the 
Euler  angles  H,  E,  0,  for  heading  (azimuth  of  X  axis),  elevation  (of  X  axis) 
and  roll  (about  the  X  axis).  For  flare  axes,  (  )^  denotes  aerodynamic  axes 
and  (  )  denotes  velocity  axes.  0C  and  H  are  used  to  denote  the  tilt 

of  the  spin  axis  from  the  vertical  and  the  azimuth  of  the  axis  about  which  it 
is  tilted,  respectively. 
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2. 

TRANSLATION  EQUATIONS 

U 

m 

-  /A\  cos  i |u+  g  a31  -  Q  W 

Axial  Acceleration 

• 

W 

V 

m 

-  /A\  sin  g  a32  +  P  W 

Lateral  Acceleration 

w 

w 

u 

-  N  ♦  g  «33  +  QU-PV-T 

Normal  Acceleration 

m  f  .  m 

u 

s 

U0  +  /  U  dt 

Axial  velocity  components 

V 

s 

V0  +  /  V  dt 

Lateral  " 

w 

X 

W0  +  j  W  dt 

Normal 

VT 

s 

Vu"  +  vJ  +  Wz 

Total  Velocity 

*0 

s 

VT  COS  Y  COS  Hu 
•o  0  vo 

Initial  velocity  components 

a 

Vjo  COS  Y0  COS  Hyo 

*0 

« 

-  Vt  slriY 
'o  'o 

• 

X 

s 

U  a,!  +  V  a12  +  W  a15 

• 

Y 

a 

U  i2i  *  V  8jj  +  W  a2 j 

Z 

s 

»3*  *  “ 

X 

■ 

ft  dt 

Down  range 

Y 

■ 

f  * 

Cross  range 

h 

■ 

h  -  Ji  dt 

Altitude 

AX 

* 

/ (v,  -  *>  « 

Position  relative  to  launch 

AY 

m 

Y 

Ah 

m 

-  y  z  dt 

q 

m 

%  PA  VT2 

Dynamic  pressure 
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3. 


AERODYNAMIC  ANGLES 


W 


'll2  +  V* 


tan"1  V 


where  Cj  ■  2.378  x  10* 

C,  -  3.2  x  10"* 


Total  angle  of  attack 


Yaw  of  aerodynamic  plane  with 
respect  to  X  flare  axis 

Atmospheric  density 
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4.  AERODYNAMIC  FORCES  AND  MOMENTS 


A  ■  q  S  CA 

/  \ 

Axial  force  parallel  and  opposite 
to  aerodynamic  X  axis 

cn  "  f  rc/r  J  * 

Seo  Input  Data 

N  -  q  S  CN 

/  \ 

Normal  force  parallel  and  opposite 
to  spin  axis 

CN  (“T*  rc^rJ  * 

See  Input  Data 

Hr  -qSdCwT 

/  \ 

Pitching  moment  about  aerodynamic 

Y  axis 

Cm  ‘  f  («r-  rc'r)  * 

See  Input  Data  (CM) 

r  -  d/2 

Radius 

S  *  ird2 

Reference  area 

Uu  "  *Sl  CMT 

Where  K  Is  Input  data 

L(u»)  ■  q  S  d  /  d  k,  (o 

KT 

Magnus  moment  roll  moment  about 
aerodynamic  X  axis  due  to  spin 

flg-qSd  /  d  XCj^q 

w 

Pitch  damping  moment  about  aero¬ 
dynamic  Y  axis. 

chq  ■  f  (®r) 

See  Input  Data 

Hp  -  q  s  d  -  Cmq^ 

Dlff  between  roll  and  pitch  damping 

C£p  •  f  (oT) 

Roll  damping  coefficient  about  aero¬ 
dynamic  X  axis.  See  Input  Data 

-  q  s  d  /d VC„ 

l2VfJ 

Spin  damping  moment 

cn  «  f  (<*t)  0 

"u 

Spin  damping  coefficient 

♦For  Ca,  Cn,  and  Cmt  provision  Is  made  for  choice  between  1  or  4  values  of  rc/r. 
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5.  INERTIAL  AND  THRUST  TERMS 

Wj  «  Up  +  wc  +  AWj  Total  Weight 

(1)  Before  Ignition  t  <_  tj 
m  ■  0 

rc/r  ■  0 

AW!  -  AWj,  AI$I  «  AI$I,  AItI  •  AIU 
T  -  0 

(2)  After  Ignition  tj  >  t 

;  -i  !!eo 

9  tB 

rc/r  -  ri/r 
AWj  ■  Aljj  ■  AI^j  ■  0 

T  ■  Thrust  ■  Constant 


*s  *  *sp  +  *sc  +  AIs1 
lt  "  *tP  +  !tc  +  Alt1 


Spin  Inertia 
Transverse  Inertia 
Spin  Inertia  of  Pyrotechnic 
Transverse  Inertia  of  Pyrotechnic 

Spin  Inertia  of  cup 
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Transverse  Inertia  of  cup 


Itc  ■  me  r2  1  +  4  (hpA)  *  (hp/y  *  r  ( 

r[  '  *  *  (s/-) 

«,i  ■  ^r’ 


Inertia  loss  at  Ignition 


Rate  of  change  of  Inertia 
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6.  ANGULAR  RATES 


(1)  Axes  Rates 


-  %  |  COS  i >A 


+  sin  1|»A  At  -  mq\-  (\ _ j 

L\  h  w  /  V1,  “/J 


Yaw  Rate 


Roll  Rf’3 


-  |  sin  4»a  -  cos  / *t  "  Aq 


-  Qj  sin  \l>n  -  cos 


Pitch  rate 


+  | cos  i <A  ♦  sin  (pA  /  *t  '  ^  [ 

\  !s  “  / 


P  *  P1  KP 
Q  -  Q,  KP 


•  /(coot) 


where  KP  Is  Input  data,  1  for  free>precess 

0  for  non-precess 


Aerodynamic  precession  rate 


(2)  Spin  Rate 


Is  «  *  u 


Magnus  precession  rate 


Spin  acceleration 


%  +  fi 


Spin  rate 
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7.  AXIS  TRANSFORMATIONS  USING  EULER  PARAMETERS 

EP  (1)  •  h  (-P  EP(2)  -  Q  EP(3)  -  R  EP(4)] 

EP  (2)  -  >s  I  P  EP(1 )  -  Q  EP(4)  +  R  EP(3)] 

EP  (3)  -  *s  I  P  EP(4)  +  Q  EP(1)  -  R  EP(3)] 

EP  (4)  -  H  [-P  EP(3)  +  Q  EP(2)  +  R  EP(1)] 

EP(1 )  -  EP(1)o  +  /ep«)  dt  1  *  1,  ...4 

(1)  Initialize 

EP(1 )  -  C(^°/2)  C(E°/2)  C(*°/2)  +  S(Ho/2)  S(Eo/2)  $(*o/2) 

o 

EP(2)C.  -  (c(Ho/2)  C(Eo/2)  S(*o/2)  -  S(Ho/2)  S(Eo/2)  C(*o/2) 

EP(3)0  -  [c(Ho/2)  S(Eo/2)  C(*o/2)  +  S(Ho/2)  C(Eo/2)  S(*o/2) 

Ep(4)c  -  |-C(Ho/2)  S(Eo/2)  S(*o/2)  +  S(Ho/2)  C(Eo/2)  C(*o/2)] 

where  C  *  cos,  S  ■  sin 

(2)  Normalize 

EP(1)  -  EP(1)  1  -  1,  ...  4 

(EP(J))2 

j^l 

(3)  Transformation 

A2  +  B2  -  C2  -  D2  2(BC  -  DA)  2(BD  +  CA) 

af j  -  2(BC  +  DA)  A2  -  B2  +  C2  -  D2  2(CD  -  BA) 

2(BD  -  CA)  2(CD  +  BA)  A2  -  B2  -  C2  +  D2 

m 

where  A  *  EP(1) 

B  -  EP(2 
C  «  EP(3 
D  •  EP(4) 

1  ”  1,2,3  denotes  Inertial  X,  Y,  Z  axes 
respectively. 

j  ■  1,2,3  denotes  moving  flare  X,  Y,  Z  axes 
respectively. 

Similarly  ajjfl  denotes  transformation  for  Aero  axes, 
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8.  AUXILIARY  PARAMETER  TRANSFORMATIONS 


L  ■  N  cos°L  -  A  sin**!- 
0  *  A  cosof  +  N  sln0^ 

CL  =  Cj«|  cos°^  -  Ca  s1n<*j- 
Cq  *  CA  cosoij.  +  CN  slno^ 

Y-t4n''fr) 

where  Uy  *  J  X2  -  Y2 

'  ■  (f ) 

H«.  -  tan"1  /  a13  \ 

V^) 


Lift 

Drag 


Lift  and  Drag 
Coefficients 


Flight  path  angle 


Flight  path  heading 


Azimuth  of  plane 
containing  spin  axis 


hA  *  tan_1  (  a21A  /  allA  ) 
*A  *  tan_1  (  a32A  /  a33A  ) 


where  a2^A  *  a2j  cos^  +  a22  s^n^A 

al  1 A  =  all  C«*A  +  a12  s1n*A 

a32A  *  a32  C0*h  "  a31  s1n^A 

a33A  =  a33 

a31A  s  a3i  cos^  +  a32  s1ni|»A 


Tilt  of  spin  axis 

Azimuth  of  Aero  axes 
Roll  of  Aero  axes 
Elevation  of  Aero  axes 
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\  ■  (  *32v/a33v) 

where  a32#  •  «32A 

a33v  *  a33A  cos  aT  "  *31 A  s*n  “T 

PA  *  P  cosily  +  Q  s1n^A 
Qa  ■  Q  cost|/A  -  P  s1ntpA 

R  «  g  *32A  +  PA  tar^j 

A  VT  cos<*j 

Pv  *  PA  C0SV  RA  Sin°T 

°v  *  Tt  (  ff  "  9  *33v  ) 

“  “  VQv 

rv  -  9  y \ 

®v  ■  pv  +  tan  Y  Jrv  cos  *y  +  Qy  sin  *v| 

Hy  -  secy  |  Rv  cos*v  +  Qv  s1nfy| 

Y  *  Qv  cos*v  -  Ry  s1n*y 


Roll  about  velocity 
vector 


Aero  Axis  rates 


Velocity  axis  rates 


Flight  path  rates 
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9.  INITIAL  CONDITIONS/INPUT  CONSTANTS 


DESCRIPTION 

Initial  velocity  -  Vq  (feet  per  second) 

Flight  path  angle  -  yq  (degrees) 

Flight  path  heading  -  HvQ  (degrees) 

Euler  angle  -  Hq  (degrees),  heading 

Euler  angle  -  Eq  (degrees),  elevation 

Euler  angle  -  *  (degrees),  roll 

Initial  spin  rate  -  u)Q  (rpm) 

Initial  altitude  -  H  (feet) 

o 

Aircraft  velocity  -  (feet  per  second) 
Diameter  of  flare  -  d  (inches) 
Thickness  of  flare  -  hp  (inches) 

Initial  weight  of  flare  -  WPQ  (pounds) 
Density  -  Op  (slugs  per  cubic  foot) 
Weight  of  case  -  (pounds) 

Weight  lost  at  ignition  -  AWj  (pounds) 
Burning  time  -  tg  (seconds) 

Ignition  time  -  t.  (seconds) 

Final  time  -  (seconds) 

Final  altitude  -  *  h^  (feet) 

Scale  factor  for  P  and  Q,  KP 
Scale  factor  for  magnus  moment,  KM 
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T 


10.  NOMENCLATURE 


HATH  SYMBOL 

FORTRAN  SYMBOL 

DEFINITION 

A 

AXIAL 

Axial  force  parallel  and 
opposite  to  aerodynamic 

X  axis.  (1b.) 

a<j 

A(3,  3) 

Euler  angle  matrix. 

cA 

CA 

Axial  force  coefficients. 

CMT 

CMT 

Pitching  moment  coefficients 

cN 

CN 

Normal  force  coefficient. 

d 

DI 

Diameter  cf  flare  (ft.) 

D 

DRAG 

Drag  (1b.) 

9 

G 

Gravity  (ft/sec.) 

h 

HA 

Altitude  of  flare  above 
sea  level,  (ft) 

hp 

HP 

Thickness  of  propellant  (ft.) 

H.E.* 

H.E.PHI 

Euler  angles  of  flare  axes 

_ J _ 1.L  _  1 _ a.  A _ _ -1 _ in 

azimuth,  elevation  and  roll 
attitude  taken  In  that  order 
of  rotation,  (deg.) 


Ha.Ea.*a 

Euler  angles  of  aerodynamic 
axes  azimuth  elevation  and  roll 
attitude  taken  In  that  order 
of  rotation  (deg.) 

hv*y  *v 

Euler  angles  of  velocity  axes 
azimuth  elevation  and  roll 
attitude  taken  In  that  order 
of  rotation  (deg.) 

Is 

AIS 

Moment  of  Inertia  about  spin 

Z  axis,  (slug  ft2) 

It 

AIT 

Moment  of  Inertia  about  trans¬ 
verse  X,  Y  axes,  (slug  ft2) 

is 

DIS 

Rate  of  change  1$  (slug  ft2/sec). 

• 

it 

DIT 

Rate  of  change  of  1^  (slug  ft2/sec) 

L 

ALIFT 

Lift  (lb.) 

M 

AMASS 

Total  mass  (slugs) 

M 

DM 

Mass  flow  rate  )slugs) 

MC 

CMASS 

Mass  of  the  cup  (slugs) 

Rp 

BARMP 

Cross  damping,  pitching  moment 
per  unit  roll  rate.  (ft. lb/rad/sec) 
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f 

V. 

* 

\ 

\ 

\ 

\ 

V 

\ 

\ 

\ 

m 

velocity  along  flare  X.  Y,  Z  axis 
(rad/sec). 

q 

QBAR 

Dynamic  pressure  (lb/ft4) 

a 

r 

R 

Radius  of  flare  (ft) 

rc/r,  rl/r 

RCOR,  RIOR 

Radius  ratio 

S 

S 

Reference  area,  for-  forces 
and  moments,  (ft4) 

t 

T 

Time  (sec). 

tB 

TB 

Burning  time  (sec). 

*1 

TI 

Ignition  time  (sec) 

u.v.w 

u.v.w 

The  component  of  toatl  velocity 

along  the  flare  X,  Y  and  2 
axes  respectively,  (ft/sec). 

u.  v,  w 

UOOT,  VDOT,  ZDOT 

The  derivatives  of  U,  V.  and 

W  respectively  (ft/sec4). 

VT 

VT 

Total  velocity  (ft/sec). 

wc 

WC 

Weight  of  cup  (lb) 

WT 

WT 

Total  weight  (lb). 

x,  Y,  Z 

X.  Y.  Z 

The  components  of  position 

A 

of  the  flare  along  Inertial 

X,  Y  and  Z  axes/  (ft). 

X,  Y,  Z 

XDOT ,  YDOT,  ZDOT 

The  components  of  total  velocity 
along  Inertial  X,  Y  and  Z 

I 

Pa 

PHOA 

axes  (ft/sec) 

Density  of  air  (slug/ft3) 

Pp 

PHOP 

Density  of  pyrotechnic  (slug/ft3). 

aT 

ALPHAT 

Angle  of  attack  (deg). 

h 

PHIA 

Yaw  of  aerodynamic  plane  with 
respect  to  X  flare  axis.  (deg). 

k,  nn*"1 

SIM 

OMEGAM 

Magnus  precession  rate. 

% 

OMEGAT 

Aerodynamic  precession  rate. 

to 

OMEGA 

Spin  rate  (rpm) 

AWj 

OWI 

Weight  lost  at  Ignition  (1b). 

AISI 

DISI 

Spin  Inertia  lost  at  Ignition. 

AItI 

DITI 

Transverse  Inertia  lost  at 

* 

Ignition. 
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| 

Mr  AMT 

N  AN 

P.Q.R  P.Q.R 

</ 

c4j 

C-Cp»  CMq 

"71 

Cn 

cnu 

EP(1) 


Total  static  pitching  moment 

The  normal  force  parallel  to 
the  flare  Z  axis  (lb.) 

Component  of  flare  axes  angular 

Rolling  moment  about  aerodynamic  X 
X  axis 

Rolling  moment  coefficient 

Magnus  moment  coefficient 

Roll  and  pitch  damping  moment 
coefficients  respectively 

Yawing  (spin)  moment  about 
Spin  Axis 

Yawing  Spin  moment  coefficient 

Spin  damping  coefficient 

Euler  parameters  used  to  generate 
axis  transformation  matrices 
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APPENDIX  C 

EQUATIONS  FOR  STABILITY /PRECESSION  ANALYSIS 

Transforming  the  flare  axis  equations  of  Appendix  B  to  the  velocity  axes 
results  in  the  following  equations  in  terms  of  velocity  and  aerodynamic  axis 
parameters  which  were  then  linearized  for  stability /precession  analytical 
solutions. 

1.  ANGLE  OF  ATTACK 

aT  =  QA  "  QV* 

2.  ROLL  ANGLE  OF  LIFT  VECTOR 

t  V  =  Py  +  tan  Y  (Ry  cos  i  +  Qv  sin  *). 

3.  AZIMUTH  OF  VELOCITY  VECTOR 

Hy  =  sec  Y  ^  sin  $  . 

4.  ELEVATION  OF  VELOCITY  VECTOR  (FLIGHT  PATH  ANGLE) 

V  =  Qv  cos  I  -  Rv  sin  *y  •  v  [(k)  cos  V  '  e  COS  Y]’ 

a.  Roll  Rate  of  Velocity  Vector  Axes 

Pv  =  PA  cos  aT  +  RA  sin  a^. . 

b.  Yaw  Rate  of  Velocity  Vector  Axes 

Ry  =  RA  cos  aT  -  PA  sin  aT 
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c.  Pitch  Rate  of  Velocity  Vector  Axes 


%  =  V  [k  _ 


g  cosy  cos  t 


Velocity 


v  =  "  m  "  gsin  Y* 


d.  Roll  Rate  of  Aerodynamic  Axes 


pa  -  V*  +  r  ’  '  ra1  qa  -  aop>' 


Pitch  Rate  of  Aerodynamic  Axes 


QA  "  MQQA 


‘ft-  *  ‘  Ra] 


*A  *  "TJT  (a  )  =  Ma  (aT  ‘  aOP)< 


e.  Yaw  Rate  of  Aerodynamic  Axes 


[cos  y  sin  tv 


cos  a 


+  PA  tan  aT. 


where 


<Xqp  =  angle  of  attack  for  zero  pitching  moment 
a  OM  =  aOP  =  anS^e  attac^  for  zero  Magnus  moment. 

M_  =  T  C.. 
a  \  Ma 

_  erf  H 

Mq  =  q  -j —  CM  "W  =  Pitch  Damping  Derivative, 

t  Q 
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Magnus  Derivative. 


cL„„  •  T  It  c 


T  "V57 


qj  j 

=  q  clp  W  =  Roll  Damping  Derivative. 


aCM  3Cl 

and  where  CM  ,  Cj  =  — ,  — -  -  ,  respectively,  the  slopes  of  the 
Ma  iuja  ®  ° 


pitching  and  Magnus  moment  coefficient  curves. 


—  = C.  =  -3-^-  CT  a_  =  Lift  acceleration, 
m  m  L  m  T 


CL  =  CN  "  CA  =  ~ST 
a  a 


=  Slope  of  the  lift  coefficient. 


Cn  =  =  Slope  of  the  normal  force  coefficient. 

oi 
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APPENDIX  D 

TYPICAL  CONFIGURATION  DATA  INPUT  FOR 
TRAJECTORY  SIMULATION 


TYPICAL  GEOMETRIC  AND  INERTIA  DATA 


.  Diameter,  d:  As  applicable 

.  Thickness,  hp:  Flares,  2  in,  others  as  applicable 

.  Pyrotechnic  Density,  pp:  2.  86  slugs  per  cubic  feet 

.  Case  Weight,  Wc?  Live,  0.  5  pound;  inert,  0. 

.  Initial  Pyrotechnic  Weight,  Wp~:  Live,  Actual /intended;  Inert, 
Total  Weight  (used  3  pounds  fornB  inch  flare,  12  pounds  for  12  inch 
flare,  etc.  ). 

.  Burning  Time,  t^:  Live,  as  desired;  Inert,  any  finite  constant  >0. 

.  Ignition  Delay  Time,  t.:  Live,  as  desired;  Inert,  any  arbitrary 

number  >  flight  time. 1 

.  Weight  Dropped  at  Ignition,  AWj?  As  required,  usually  used  0. 

.  Thrust,  T ;  As  desired. 
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STATIC  AERODYNAMIC  COEFFICIENTS  (Concluded) 


